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The Story of Shimmy 


WILLIAM J. MORELAND i 
Wright Lir Development Center 


SUMMAR\ 


A theory of landing-gear shimmy is developed in which the basic 
assumptions involving the part played by the tire, the air frame, 
and the gear proper have been individually submitted to direct 
test 

Analysis and analog studies of the theory are used to establish 
a design procedure for obtaining a stable gear with a minimum of 
artificial damping 
been tested on a model 


The predictions of the theory have 


whose parameters can be varied to simulate nosegear structures 
of light high-speed fighters, as well as heavy bombers. A com 
plete description of the model and a summary of the test results 
are presented 

Full-seale field test results are given for instrumented fighter, 


bomber, and cargo-type aircraft 


INTRODUCTION 


ames IS THE term used to describe the self-induced 
swiveling of the airplane. Its 
appearance is unmistakable, and is characterized by 
the rapid build-up of a high-frequency oscillation ac 
companied by a shudder which progresses along the air- 
In some in- 


nosegear of an 


craft from the nose to the tail surface. 
stances the action is too rapid to permit the pilot to 
take the necessary corrective measures and the destruc- 
tion of part of the plane may result. 

The problem has been studied intensively for more 
than 20 years, yet, even since 1950, millions of dollars 
have been spent because of shimmy troubles. There ts 
an interesting paradox about shimmy. It is one of the 
most widely known of dynamic phenomena. 
aside from the esthetic meaning of the term it is almost 


Even 


Presented at the Acroelasticity Session, Twenty-Second Annual 
Meeting, IAS, January 25-29, 1954. 

* The following individuals gave valuable assistance in the 
Sanborn and 


several phases of the program: computations —T. G 
Clarke, 


W.C. Bennewitz; analogue studies——-W. L. Cyrus, J. G 
and R. T. Harnett; model tests—N. W. Carlson, P. E. Gies, and 
kK. A. Erfurth; field tests—L. T. Waterman and D. W. Weller, 
Fairchild Engine and Aircraft Company C. Minch 

t Chief, Mechanics Research Branch, Acronautical Research 


and F 


Laboratory, 


a household word because of its common occurrence in 
automobiles, trailers and tea carts. But as recently as 
1952, three mutually contradictory theories on this 
subject have appeared. True, means for controlling 
shimmy are available but they are restricted in their 
application and provide little toward an understanding 
of its cause. The explanation for this anomalous situ 
ation is not that insufficient attention has been focused 
on the problem. Some 300 papers*® bearing directly or 
indirectly on this subject have been published since 
and important contributions have been made by 
Broulhiét; by the 


1922 


the French workers, Lavaud and 
Germans, von Schlippe, Dietrich, Dietz, Harling and 
Fromm; and by the British and Americans, Lanchester, 
Temple, Bull, and Evans. The primary purpose of 
this paper is to describe the program which was pursued 
at the Wright Air Development Center to study the 
problem of shimmy. In addition, some attention will 
be given to learning the cause of the inadequacy of ear 
The program was initiated in 1951, follow 
The 


lier efforts. 
ing a series of accidents traceable to landing gears. 
work was planned to cover four distinct areas: 

(1) Development of a theory of shimmy. 

(2) Electronic computer studies of the theory. 
(3) Experimental research on laboratory model. 

1) Full-scale field testing of fighter, bomber, and 
cargo-type aircraft. 

SHIMMY 


DEVELOPMENT OF THE THEORY OF 


Earlier Efforts 

In 1925, Broulhiét in France, published a paper in 
which he succeeded in showing that the mechanical 
properties of the pneumatic tire could account for a 
sustained oscillation 1n a swivelling structure. Follow- 
ing this work the notion became widespread that the 
problem of achieving stability in a castering wheel 
depended in great measure on the properties of the 
This viewpoint is still adhered to by the majority 
The philosophy 


tire. 
of analytical workers in the subject. 
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of this school is expressed in the paper by Bourcier de 
Carbon (translated from the French and published 
1952 by NACA, TM 1337). It will be helpful to 
briefly review the main points of this paper in order to 
obtain a clearer understanding of the reasons for the 
radically different conclusions obtained from the WADC 
approach. The mathematical model, Fig. 1, chosen 
by de Carbon to represent the nosegear system, con- 
sists of a rubber-tired wheel whose axle is carried in a 
fork which in turn is permitted to swivel in a fixed ver- 
tical pivot after the fashion of an ordinary furniture 
caster. The significant property of this simple system 
is based on the observation that a rolling wheel which 
has no force acting normal to its plane will describe a 
straight line over the ground parallel to its plane. 
But, if a steady force is applied normal to the wheel, 
it will drift sidewise while rolling and thus trace out a 
line on the ground which makes an angle (called the 
drift angle) with the wheel plane. This angle is di- 
rectly proportional to the applied force. The action 
is still one of pure rolling—i.e., does not involve skid- 
A theoretical explanation of the phenomenon is 
Employing this experimental 


ding. 
given by von Schlippe.* 
law, the simple two-degree-of-freedom system is readily 
described by a third order differential equation which 
yields the stability criterion” 


ap > Tt (1) 
where 
a = trail distance, in. 
D = drift angle divided by the force causing the 
drift, rad./Ib. 
T = the lateral deflection of the tire per pound of 


force. (Found by measuring the sidewise 
deflection of the nonrolling wheel under the 


action of a side force) 


Eq. (1) states that if the trail a is made sufficiently 
large so that the product of trail times drift coefficient 
exceeds the lateral compressibility of the tire, then the 
This 


law is, of course, rigorously true for the system under 


wheel will not experience a self-excited vibration. 


discussion. 
qualitatively, to an aircraft structure. 
this will be discussed in detail in the remainder of the 


But, this simple rule does not apply, even 
The reasons for 
article. However, the most convincing explanation is 
quite simple. 
to have a rigid tire for which 7 is zero while D has a 


Suppose the model, Fig. 1, is assumed 
value approaching zero. The system would then be 
stable according to Eq. (1). 
assuming that the pivot is fixed laterally, it is per- 


However, if instead of 


mitted to have a lateral degree of freedom against a 
restoring force whose coefficient of elasticity is Aj, 
then, as shown in reference 1, page 9, this undamped 
system is inherently unstable no matter how large the 
lateral rigidy (K,) may be. The assumption that the 
rigidity of the pivot is sufficiently great to permit 
describing it as “‘fixed’’ is tempting, but untenable. 


RARER ONAUTICAL 


SCIENCES DECEMBER, 1954 

Even if we permit ourselves the luxury of assuming a 
heavily reinforced structure, the “‘fixed’’ pivot theory 
will not apply. 

Fig. 2 may be used to illustrate this point. Sup 
pose it is required to determine the stiffness of the 
spring, one end of which is attached to the finite mass 
mt. If point A is moved slowly 1 in. toward the right, 
no resistance is felt at point A as m7 slides over the fric 
tionless support. Under these circumstances, it must 
be concluded that the rigidity of the system is zero since 
a l-in. deflection without This 
result is quite independent of the A of the spring itself. 


was obtained force. 
If point A is caused to move with simple harmonic 
motion of l-in. amplitude, it will be found that, at low 
frequencies, when A has moved | in. to the right the 
mass #7 will actually have traveled more than one inch 
in the same direction. The spring will be stretched and 
in tension, and the force at A will be directed toward the 
left. 
of requiring a push, actually demands a pull, it must 
be concluded that the rigidity of the system is negative. 
Stiffening the spring itself will not change this result. 
Next, if the frequency of oscillation of point A is gradu- 
ally increased, the negative rigidity will increase until 
the natural frequency of the system is reached. Be- 
yond this frequency it will be found that when point A 
is at the extreme right of its travel, the mass mm will be 


Hence, since a l-in. motion to the right, instead 


at the extreme left of its excursion with a resulting high- 
compression in the spring and a positive value for the 
rigidity of the system. As the frequency is further in- 
creased, the mass will approach a condition of rest and 
under these circumstances the rigidity of the system 
will approach the A of the spring itself. 

If the above experiment is repeated while another 
spring and mass are attached on the right-hand side of 
m, the variations in the apparent rigidity of the system 
will be found even more significant. In fact, a fre- 
quency will always be found at which the motion of 
will be exactly the same as that of A. Under these 
circumstances, the rigidity of the system is actually 
zero regardless of the stiffness A of the spring itself. 
This same conclusion applies to the airplane. Even 
when the strut itself is made quite rigid the effective 
rigidity of the system may be small. 

Fig. 3 shows the relationship between the effective 
rigidity and frequency for this idealized model. The 
relevance of this discussion to the problem of deter- 
mining the ‘‘true rigidity’’ of the nosegear structure is 
it cannot be found by static deflection of the 
strut. As shown by instrumented taxi tests on a shim- 
mying airplane, stiffening the strut cannot suppress the 
Hence, a theory of shimmy 


obvious 


lateral motion of the pivot. 
based on the elastic properties of the tire alone is in- 
adequate. A nosegear structure designed on the basis 
of the rule aD > T would, if not otherwise stabilized, 
be destroyed on its first taxi run. But this is not all. 
An even more significant degree of freedom than the 
lateral motion of the pivot is that due to the finite 





torsit 
shim 
of lir 
oil tl 
to Cc 
size i 
syste 
finite 
this 
that 
ticity 
syste 
by tl 
siona 
ertia, 
latin: 
Tt 
stabi 
fairly 
part 
Ti 
freed 
expe 
term 
such 
made 
of th 
fight 
show 
by E 
woul 
possi 
Ti 
in th 
of th 


wher 
I 
i 


Since 


escay 


In ed 
is pr 
porte 
later 
assul 
reimc 
show 
of th 
Se = 


torsi 


ing a 


heory 


lass 
‘ight, 
fric 
must 
since 
This 
tself. 
lonic 
low 
the 
inch 
and 
| the 
tead 
ust 
ive. 
sult. 
idu- 
ntil 
Be- 
tA 
| be 
gh- 
the 
in- 
ind 


ell 


her 
of 


em 





THE Stoxgy 


torsional elasticity of the damper system. During a 
shimmy, the wheel will swivel in its pivot and by means 
of linkage this motion forces a damper piston to drive 
oil through an orifice which offers sufficient resistance 
to cause the oscillation to attenuate. If the orifice 
size is reduced, the damping action is increased and the 
system becomes still more stable. However, if the 
finite elasticity of the linkage is taken into account, 
this conclusion may be reversed. It can be shown! 
that below a certain critical value of the linkage elas- 
ticity, no orifice size can be found that will stabilize the 
system. This 
by the nondimensional number composed of the tor- 
sional and lateral rigidities, the wheel moment of in- 
The method for calcu- 


critical elastic coefficient is disclosed 


ertia, and the strut weight. 
lating this number will be given later. 

This does not imply that the influence of the tire on 
stability may not be important. Only by analysis of 
fairly complete models of the whole structure can the 
part played by the tire be properly evaluated. 

This is the key to the long controversy; no degree of 
freedom can be neglected until numerical analysis and 
experiment have demonstrated the part it plays in de- 
termining the system stability. Over two thousand 
such analyses of various types of aircraft have been 
made at the Wright Air Development Center. In one 
of these, using experimentally determined data from a 
fighter-type aircraft and published‘ tire data, it was 
shown that in order for the stability criterion given 
by Eq. (1) to apply, the lateral rigidity of the pivot 
would have to be increased threefold—a physical im- 
possibility. 

The author next examines the effect of transient drift 
in the tire and from this study obtains a modified form 
of the original stability criterion: 


(URI? — a) (aD — T)>0 (2) 
where 
{/ = moment of inertia of the oscillating system 
about the pivot 
V = velocity of the aircraft 


Since stability exists if both terms are positive, the in- 


escapable conclusions to be drawn from Eq. (2) are: 


(1) Increasing the system inertia, /, is stabilizing. 
(2) Increasing the taxi speed, J’, is stabilizing. 
(3) Sufficiently great trail, a, is destabilizing. 


In each of the above conditions the opposite conclusion 
is predicted by WADC Theory and the latter is sup- 
ported by both model and full-scale tests to be described 
The erroneous conclusions are due in part to the 
With this restriction 


later. 
assumption of a “‘fixed”’ pivot. 
removed, the above rules are actually reversed as 
shown in reference 1. (Because of a misinterpretation 
of the “law of turning,” 
Se = Dis obtained, where S = 6/1/ is the coefficient of 
torsional elasticity of the tire, and « = J Fis the pneu- 


reference 2, p. 23, the relation 
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matic trail. 
equation, some experimental data were substituted in 


In what appeared to be a check of this 


what closer examination reveals to be a mere identity, 
viz., (6/17) (17 F) = (6 F). Of course, substitution of 
any arbitrary numbers will give an ~check.”’ 
Published data‘ giving independent measurements of 
S, «, and D do not verify the given equation.) 


exact 


Comments on the Wright Air Development Center Theory 
of Shimmy 
Details of the mathematical development of the 
WADC theory of shimmy! will not be repeated here. 
The basic assumptions employed, the experimental 
evidence in support of these, and the significant con- 
clusions are discussed. 
Degrees of Freedom.—The degrees of freedom of the 
mathematical model used in the development of the 
theory are (1) the lateral displacement of the wheel 
while the tire footprint remains fixed on the ground. 
This is due to the lateral flexing of the tire; (2) the 
swiveling motion of the wheel about the strut axis; (8 
the lateral displacement of the pivot or strut against the 
elastic restraint of the strut and airframe; (4) the rela- 
tive motion between the swiveling of the wheel and the 
movement of the damper piston. This relative mo- 
tion is due to the elasticity of the damper linkage; (5) 
the relative motion between the strut and the airframe. 
(The airframe itself is in turn represented by a lumped- 
parameter system of one or more degrees of freedom. 
Taxi tests on instrumented bomber and cargo-type 
aircraft have been conducted to obtain quantitative 
data on the relative magnitudes of the motions in the 
various degrees of freedom. Velocity and acceleration 
type instruments were placed on the nosewheel, along 
the fuselage, in the wings and stabilizers. Oscillograph 
records of the instrument readings were taken for vari 
ous taxi speeds and for shimmy conditions which varied 
from mild to severe. These records disclose the fre- 
quency, amplitude, and phase of the several parts of the 
airframe. Instruments located as far as 50 feet from 
the nosegear recorded vigorous vibrations of shimmy 
frequency usually within one second of the onset of the 
disturbance at the gear. The lateral motion of the 
pivot (or projected strut) is always present and fre- 
quently exceeds the lateral movement of the wheel. 
Lateral Rigidity. 
rigidity of the strut is frequency-sensitive. 


As previously indicated the lateral 
Experi- 
mentally, this quantity is obtained by applying a har- 
monic force at the axle in a direction perpendicular to 
the axis of the airplane and, with the wheel removed, 
record the lateral deflection as a function of the fre- 
quency of the vibration. A fighter-type aircraft under 
test at WADC and Cornell Aeronautical Laboratory 
showed a variation in lateral rigidity from 500 Ibs. /in. to 
12,000 Ibs. ‘in.—a ratio of 24 to 1—as the frequency was 
varied from below to above normal shimmy ranges. 
Hence, the lateral rigidity cannot be known until the 
shimmy frequency is given, and the shimmy frequency 
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is itself dependent on the lateral rigidity. A graphical 
solution of these simultaneous functions is described in 
Reference 1. 

How the Tire Influences Shimmy. 
controversial subjects in the problem of shimmy con- 
cerns the part played by the pneumatic tire. The ap- 
proach taken in the WADC report is based on assump- 
tions which individually can be submitted to direct 
A wheel rolling over 


One of the most 


test. The theory is as follows: 
the ground with a forward velocity V (Fig. 4) will, when 
acted upon by a force F from the ground, travel along 
the path |’ making the yaw angle y with its plane. 
Thus the wheel center has a velocity |, = Vy normal 
to its plane and directed against the ground force. But 
the force F does not instantly create the yaw angle y. 
A measurable time delay follows the application of the 
force before a steady yaw angle is established and this 
condition can be described by the equation 


¥+C.v¥ = CF (3) 
where 
C = the coefficient of yaw, rad. Ib. 
C, = the time constant, sec. 


Due to the lateral compressibility of the tire, the center 
of the wheel will be displaced relative to the ground con- 
tact point whenever a force F acts between the ground 
This force is transmitted through the rubber 
The accompanying 


and tire. 
from the ground to the wheel. 
strain in the rubber is 


A = pF (4) 
where 


p = the lateral coefficient of compressibility of the 
tire, in. /Ib. 


Hence, the total normal velocity of the wheel due to the 
combined effects of yaw and tire stretch is Vy + A, 
A third source of normal velocity exists when the wheel 
is trailing the pivot at an angle @ with the direction of 
forward velocity of the airplane. This is most easily 
pictured by assuming the pivot at rest and the ground 
moving under it toward the left. As shown in Fig. 5, 
this creates a velocity normal to the wheel plane equal 
to — V0. 

Combining the three effects, the total normal velocity 


becomes 
V, = Vw¥+A— Vo 
and the angular velocity of the wheel about the pivot 
is thus 
v ly+ A— Ve R 
= §)= (>) 
L L 


This equation defines the coupling between the two de- 
grees of freedom, @ and A, and may be combined with 
the dynamic equations of the system, including the 
lateral motion of the pivot, as described in reference 
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|. In this way (as discussed under computer studies) 
a true perspective of the part played by the tire in the 
dynamics of the overall system is obtained. (Both 
the present theory and that of de Carbon employ the 
same basic law of tire yaw. Hence, the former must 
contain the stability criterion given by de Carbon. 
This can be seen by combining Eqs. (B1) and (5) of 
reference | and omitting the term that describes the 
lateral motion of the pivot.) In the derivation of Eq. 
(5) above, the critical assumption is that expressed by 
Eq. (3) and this may be subjected to direct test in the 
following way: Place the wheel and pivot assembly in 
a fixture that will permit free swivelling of the wheel 
while the latter is pressed against a rotating drum. 
A force applied to the axle normal to the wheel plane 
will cause the trail arm to swing about the pivot and 
assume a fixed angular position with the direction of 
This fixed angle is the yaw created by the 
The coefficient of yaw C is the ratio of 


free trail. 
normal force. 
the angular displacement to the applied force. In 
addition, if an oscillograph record is made of the angular 
motion following the application of a step input force, 
it will be found that the wheel will oscillate about its 
equilibrium position of steady state yaw with a decay 
The logarithmic decrement of this 
Thus, 


ing amplitude. 
motion defines the time constant C; of the tire. 
a single simple setup provides the means for determining 
C and C, simultaneously. Tests conducted at the 
Wright Air Development Center have fixed this latter 
quantity at about 0.01 sec. for the tire used on the 
model nosegear. The significant point to observe here 
is that direct dynamic measurements involving no un- 
substantiated assumptions have established the co- 
efficients of Eq. (3). It is important that a theory of 
shimmy provide for the concept of the time constant C;. 
This is apparent when it is recalled that for shimmy fre- 
quencies encountered in practice, the yaw angle is built 
up from zero to its maximum value in about !/¢ of a 
second. For this reason, a theory of shimmy which 
assumes that the full steady state yaw is associated 
with a given side force will be considerably in error. 
Replacing the airframe by a two-mass lumped 
parameter system, the complete mathematical model 
will contain six degrees of freedom. From the set of 
simultaneous equations that express the dynamics of 
the parts of the model a single tenth-order differential 
equation for the complete system is obtained. De- 
tails of the process differ according to the coordinate 
system used to describe the motion of the strut. If the 
strut is assumed to oscillate laterally about a fixed point, 
polar coordinates apply; but if the point about which 
the oscillation occurs is taken infinitely far above the 
wheel axis, rectangular coordinates may be used. The 
choice is not arbitrary and its significance will be 
brought out when discussing the theory of the experi- 
mental model. Rectangular coordinates are employed 
in presenting the theory in reference |, while polar co- 
ordinates are used to describe the experimental model. 
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THE STRUT MOVES FORWARD IN 
A STRAIGHT LINE WITH NO LATERAL 
VIBRATION (FIXED PIVOT). THE 
WHEEL SWIVELS AROUND THE 
STRUT BY VIRTUE OF THE ELAS- 
TIC PROPERTIES OF THE TIRE. 
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The problem now at hand is to determine whether re 
gions of stability for the nosegear can be mapped from 
the dynamic equations. The problem is considerably 
simplified if the system is assumed divided at the junc- 
tion of the self-excited portion with the passive airframe. 
Analysis of the two parts may be carried out separately 
and a simultaneous solution obtained by an impedance 
match at the junction. In this way the search for sta- 
bility criteria is directed to the portion of the system 
which contains the most significant parameters. In 
addition, since normal ranges of tire parameters have 
been found to have remarkably little influence on the 
damping required for stability, the tire may at first be 
taken asrigid. This simplified system can now be com- 
pletely described by using only four nondimensional 
parameters.? As shown below, the ratio of two of these 
numbers, called the inertia ratio, is a useful index of 
stability. -the 
parameter which plays so prominent a part in tire mode 


This number does not contain the trail 
theories. The principal conclusions from the analysis, 
applicable to the design of a stable nosegear structure, 
are repeated here. Let 
QO V2V 12/LV KK, 
R = Cy 2 VAT; 
= nondimensional strut mass = NL?/T, 
b/s H/ R 


nondimensional ‘‘inertia ratio” 


nondimensional velocity = 


nondimensional damping = 
Bs 

T = nondimensional trail = 
T?/S = 

C, = coefficient of viscous damping, Ib. in. sec. 
I; = moment of inertia of nosewheel about a di- 
ameter plus moment of inertia of spindle 
about its axis, Ib. in. sec.” 


K, = lateral rigidity of strut, lb./in. 

aa = torsional elastic coefficient of the damper 
linkage (obtained by applying a couple in 
a horizontal plane to the wheel while the 
damper piston is locked), Ib. in. ‘rad. 

i = tral, im. 

N = equivalent mass of strut (the lateral moment 
of inertia of the strut about its nodal 
point divided by the square of the dis- 
tance between this point and the ground), 
Ib. sec.?/ in. 

V = ground speed, in./sec. 


(1) The complete stability criterion for the fourth 
order rigid tire system is 


QS(1 — T?/S)R? + [1 + S(2 — T?/S)]R?2 + 
(1 + S)/Q — QS/4JR- (1+ S)/4>0 = (6) 


Examination of this inequality will reveal all of the 
essential phenomena pertaining to the dynamics of the 
fourth order system, and, as shown below, will give the 
damping required to stabilize the complete seventh order 
system for all normal parameter ranges. It should be 
particularly noted that only three numbers, Q, S, and 7 
are required for analysis of any given airplane. Fami- 
lies of curves can be constructed! from this expression 
to portray the effects of variation of any one of 9 system 
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parameters. However, simple inspection of Eq. (6 


yields significant information. For example, it is 
apparent from the first term of Eq. (6) that if 7° S§ 
exceeds unity, then sufficiently large values of damping 
(R) will cause the system to become unstable. In 
effect, the main gear has infinite damping since it is 
not permitted to swivel, and in this case the system 
S (the 


Hence, the condition that must 


cannot tolerate any advance of 7° “inertia ra- 
tio’’) beyond unity. 
be satisfied to insure stability without the danger of 


overdamping is that the inertia ratio 


7? I ma ¢ 
S NK, ” 
For example, measurements of a fighter-type jet nose- 
gear gave /,, = 3; K, = 0.5 X 10°: N = 0.033. Hence, 
I? Ih, a Ky 
S NE, 00383 X 035 x 10° 


or AK; must not exceed 5,500 Ibs. in. when measured at 
the shimmy frequency of the nosegear. If A, is larger 
than this amount, the first term of Eq. (6) becomes 
negative and, for large values of R, the left side of the 
inequality is negative indicating instability. 

(2) The damping required to maintain neutral sta 


bility when 7*/.S = 1 is found from Eq. (6): 
QS | QS i- 
“oe _ + | 
$H(1 + S) @) \ ca oe) O 
R = : | (GA 


») 


Eq. (GA) applies quite closely for all values of 7°) S 
from 1 to zero. Ample attenuation of wheel disturb 
ance is provided when C, is given a value between two 
and two and one-half times that found from Eq. (6A). 
It is significant to observe the values of C, required 
for stability when extreme values of the parameters 


are substituted in the above equations. For example, 


when 
V = 0, C,=0 
: ; Kita : 
V = o, C, = — [from Eq. (6) ] 


js 7 


Hence, C, must always be made greater as J” increases, 


yet will remain finite if 7°/.S < 1. (Gyroscopic effects 
reduce the required damping at high speeds but this 
favorable result is approximately cancelled by tire yaw 
and flexibility.) 


=e C, = VK1;z, ie., stability is possible 
with zero trail 

L = 2, C,=0 

K, = ™, C,= VI;/L 

I; = 0, C,.= 0 

Iz; = @%, C, = « 


It is apparent that J; is the most important factor in 


determining C;. It should be noted that increasing /; 
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is here predicted to have the opposite effect to that 
found on the basis of tire effects alone. As will be 
shown later, both model and full-scale tests support 
the conclusion that increasing the wheel moment of 
inertia is destabilizing. Returning to Eq. (6), it can 
be seen that for 7°*/.S < 


tive, the third may be positive or negative, while the 


1, the first two terms are posi 
fourth 1s negative. ‘That is, there 1s only one change in 
sign and hence no more than one positive root. All 
values of R greater than this root are stabilizing. But 
if 7*/S > 1, there can be two changes in sign indicating 
a possibility of two values of R for neutral stability be- 
tween which the system is stable. Hence, stability 
may exist for values of the inertia ratio greater than 
unity. But for this to be a practical solution the two 
roots must be spread sufficiently far apart so that inter- 
mediate values of R will provide sufficient attenuation 
and also so that the risk of accidental overdamping will 
be avoided. Analog studies of the system have shown 
that if the larger root R is taken as 2.5, satisfactory sta 
bility can be obtained over the complete range of system 
parameters. If this value of R is substituted in Eq. 
(6) and the inequality solved for 7° .S, it can be shown 
that when 7° Sis less than (1 + S$)/2.5 QS the stability 
will be satisfactory. Substituting for the dimensionless 
terms gives 
; [1 + (NL?/Iz)] K, 
A, < (S 


SVLV 1; 

K, should be computed from Eq. (7) and Eq. (S) and 
the larger value retained. Further analysis ef Eq. (6 
shows that when the inertia ratio is continuously in- 
creased beyond unity the two real roots approach each 
other and finally coalesce and beyond this point sta- 
bility is impossible regardless of the amount of damp- 
ing employed. 

(3) ‘‘Middle”’ values of trail are destabilizing. This 
conclusion was reached by a number of experimenters 
(Tonnies and Gollung, German View of Tricycle Under- 


carriage, Aircr. Engrg., May, 1947). Also (Seifert, 
Stability in Taxiing, ME-GE 1940) and (Fairchild 


C-119 Taxi Tests 1953). No explanation had been 
given for these observations, but as shown in reference 
1, page 32, this phenomenon has a simple theoretical ex- 
planation due directly to an excessive value of the inertia 
ratio. 

(4) For the values of strut mass and trail commonly 
encountered, stability is impossible if the inertia ratio 
exceeds a calculable upper limit. This conclusion is 
valid regardless of the amount of viscous damping em- 
ployed. (See Fig. 7a, reference 1). 

(5) The above conclusions apply strictly to the case 
of a rigid tire with negligible gyroscopic effects, but re- 
main valid when normal tire parameters‘ are inserted 
in the general equations of the system (see Appendix B, 
reference 1). Pneumatic trail—i.e., the additional trail 
induced by the shifting aft of the center of pressure in 
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the tire footprint, has little eect on the damping re- 
quirement. 

6) No theory of shimmy can be complete. Factors 
such as nonlinearity of parameters, backlash, fore-and- 
aft strut motion and extreme tire parameters play a 
part in the phenomenon. Results of studies to date 
at the Wright Air Development Center indicate that 
conditions these effects are not pre 


under normal 


dominant. 


Electronic Computer Studies 


The purpose of these studies was to examine the 
effects on stability due to the pneumatic tire. A system 
of four degrees of freedom containing 15 parameters 
was selected. Twelve of the 15 parameters were ob 
tained by direct measurements from a fighter-type air 
craft. The tire characteristics C and p were obtained 
from published data* and C; was determined from model 
tests at the Wright Air Development Center. The 
following procedure was used in setting the problem 
on the machine: All of the tire parameters were set 
equal to zero (simulating a rigid tire) and the wheel was 
given an initial displacement corresponding to a lateral 
shock. 
by trial until a neutrally stable shimmy was obtained. 


The damping coefficient C, was then adjusted 


The same procedure was repeated with different taxi 
Next, 
combination 


speeds and a graph of C, vs. speed was plotted. 
with the same structural parameters, a 
of tire parameters, C, C), p, were selected so that accord- 
ing to the stability rule obtained from the fixed pivot 
theory,” the system should be strongly unstable—.e., 
LC < . 
computer and data were again obtained to plot C, for 
The results are shown 


The complete system was then placed in the 


neutral stability vs. taxi speed. 
in Fig. 6. 

It is to be particularly noted that at every point the 
system that included the elastic tire is more stable, 
i.e., a greater region of stability is enclosed by the graph, 
despite the fact that the tire system by itself is un- 
stable. It might be reasoned that although the tire 
system by itself is unstable, the presence of C, serves 
to hold it stable. But this does not explain the fact 
that the combination of the two unstable systems re- 
quired less damping than that demanded by the original 
Fig. 6 is typical of the results obtained over 
Only for abnor- 


system. 
a wide range of system parameters. 
mally large values of tire parameters does the desta 
bilizing influence of the rubber become comparable to 
that of the structural parameters. 

A second series of studies was made using a com- 
bination of tire parameters that would be stable on the 
i.e., LC > p (or, in de Carbon’s 
In every case, this system was 


fixed pivot theory? 
notation, aD > T). 
unstable for normal ranges of structural parameters 
unless artificial damping was employed. Moreover, 
when the lateral rigidity of the strut was increased from 
its normal value of 7 X 10° Ibs./in. to 10 K 7 X 10° 


Ibs. in., the system was still unstable, again demon- 
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strating the impracticality of attempting to base a 
stability criterion on the mechanical properties of the 
tire alone. 

Fig. 7A is an oscillograph record of the lateral move- 
ment of the wheel during a neutrally stable shimmy. 
The system is composed of a laterally soft (unstable) 
Following the 
low 
the 
the 
the 
strut mode, the tire component is strongly attenuated. 


tire and an abnormally rigid strut. 
initial displacement, about three oscillations of 
frequency (the tire mode) can be observed, while 
With 
damping increased (Fig. 7B) to slightly attenuate 


high-frequency strut mode is sustained. 


In a sense, these diagrams visibly portray the grounds 
on which the opinions of the two schools of thought are 
based. Both modes of motion that due to the tire 
and that due to the structure 


In this particular case where 


are ever present and in- 
extricably interwoven. 
sufficient damping is provided to control the tire mode, 
it is obviously not sufficient to give stability to the en- 
tire structure. This again demonstrates that generali- 
zations from particular cases may be misleading espe- 
cially if the system is oversimplified. 


The Experimental Model 


Few problems in the field of experimental dynamics 
present greater difficulty than that of devising a suit 
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Landing gear model. 
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able procedure for testing a theory of stability of a me 
chanical system. If the theory predicts stable oper 
ation for a given configuration and observation shows 
that an initial disturbance will die out, this is merely 
agreement the theory. 
Again, if the tests are conducted at a condition of neu- 


with—not confirmation of 
tral stability, this delicate balance may be dispropor- 
tionately affected by the ever present nonlinearities 
and other unknowns. For example, a | per cent error 
in the measured damping might well correspond to a 
20 per cent error in the neutrally stable speed. A 
summary of possible test procedures follows: 

(1) Measure all the system parameters and, applying 
the Hurwitz® criterion to the characteristic equation, 
determine whether the system is stable or not. Ob- 
servation will agree or disagree with the theory. At 
best, this method provides only a rough qualitative 
check on the theory. 

(2) Plot a graph showing the relation between the 
theoretical damping coefficient C, and the taxi speed |’ 
for neutral stability. Then, with the damper adjusted 
for a given value of C,, adjust the speed | until a sus- 
tained oscillation is developed. This method supplies 
a quantitative check on the theory. 

(3) A method presently under investigation at the 
Wright Air Development Center is analogous to the 
frequency response technique employed in servomech- 
anism With the 
justed for stable operation, a harmonic excitation sup- 
plied by a mass attached to the rotating wheel induces 
a steady harmonic response at each speed which may 


research. system parameters ad- 


then be compared with the corresponding theoretical 
value. 

Most of the tests on the WADC model were made 
using the second method. In order to reduce the effects 
of possible compensating errors in observing stability, 
precise measurements of the sustained frequencies were 
obtained and compared with the theoretical values. 
In addition, the phase relation between the wheel mo 
tion and the latteral motion of the pivot provides an 
exacting check on the theory. 

As shown in Fig. 8, the model consists of a wheel and 
rubber tire on an axle which in turn is supported by a 
This assembly is carried on a vertical spindle 
The lower 


trail arm. 
which passes through two ball bearings. 
bearing is fastened to the spindle and transmits a spring 
load to the wheel forcing it against the driving drum, 
thus simulating the airplane weight. The upper bear- 
ing is fixed to a sleeve (the airplane strut) which in turn 
is pivoted to the frame of the machine. The lateral 
rigidity of the nosegear is created by means of helical 
the the 
springs are carried on a collar which surrounds the 


springs between frame and sleeve. These 


sleeve and which can be adjusted relative to the sleeve 
pivot. In this way the lateral rigidity A, can be ad- 
justed to any predetermined value by simply selecting 
the appropriate position of the collar relative to the 


pivot. The spindle passes up through the sleeve and 
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Model and test equipment 
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is attached by means of a spring to a viscous damper 
This spring supplies the necessary degree of freedom in 
torsion. The base of the damper is mounted rigidly to 
the oscillating strut, thus permitting damping action 
only against the swivelling motion of the spindle. 
When in operation, the sleeve which supports the 
spindle oscillates on its pivot thus providing the de- 
gree of freedom corresponding to the lateral motion of 
the nosegear pivot. To simulate the airframe, the os 
cillating motion of the pivot drives a massive ring 
through an elastic shaft. Smaller masses, attached 
by flat springs to the main ring, simulate higher modes 
of motion in the airframe. Altogether, the model is 
provided with six selective degrees of freedom. Fig. 9 
shows the complete test setup. Fig. 10 shows the two 
degree-of-freedom airframe simulator. 
the accurately adjustable silicone-fluid viscous damper. 

The following adjustments can be made: Various 
wheels and tires; 0- to 4-in. trail; vertical load; lateral 
pivot rigidity; spindle-damper torsional rigidity; damp- 
ing from zero to several times critical; infinite damping 


Fig. 11 shows 


coefficient; strut mass; drum speed from 0 to 400 in. 
sec.; airframe simulator modes. In addition, any of the 
degrees of freedom can be suppressed at will. 

An autosyn attached to the top of the spindle is used 
to measure the swivelling motion of the wheel. A simi- 
lar instrument attached to the sleeve pivot measures 
the lateral motion of the nosegear. These two readings 
are recorded on an oscillograph thus permitting accurate 
measurements of stability, frequency, amplitudes, and 
phase relations of the two principal motions. The 
drum speed is recorded separately. 


Model Theory 


A precise mathematical description of the model 
requires the use of polar coordinates to specify the 
lateral motion of the strut since this element is carried 


on a fixed pivot in the model. In addition, in order to 
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permit accurate prediction of shimmy frequencies over 
the wide range of operating conditions for which the 
model can be adjusted, two parameters not included in 
the published report are given here—gyroscopic and 
vertical load effects. 


Notation 

3 = coefficient of yaw, rad. Ib 

C = tire time constant, sec. 

C = coefficient of viscous damping, in.-Ib. sec 

D = differential operator 

F = ground reaction normal to wheel plane, 
lbs. 

F, = ground reaction parallel to strut, Ibs 

I] = distance from trunion to axle, in 

l = polar moment of inertia of rotating wheel, 
Ib. in. sec. 

I, = moment of inertia of the system about the 
strut axis, Ib. in. sec.’ 

i, = moment of inertia of the swivelling parts 
about the spindle, Ib. in. sec.* 

K = torsional coefficient about the strut axis, 
in. Ib. /rad. 

Ko = K+ Wolo + Wil, + Well 

K, = torsional coefficient of the damper link 
age, in. Ib. ‘rad. 

L = trail length, in. 

Lo = distance from strut axis to c.g. of strut, in. 

Ly = distance from strut axis to c.g. of spindle, 
in. 

m = wheel mass, Ib. sec.” in. 

p = lateral compliance of tire, in. Ib. 

r = wheel radius, in. 

S = spring load, Ibs. 

1” = taxi velocity, in. sec. 

i" = S+ W, 





Airframe simulator 


Fic. 10. 
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Inspection of Fig. 12 reveals the presence of the fol 
lowing moments about the strut axis: 

(1) F,U7 + r), due to the ground reaction acting 
perpendicular to the wheel plane. 

(2) F,L6, F, is in the wheel plane but does not pass 
through the strut axis because of the trail angle 6. 

(3) 1, & the inertia torque of the whole system if the 
trail arm did not pivot about the spindle axis. 

(4) mLI/ 6, the c.g. of the wheel mass feels an acceler 





Fic. 11. Model damper. 


ation due to 6 equal to L6. The resulting inertia force 


Wo = strut weight, lbs. 
W, = spindle weight, Ibs. acting on a moment arm // about the strut axis produces 
W. = wheel weight (mg.) Ibs. a torque superimposed on that due to @. 
sie = angular velocity of wheel, rad. ‘sec. (5) w@, due to the wheel spin axis precessing in the 
oa = angular deflection of strut, rad. horizontal plane creating a gyroscopic couple in_ the 
A = lateral deflection of tire, in. vertical plane. 
6 = angular position of trail arm, rad. (6) Wolo + Wil, + We (Ha + Lé), the pendulous 
A = equivalent angular position of damper, restoring moment due to the weight of the parts. 

rad. (7) Ka, the restoring moment due to the elasticity 
(@ — 6) = angular strain in damper linkage, rad. of the strut and its support. 


Summation of these moments about the strut axis gives: 


F,(H +r) + F,LO — La — mLH6 — [06 — (Wolo + WiLia — WA(Ha + L6) — Ka = 0 i) 
Summation of moments about the spindle axis gives: 
F,L — 1,6 — mLH& — C6, + Iwa — WeoLa = 0 (10 
Equating damper torque to damper strain times the torsional elastic coefficient : 
CO, = K, (6 — 6) (11 
In Eq. (5) above, replace y from (1) by 
¥ = —CF,/(1 + GD) 
and with A = — pF, obtain (see also reference 1): 


(CipD? + pD + CV)F, + (7 + r)CD? + (1+ r)Dia+ [GLD? + (VC, + L)D + Vie =0 (12 


(It should be noted that this polar form of the system equations reduces to the coplaner form given in reference | 
when H approaches infinity. As previously remarked, the effect of this transformation is to eliminate the gyro- 
scopic stabilizing action of the wheel. For normal configurations, and at low speeds, the gyroscopic action is negli- 
gible.) 
The determinant of these four simultaneous equations yields the seventh-order characteristic equation of the 
model which may be written 
A;D'0 + AsD*9 + A;D*0 + AyD*O + A;D*0 + A.D°O + A,DO + Aid = 0 (13 
where, with lV. = 0 
Ao = KoK,LV 
K,[KoL? + IlwV(H + r) + WL?(H + r)] + C.Ko V(L + CK,) — CWLIwVK, + C.KoK,LV 
K LVI, — mH(H + r)) + C,[KoL? + JwoV(T + r) + WL +r) + KA + 7r)?] 
+ pIC.AoK, — WLIwK,| + CG[C,KoLV + Ay KL? + lwVK,(H + r) + WL*(H + r)K;] 
+ CV[|K,7,Ky + Pw? + mL?WH) — C,WLIe) 
A; = K,[/,(71 + r)? + 1,L? — 2mL?H(H + r)] + (C, + CA, [ly — mH + r)| LV 
+ C,CVU,K, + 1;Ko + [?w?) + C,CQi[AoL? + AH +r)? + lwV(H +r) + pRKoK,]) 
+ pK,I,;Ko + Iw?) + WL*[pK mH + C,CmHV + C,C\G7 + 7))] — (C; + C.K.) WLIwp 
A, = (C, + GK,) 7 + 1)? + T,L? — 2mL?H(H + r)| + CK,VU,1, — m?HP?L?*) 
+ Cipl,K, + (C, + GK) pUsKo + Pa? + mL?HW)] — C,ClwpWL + C,CQLV[l, — mH(A + 1r)] 
As = C,C\[[,47 + r)? + 1,1? — 2mL?H(H + r) + pU.Ko + 1K, + Pw? + mL2WH)] 
+ (C,CV + pK,) Ud, — mL?) 
Ae = (C, + CK, pd, — m?H°L?) 
A; = C,CipU ls, — m*H?L*) 
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Eq. (13) describes the nuances in the performance of the model with precision over the widest ranges of system 


parameters. 


reached for this case are actually reversed with more realistic configurations. 
the damping coefficient computed by Eq. (6A) will stabilize the system represented by Eq. (13). 


It contains the stability criterion found by assuming a ‘‘fixed’’ pivot; yet shows that the conclusions 


For all normal ranges of parameters 
For each valu 


of the lateral rigidity used in this equation there will be associated a sustained natural frequency provided the 


damping is adjusted to maintain neutral stability. 


may be superimposed on a similar graph of frequency vs. rigidity for the airframe. 
represent neutrally stable operating points for the complete system as described in reference 1, p. 47. 


A graph of this frequency vs. the lateral rigidity of the strut 


The points of intersection 
In order t 


determine whether the coefficients in Eq. (13) will provide stable operation, Routh’s criterion® (as given by Hur 


witz) should be used. 


For the purpose of electronic analog studies, it is simpler to use the four simultaneous equations rather than the 


single Eq. (13). 


But the latter is more suited to the study of particular cases. 


For example, the particular case 


in which ~p approaches zero and the motion of the damper piston is suppressed, (C, = ©) is described by 


Ay = As = A; = 0 
A, = CVA)K, + AyLV 


A, = Aol? + Af +7)? + WL +r) + JwV(A +7) — CVWIeLk + GLVA, 


A; = LVI, — mH(H + r)|) + CV[L,K, + [Ko + Pw? + WHmL?| + C[AoL? + A(T + r)? 


+ WL°(H + r) + lwV(H + r)] 


Ay=1(H+ 7)? + 1,0? — mL’?H(H + r) + GLV{[, — 
As = CVI), — m?H?L?|] + C77 + 7)? + TL? - 


where W = $+ WW 


NOTE: 
system of the main landing gear because in this case the 
torsional motion of the wheel is opposed by the tor- 
sional elasticity of the strut. There is no damper. It 
is important to note that the main gear, although not 
freely swivelling, may shimmy if the stability criterion 


These equations now describe the physical 


is not satisfied. 
The criterion for neutral stability for this case is! 


A oA 3A {= A 1A ‘7 + A»*As (14) 
and the associated sustained frequency becomes! 
w* = A>» Ay ( 15) 


With the model damper locked and the surface speed 
adjusted until a neutrally stable oscillation is estab- 
lished, these expressions can be used to obtain two sep- 
arate checks on the accuracy of the theory. Other spe- 
cial cases may be tested in a similar manner. 

Table 1 contains a summary of the various nose- 
gear configurations for which over 2,000 analog studies 
and 100 model runs have been completed. For ex- 
ample, Case 1 describes a structure for which the tor- 
sional rigidity is infinite and the lateral motion of the 
strut is completely suppressed. Only the degrees of 
freedom described by the spindle swivel (@) and the tire 
deflection (A) are present. For this simple case the 
coefficients of the general Eq. (13) reduce it to a third 


order equation, 
pl.6 + CVI6+ L9+ LV@ = 0 


and for stability the product of the coefficients of the 
middle terms must exceed the product of the end terms, 
i.e, LC > p. This is the stability criterion advocated 
by Riekert and also by de Carbon. Case 2 differs from 
Case 1 only in the presence of the tire time delay factor 


mH(H + r)] 
2mL?H(H + r)] 


C;. It is obvious that the stability criterion for this 
system may upset the predictions for the simpler system. 
In Case 6 the pivot is allowed to have lateral motion 
(aw) while the tire is assumed rigid (pb = 0). The sta- 
bility criterion is now markedly different from that 
given for Case 1. In particular, it is sensitive to 
changes in Ay, /,, and V, factors which do not appear 
The influence of these factors 
can be demonstrated by model tests. Case 13 is a close 
approximation to the main gear structure. The sta- 
bility criterion may be readily formed from the general 
equation by substitution of » = 0 and retaining only 
those terms which contain C,,. 


in the simpler criteria. 


TEST RESULTS 

Since the equations given in the report! are expressed 
in nondimensional form, the problem of selecting the 
model parameters to represent any given prototype 
gear is simple. For example, the proper operating 
speed for the model is found by equating its nondimen- 
sional speed to that of the prototype. As shown on p. 
17, reference 1, speed is nondimensionalized by multi- 
plying it by C./A,L. If the prototype and model di- 
mensions are, respectively, 2 = 5 in. and 2 in., A, = 
0.5 X 108 and 50, 7; = 3 and 0.005, then for a proto- 
type speed of 3,000 in./sec., the model speed is 295 in. 
sec. In this way over 100 model configurations simu- 
lating various airplane gears were tested. Typical re- 
sults are shown in Figs. 13 and 14. 

Fig. 13 compares the predicted values of shimmy fre- 
quency, Eq. (15), with the measured values on the model 
for several values of the torsional rigidity. It should 
be noted that the close agreement between the two 
extends over a range of parameters several times 
greater than that covered from the lightest type fighter 
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to the heaviest cargo planes (Runs 55, 111, 113, 114, 
Lh¢, FES). 

An even more searching test of the theory is its ca- 
pacity to predict whether a given combination of param- 
eters will produce a stable system. The ordinate of the 
plotted point of Fig. 14 is the ratio of the left to the 
right side of the stability equation, Eq. (14). The 
test was made by adjusting the drum speed of the model 
until a neutrally stable vibration was established. The 
measured model parameters were then used to com- 
pute the stability ratio. If both the theory and meas- 
urements were exact, this ratio would be unity. The 
absolute size of the error carries no meaning unless it is 
related to some significant parameter. For example, 
in Fig. 14, if a plotted point shows a stability ratio of 
1.01, this | per cent error corresponds approximately 
to an error in damping of | per cent of the critical 
i.e., the error could be corrected by a rela- 
Equally 


damping 
tively small change in the measured damping. 
accurate results have been obtained with the variation 
of other parameters such as trail, wheel mass, strut mass, 
and wheel load for both rigid wheels and rubber-tired 
wheels. 

When the results shown in Fig. 13 and Fig. 14 are 
considered together, they provide conclusive evidence 
that the theoretical equations describe the model with 
satisfactory precision. 


The Four-Degree-of-Freedom System 


When an analog or digital computer is available, the 


complete 15-parameter nosegear system should be 
examined. Table 2 contains all the parameters for 
nine airplane configurations tested on the model. For 
each test the coefficient of damping (C,) was held con 
stant and the drum speed adjusted until a precisely 
neutral oscillation was established. The coefficients of 


Eq. (13) were then computed on an IBM Card Pro 
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grammed Calculator by causing it to search auto 
matically for the value of C, which would satisfy 
Routh’s criterion for neutral stability. As a further 
check on the accuracy of Eq. (13), the roots of the 
characteristic equation were computed, and from them 
the neutrally stable frequencies were found. The 
agreement between the experimental values and the 
theoretical values is within the limits of accuracy of the 
measurements. 

By assigning various velocities and computing C;, a 
complete graph similar to Fig. 6 can be constructed for 
each configuration. In this way, the most desirable 
proportions for the gear to achieve any degree of sta 
bility may be readily determined. 

Following is a summary of the model test results: 

(1) A nosegear can be made to operate in a completely 
stable manner with no artificial damping. This result 
was achieved with both the rubber-tired wheel and an 
aluminum wheel without a tire in complete accord with 
the theory. The stabilizing action is due to the gyro- 
scopic and vertical load effects. 
measurable part in this remarkable action since both 


Friction plays no 


the spindle and pivot are mounted in nearly friction-free 
ball bearings. To further test this theory, the gyro 
scopic effect was increased by mounting a small high 
speed gyro on the spindle with the result that the sta 
bility became greater as the speed of the auxiliary gyro 
was increased. The effect was even more dramatically 
demonstrated when, with no other changes in the sys 
tem, the direction of rotation of the gyro was reversed 
resulting in an instantaneous and violent instability. 
Possible applications are being explored. 

(2) Stabilizing effects are obtained by each of the 
following parameter changes: Decreasing the wheel 
moment of inertia; decreasing the lateral rigidity of the 
pivot; increasing the mass of the strut; increasing the 
torsional elasticity ; 


rigidity of the damper-system 


TABLE 2 
Comparison of Observed and Theoretical Damping and Frequency for Nine Nosewheel Configurations 


Parameter —- ee : A 











10 11 12 13 

c 0.01 

P 0.006 

Ci 0.025 

W 20 

I 0.01 } constant for all runs 

H 15.3 

r 2.7 | 

Ky 2,700 

pa 0.5 } 

I; 1.45 2.06 2.68 2.68 
Be 0.0169 0.0245 0.0322 0.0322 
m 0.0035 0.0061 0.0087 0.0087 
Ko 12,200 12,200 12,200 9,680 
V (stable) 59 38 30 
Ce 0 84 0.84 0.84 0.84 
CA 0.88 0.88 0.83 
fs 6.6 5.0 4.0 3.5 
fs 6.8 5.0 3.9 3.3 
' Measured damping coefficient. 


Theoretical damping coefficient given by Eq. (13) 


3’ Measured shimmy frequency, cycles per sec. 


4 Theoretical shimmy frequency given by Eq. (13). 


—Run Number— ieee _aaamiineiads. 
14 15 16 17 18 
2.06 1.45 1.45 2.06 2.68 

0.0245 0.0169 0.0169 0.0245 0.0322 
0.0061 0.0035 0.0035 0.0061 0.0087 
9 , 680 9 ,680 5,830 5,830 5,830 
50 127 187 63 40 
0.84 0.84 0.84 0.84 0.84 
0.85 0.91 0.88 0.88 0.91 
4.5 4.2 6.5 3.8 3.0 
4.3 6.8 6.4 3.0 2.9 
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THE Story 


and by increasing the vertical load. Finally, increasing 
the drum speed is destabilizing except for artificially 
high gyroscopic parameters. All of these results are 
in accord with the theory. 

3) Qualitative check runs on de Carbon’s stability 


criterion 
URV? a) (aD jr) > 


were made as follows: 

A wheel and tire combination for which aD 
IR\* > a were placed on the model, and the system 
proved to be strongly stable in complete accord with de 
Then with no other changes in the sys 


S7 and 


Carbon’s rule. 
tem except a slight increase in the lateral rigidity of the 
pivot, the operation became unstable. This observa 
tion cannot be explained by the above stability criterion 
Next, 


the lateral rigidity was adjusted to its original stable 


since it does not contain the lateral rigidity. 


value, whereupon it was found that by tncreasing either 
/ or V the oscillations became divergent in direct contra- 
diction to the given stability criterion. These observa 
tions have been accepted by the research people who 
have witnessed them as proof of the inadequacy of a 
theory of shimmy based on the properties of the pneu- 


matic tire alone. 


Full-Scale Field Tests 

The aircraft selected for these tests was the Fairchild 
C-119 cargo plane. Cumulative operating experience 
with this aircraft has demonstrated that its nosegear, 
though structurally sound, is only marginally stable. 
The particular plane on which the tests were conducted 
was susceptible to shimmy at speeds above 60 knots, 
and sustained oscillations could be readily induced at 
lower speeds by striking a 2'/2-in. obstruction placed at 
a 45° angle to the runway. Examination of the hy- 
draulic damping and steering system and a systematic 
check for backlash failed to reveal any abnormal con- 
ditions. An analysis of the C-119 gear on the basis of 
the WADC theory of shimmy indicated that the cause 
of the poor performance of this gear was to be found in 
the unfavorable proportions of four structural features, 
the most severe of which was the excessive diametral 
moment of inertia of the wheel proper. The nondi- 
which, in general, should not 
exceed unity This 
high value together with the unfavorable trail ratio 
describes, as shown on p. 32, reference 1, a system 
which cannot be stabilized regardless of the amount of 


mensional inertia ratio 
is approximately 6 for the C-119. 


viscous damping employed. In practice, of course, 
the presence of Coulomb damping does permit stable, 
but generally undependable, operation. 

The experimental program was planned to provide 
a series of modifications to the gear and to conduct 
instrumented taxi runs to determine the influence on 
stability and frequency of each of the alterations. The 
first instrumented run was made with the normal con- 


figuration and resulted in an accident. The oscillo- 
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graph records, however, were intact. The divergent 
oscillations began at a speed of 60 knots. The runway 
obstacle was struck two seconds later, resulting in an 
increased rate of divergence. In less than two addi 
tional seconds, the amplitude built up sufficiently to 
shear off the damper anchor bolts. The entire air 
plane then shook violently until its forward motion was 
arrested. As indicated above, theoretical analysis 
shows that the most significant factor in determining 
stability is the diametral moment of inertia of the wheel. 
To test this point a 34-in. wheel was substituted for the 
original 44-in. wheel which represented a reduction in 
moment of inertia of 65 per cent. The first run was 
made at a speed of 60 knots and when no shimmy was 
observed the speed was increased to 120 knots with no 
tendency toward instability. The next series of runs 
was made striking the runway obstacle at speeds up to 
100 knots. Shimmy could not be induced, and the 
oscillograph records showed a rapid attenuation of the 
oscillation following the impact. 

The actual trail on the C-119 is 5.5 in. 
Under these circumstances 


Its nondi 
mensional value is 0.41. 
increasing the trail should be destabilizing [p. 32, 
Ref. (1)]. The results of tests conducted with a 7.91- 
in. trail showed shimmy starting at lower taxi speeds 
and rapidly reaching greater amplitudes than with the 
5.5-in. trail. (It may be noted that these results are 
again in conflict with the predictions based on the rule 
given by de Carbon which indicate greater stability for 
moderate increase in trail.) 

When the inertia ratio is high, it is predicted that 
the system will be unstable for very large values of C, 
as well as for values below a given amount. To exam- 
ine this point a series of tests were conducted using 
damper orifice sizes of 0.042 in., 0.0625 in., 0.086 in., 
0.125 in., 0.156 in., and 0.1SSin. For both the smallest 
orifice (high damping) and largest the shimmy was un- 
controllable; it occurred at lower taxi speeds and di- 
verged more rapidly, even without striking the obstacle, 
than was the case with the intermediate sizes. 

A nosegear structure which is unstable because of a 
high inertia ratio can be stabilized by increasing the tor 
sional rigidity (K,) of the damper system. This is not 
the most efficient method because on an existing struc- 
ture it is accompanied by an increase in weight and, as 
shown by Eq. (6A) above, requires somewhat greater 
damping. For the C-119 tests, a 20% increase in tor- 
sional rigidity was obtained by replacing the production 
scissors links with links of considerably greater weight 
and cross section. Results showed a slight but definite 
improvement in stability indicated by the higher speeds 
reached before the onset of self-sustained oscillations 
and greater attenuation of shock disturbances below 
these speeds. 

No attempt was made to decrease the lateral rigidity 
of the strut to improve the stability as indicated by the 
theory. But in one test, the lateral rigidity was in- 
creased by restraining the motion of the top of the strut 
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by angle irons riveted to a crossbeam. This structure 
was destroyed during the first taxi run. Whether this 
was due to the added vigor of the shimmy or to struc- 
tural weakness is not known. However, it was felt 
that it would be unsafe to pursue this particular type 
of test on an already unstable structure. (A 1948 Boe- 
ing investigation (D-9361) confirms in considerable de 
tail the conclusion that stability is increased by decreas 
ing the lateral rigidity.) 

Fig. 15A (Run 4) is a trace of the oscillograph record 
taken at the damper link during the first of the tests 
The nosegear configuration is nor- 
The 
system is obviously unstable and, in fact, the motion 


described above. 
mal, with its inertia ratio of approximately six. 


was divergent before striking the runway obstacle. 
Failure occurred shortly beyond this point. 

Fig. 15B (Runs 17, 18, 19) is the oscillograph record 
taken at the same point as in 15A, but here the inertia 
ratio has been reduced to two. In this case, the run- 
way obstacle was struck in an attempt to induce shimmy 
but the displacement was rapidly attenuated thus dem- 
onstrating the strong influence of the inertia ratio on 
stability. 

Fig. 15C (Run 21) is the record obtained for the same 
configuration as in Fig. 15B but with the damping co- 
efficient increased by reducing the damper orifice di- 
ameter. This is in accord with the theory wherein it is 
shown that for inertia ratios in excess of unity there is 
an upper limit of damping coefficient beyond which 
stability is not possible. As this condition is ap- 
proached, stability is reduced as is apparent from the 


record. 


CONCLUSIONS 


(1) A stable landing gear that will show satisfactory 
attenuation of disturbances, without the risk of over- 
damping, can be obtained provided that the inertia ratio 
[,,K,/NK, does not exceed unity. 

(2) When the above condition is satisfied, the damping 
ratio required for neutral stability is closely approxi- 
mated over a wide range of system parameters by 


QS | S 1]? 
411 +S) Q . Vlaa Ss) | 7 
R= : (GA) 
(3) A precise description of all the essential phe- 
nomena of shimmy requires a mathematical model of 
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5 degrees of freedom, namely, tire deflection, swivel 
angle, strut deflection, damper-linkage strain, and air 
frame motion. Eq. (13) is the most accurate represen 
tation of the dynamics of this system that can be ex. 
pressed in linear coefficients. It establishes a general 
ized concept of inertia ratio and supports the conclu 
sions given in (1) and (2). This latter result is due in 
part to the cancellation of the stabilizing gyroscopic 
action by the yawing action of the pneumatic tire. 

(4) Eq. (13) predicts the conditions under which co 
Although 
not usually susceptible to large-angle shimmy’ this 
construction suffers severe fatigue loading due to wheel 


rotating dual wheels are inherently stable. 


unbalance when the stability is marginal. 

(5) A completely stable free-swivelling landing gear 
can be designed with no artificial damping when pro 
vided with sufficient gyroscopic momentum. Whether 
this goal is important, or even desirable, depends on the 
application. For example, if power steering is to be 
used, very high damping is a simple adjunct and the 
gear will be stable provided its inertia ratio can be made 
sufficiently small. For this case and also for the main 
gear, Eq. (13) should be used to check the stability of 
high-inertia systems, using an analog computer or Hur 
witz criterion. 

(6) The stability of a landing gear is influenced in 
varying degrees by 15 parameters. Eq. (13) predicts, 
and experiment confirms, the quantitative influence 
of each of these parameters. 
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Development of Charts for Downwash 
Coefficients of Oscillating Wings of Finite 
Span and Arbitrary Plan Form 


MAX A. DENGLER* 
Northrop Aircraft, Ine. 


SUMMARY 


An outline is given for the development of downwash coefficient 
charts to be used for the theoretical prediction of circulatory 
spanwise loadings of oscillating wings of finite span. The charts 
ire based on lifting line considerations which could be thought 
of as the logical extensions and generalizations of the Prandtl and 
Weissinger lifting line theories to unsteady state. The use of 
the charts should greatly simplify and facilitate the amount of 
computational effort generally involved in a flutter analysis 
They seem to be fit for routine calculations in the hand of the 
iircraft engineer 

A brief outline of the determination of lift and moment co 
efficients of a finite span wing oscillating in bending and twist as 
weli as of the complete flutter analysis is added to illustrate the 


practical use of the charts 
DISCUSSION 


i igen FUNDAMENTAL ideas of the present investigation 
are based on the nonsteady state lifting line theory 
for wings of finite span developed in reference 1. It is 
thought that this theory, which is restricted to incom- 
pressible flow and subsonic speeds, may be considered 
as the logical extension and generalization of the Prandtl 
and Weissinger lifting line theories from steady state 
to oscillatory wing motions. The basic aim of the 
theory is the prediction of spanwise circulatory load- 
ings for wings of general form, specifically for tapered 
wings of high sweep and low aspect ratio as are of great 
interest in present-day aircraft design. 

Spanwise variations of circulation are needed for cal- 
culations concerned with stability and flutter. The 
present analysis is an extension and modification of the 
aerodynamic theory of reference 1. It is a modification 
in the sense that some of the underlying aerodynamic 
principles of the wing vorticity model employed there 
have been radically changed. It is believed that this 
change constitutes a substantial improvement of the 
theory, as seen both from the mathematical and physical 
point of view. The most directly apparent effect of 
this modification is the absence of any bound vorticity 
correction factor in the new version. It is the intuition 
of the writer that the new vorticity model, designated 
as ‘‘alternative (1),’’f is aerodynamically superior over 
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the original version as contained in reference 1, desig 
nated as “‘alternative (2)’’f from now on. 

This should not mean that the latter one should now 
be completely abandoned. On the contrary, as de 
rived from past investigations concerned with agree 
ment between theory and experiment (Ref. 2) the con- 
clusion must be drawn that alternative (2) yields 
values of satisfactory accuracy for flutter frequency 
No numerical results for alternative 
The pres- 


and flutter speed. 
(1), on the other hand, are as yet available. 
ent paper is concerned with both these alternatives. 
The present investigation is an extension of the earlier 
work on the subject in two different respects. First, 
it is one of the inherent characteristics of A(1) that it is 
in no way restricted to the use of one single lifting line 
for the total wing vorticity, but permits equally well 
the use of any finite number of lifting lines spread out 
A(2) is restricted to the use of 
One 


over the wing surface. 
one single lifting line by fundamental reasons. 
of the outstanding consequences of these circumstances 
is the fact that A(1) may be applied to wings of very 
low aspect ratio, A(2) not. 

The other aspect which distinguishes the present 
work from the contents of reference | is situated in a 
change of general perspective. This study is no more 
limited to purely theoretical considerations of aero 
dynamic character but intends to pursue the theory 
all the way down to the numerical ends. The aim was 
to overbridge the gap between the theory itself and its 
application and to open the problem of the oscillating 
wing of finite span to the aircraft engineer. Seen from 
the practical side the heart of a flutter analysis is the 
determination of the spanwise variations of circulation 
due to the various modes present in the actual wing 
motion. The circulation is directly related to the dis- 
tribution of downwash and the downwash velocities re- 
sult from the vorticity of the flow. It becomes appar 
ent at this point that it is the determination of the vor- 
ticity which should be systematized and simplified to 
the greatest possible extent. Downwash coefficients 
hence play a predominant role in any flutter analysis. 

This should justify the main objective of the present 
study which is concerned with the preparation of charts 
and tables for downwash coefficients of an oscillating 
lifting surface corresponding to any assumed type of 


vibration. The fundamental idea underlying such a 
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procedure is not new. The presently developed method 
has certain basic features in common with the principles 
of F. W. Diederich (Ref. 3) with regard to downwash 
coefficient charts for finite span wings in steady state. 

The fundamental concept of the Diederich method is 
the use of a finite number of discrete rectangular horse- 
shoe vortices arranged along the forward quarter-chord 
line, having distinct intensities and vortex trails ex 
tending downward to infinity, parallel to the direction 
of flow. The downwash induced by the vortex system 
is to be equal to the geometrical downwash of the wing 
motion assumed. This condition, called downwash 
condition from here on, should be satisfied at all points 
of a lifting surface. An approximative solution of 
reasonable accuracy is obtained if the criterion is satis 
fied at a selected sequence of wing points only, the 
latter to be called downwash points. It follows from 
theoretical arguments of Weissinger, and 
others that favorable approximations should be ob- 
tained if only a few discrete points at the rear quarter- 
chord line are taken as downwash points. 

It was analyzed in reference | that these same three- 
quarter chord points should yield good approximations 
to the solution of the underlying integro-differential 
equation for the circulation also in the corresponding 
case of oscillatory motion if at the same time rectangu- 
lar horseshoe vortices of oscillatory intensity are aligned 
along the forward quarter-chord of the wing. A(2) is 
based on this scheme. There is, on the other hand, no 
reason why the total wing vorticity should be repre- 
sented by a single lifting line only. A(1) accordingly 
permits the use of a finite number of elementary vor- 
tices arranged along lifting lines (vortex lines) extend- 
Downwash 


Pistolesi 


ing in spanwise direction over the surface. 
points have now to be specified at more than one line, 
the number of downwash lines being equal to the 
number of vortex lines and the total number of down- 
wash points being equal to the total number of vor- 
tices. There are reasons that also in this case the 

4)-line should be retained as one of the vortex lines 
(the principal one) and that the (3c/4)-line should be 
kept as one of the downwash lines (the principal one). 
The question where to place the others is difficult to 
answer. The specific choice should not matter too 
much, but further theoretical clarification of this point 
The basic intuition of the pres- 
validity of the results must in- 
and the number 


would seem desirable. 
ent theory is that the 
crease with the number of vortex lines 
of vortices composing them as long as the (c/4) 
maintained as the principal vortex line and the (3c/4)- 
line is retained as the principal downwash line. 

The total induced downwash due to the whole sys- 
tem of vortices is to be equal to the geometrical down- 
wash corresponding to the assumed wing oscillations, at 
If the idea of narrow- 


-line is 


any of the downwash points. 
span horseshoe vortices is applied to harmonic oscilla- 
a new parameter related to the reduced fre- 

of the oscillation enters the picture. While 


tions, 
quency 
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one single chart is sufficient to represent downwash 
coefficients in steady state (Ref. 3), a number of charts 
for different values of the frequency parameter is needed 
to serve the same in unsteady state. How 
this may be done in the perspective of efficient design 
procedures in the hand of the flutter engineer will be 
shown. It will become apparent then that the advan- 


tages of the rectangular short-span horseshoe vortices 


purpose 


which are characteristic for steady state calculations 
are even more pronounced in unsteady state consider- 
ations. These advantages are linked to the fact that 
such vortices may be employed in case of low aspect 
ratio wings (e.g., delta wings) just as well as in the case 
high aspect ratio wings (e.g., the “Scorpion” 
feature that all of the 
permits 


and 


other low drag airfoils). The 
elementary vortices do have the same span 
a scheme of nondimensionalization which makes them 
applicable to wings of any actual dimensions. 

The bound vortex sections of the vortices are directed 
normally to the direction of flow and hence approximate 
the lifting lines in a staircase-like manner. Errors from 
this source introduced into the numerics should be of 
minor importance only in the light of other more severe 
simplifications of the theory as a whole as are the as- 
sumptions of prevailing potential flow and the neglect 
of all viscous effects in general. The use of one or two 
vortex lines, each carrying eleven vortices, should be 
sufficient in Increase of the number of 
vortex lines and the number of elementary vortices will 
increase the accuracy of the solution but, of course, 
the computational labor at the same time. Symmet- 


rical wing oscillations may be analyzed as well as non- 


most cases. 


symmetrical ones. 

The vortex schemes according to both alternatives 
Local semichord is }, root 
Wing 
plan 


are shown in Figs. 2 and 6. 
semichord is d), and wing tip semichord is b, 

semispan is s and aspect ratio is AR = (span 2s)’ 
form area. Wing and wake points are referred to a 
Cartesian coordinate system (x, y) with center at point 
of intersection of forward quarter-chord and wing axis 
of symmetry, the x-axis extending downstream and the 
y-axis being directed normally to direction of flow. A 
dimensionless spanwise coordinate is introduced ac- 
cording to y = n-s. The angle of sweep ¢ is supposed 
to be defined by the inclination of the local tangent on 
the (c/4)-line. Sweep is constant in case of a tapered 
As seen from Figs. 2 and 6, all elementary vor 
Middle points of the 
called vortex centers. 


wing. 
tices have the same span, 27. 
bound-vortex sections will be 
They are located at the same spanwise wing stations 
as the downwash points. Vortices are numbered by 
M1, i, N2, i, ..., downwash points by mm, i, Me2,%, . « 
Coordinates of the forward quarter-chord line mez waved 
at the specified wing stations are denoted by p; and qi. 
The theory of the oscillating airfoil requires that the 
wake of any of the elementary vortices building up the 
wing vorticity distribution starts at the trailing edge 
deliberately intro- 


point. The theory of reference 1 
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duces the assumption that the wakes of the vortices 
start at the respective quarter-chord points. This 
simplification considerably reduces the mathematical 
complexity. It was reasoned, on the other hand, that 
the inaccuracy or “‘‘error’’ introduced from this source 
could be compensated satisfactorily by introduction of 
a bound-vorticity correction factor n(k) which would 
be independent of the specific wing motion and de- 
pendent only on the reduced frequency & of the oscilla 
tion. Numerical results deduced on this basis (Ref. 2 


seem to justify this hypothesis. However, any “‘cor 
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whatsoever should be eliminated from 
This idea was the guiding thought 


rection term’ 
a theory if possible. 
for the formulation of an alternative vorticity model 
A(1)) which will be developed now. 

According to the theory of the oscillating airfoil in 
two-dimensional motion developed by Cicala, Glauert, 
Kiissner, Reissner, and Theodorsen the velocity distri 
bution of an airfoil in incompressible flow may be con 
ceived as being induced by a system of vorticity spread 
out over the lifting surface and its wake. The vor 
ticity moving with the surface is called bound vor 
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The fundamental properties of the vorticity distribution of a lifting surface. 
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ticity; the vorticity which is shed periodically at the 
trailing edge and is not participating in the wing 
motion is called wake vorticity. If components of the 


vorticity normal to and parallel to the direction of for 


ward motion are denoted by y = 7-e“ and 6 = 
z a . . “ ° . : 
6-e, € and n designating Cartesian coordinates as in 
Fig. 1, they obey the relation 

OVE, )/On = O4(E, n) OE (I 


at any point of the lifting surface and any point of the 
wake. The subscript B standing for bound vorticity 
and the subscript IV’ for wake vorticity, the distribution 
of 5,(g, 7) follows immediately once the distribution of 


Yx2(t, n) is known: 


ik: 
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Wing vorticity diagram corresponding to two lifting lines, based on alternative (1 


r ‘ - - ' - ‘ 
On{&, n) = 1OVa(§, n) On; dé (2) 
by) /2 


The wing circulation [ on the other hand is 


bo /T “ 


r= l-e™ = Valé, n) d&-e' = bQ-e = (->) 
} > 

where 2, as usual, is the circulation function. The 

analysis of the mechanism of vorticity shed at the trail 


ing edge leads to the further relations 


- ° nr w/V) [i 3(b0/2 
iC = 1 = he a 
I 
= Oly w/t 3(bo/ 2 
ow = é ») 
On 


and it 
satisf 
tween 
flowir 
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where 


The 
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and it is seen that also the wake vorticity components 
satisfy relation (1). Due to the absolute analogy be- 
tween the magnetic field induced by an electric current 
flowing in a closed wire on the one hand and the velocity 
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distribution induced by a field of a three-dimensional 
line vortex on the other, the induced downwash is de- 
termined by the law of Biot and Savart and for a rec 
tangular wing is therefore given by 


> ®3b0/2 _- - . - - 
i) twl __ J Ya\E, 7) | lees sv Ops, 7) (y oo ‘Tay | tw 
Ir: V-w'’ (x, y)-e* = , oe > di¢ dn-e™ + 
. LJ —i/2 ie = 2)° 4 et — ere f 
F f i iw/*) [¢—(3/2)8 (—iw/V)P(n) (x — &) + jdP(n)/dnj (y — 9 | sa 
. é ° — 7 dé . dn-e (6 
. tJ 30/2 ise = £7 + 1s = er f 


Assume now the wing vorticity consists of a single horseshoe vortex as indicated in Fig. 1 and that, accordingly, 
the bound vorticity components are 


i—-r<9< +t? 


(bo [(3/2) — r] < E< &[(B/2) — r] + (7 


fol 9 ba: wd 


] 0 


Ya(é, 7) = 


at all other points of the wing 


fool € —r<n<—rte, &l(3/2) —- tT] < & < bo(3/2)) 
bn(f, 7) = 0 .... at all other points of the wing (8 
(—V7 e....7r<n<rt+e &[(3/2) — r] < & < &(3 2) | 
The contribution to the total induced downwash velocity 47- V-e-w''’(x, y) arising in the bound vorticity then 


becomes, for e— 0, 





, ‘wi (i) 7 Pho [ (3/2 rl + € a 
tr: V-e- we (x, y) sect .- i. ” (s+ — §) { 
: =e (1/€) ; = yaya dE dn + 
(ho Vy) Jue f bo[(3/2) — rz] i(x — &)2 + (vy — 9)*}°" { 
e—y + ¢ j bo(3/2) (y + r) ) 
(1 /€) = ya de dnt 
. { J bo[(3/2 Niet — Sr SF + er f 
%, +. ®(bo(3/2 
j . (—1)-(y — r) ( 
Vl €) = _ 313/72 dé ¢ dn (9 
J» J bo[(3/2 Ais — PP + — vr} 
or 
as () bo V { jx — b[(3/2) — rj F 
w°-V°Wp (X, V)/OVY = ; s >. 3/290 T 
J+ <jx — bf[(3/2) — r]}2? + (y — 0)? >? 
3by/2 (y + 7) 3bo/2 (y — r) 
° ° 3/2 dt — -\9 : 3/2 dé (10 
TRG?) — 7) <2 = + esr reo J bof (3/2 — 1 <(x — ££)? + (vy — vr)? >" 
. l (y + ¥) Ly = 7) } 
fo? Vw, (x, y)/b Vy = J : — : . 
e— SAV (x — 8)* + (y + 1)? ve — 3s)? + YY — 7)" 
l j (x — 39/2) w@-s ' 
rs IVWix — 3bo/2)? + (y +r)? Vis —s)? + Oo+r 
] j (x — 3b)/2) (x — 2 ) 1 
< — { 
iol A ly (x — 3b)/2)? ++ (vy — 7 Vix —sz)?+(y- ry?! 
where 
z = 6 [(3/2) — 7] (12 
The downwash contribution due to the wake vorticity is of the form 
P (vy + &£ 
. ‘ i | i‘ 
for- V-wy” (x, vy)/bVy = ica , dar - 
F , 0 ; r — 3h,/2 — Xr aT ty r)*}° 
: (vy —-»? 
ices , - 73 , dr - 
J 0 (x — 3069/2 — A> + (vw — 1A’ 
*) . 9 ») 
, j " (x — 3b)/2 — AX | ? ane 
(tw / WV) é we 73 » Adi dyn (13 
. 0 (x — 36)/2 — Az + (y — 9)*t" 
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: : j (x — 39/2) (x 5309/2 
to? V-wy (x, y)/boVy = ) _ 
(y + 7) V (x — 3bo/2)? + (vy + 7)? (iy =F) Or 36/2)? ++ (y—r 
3/2 Iw (y +) (y —? 
| tiyt> y—r\) 
, Ww ; i] Ww | (y + F Ww 
in/2 + 4-3 2 = ie 3bo/2) ¢ + ci (x — 3bo/2) + I eh 
VI VI y+r V1 | 
if ua w ) 1 } . 48 = SOo/zZ 
: - f yr? + S a ae ae 
y-r (| f 3 (iy } | f 
(4 r C jx — 3h)/2 w 
“ es Ae oe 14 
y-r tiy r aw f 
where S2(t) and Cz(t) are the familiar functions ‘‘Sinus Integralis’’ and ‘‘Cosinus Integralis,’’ F(t) is the Cicala 
Function 
FO) = Gp | é kt — (Vs?+ f —t?) s{ ds 5 
70 
and S is defined by 
%a 
sa, fl) = | e (}*v (s? + 1) — 1} s) ds 16 
/ 0 


Assume now that new variables u, v, 4, and o are introduced according to 


(x — 36/2) = r-u + pP, y=rvt+q (17)* 
aw = (r/b)k = wr/V, o = (b/r)r (18)* 


The total induced downwash, obtained by combination of wing and wake contributions, may then be written in the 


form 
4en f 4 | iv te + «}° + @ + 1 V (u + o)? + (v — 1)7) 
. w' (x,y) = ; — 
y bh : (u +o) t (v + 1) (vy — 1) f 
: diate (= 2 4 f(v + 1) eo ne ' Ss Cil , 
lu-e : ‘ _ r ptw/l + 4-dt(uu) + Coluu); + 
nw (v? — 1) Uo +1 v— 1 
(v _ @— i] . (vy + 1) u 
J -F(u}v + 1)) - -F(ulv — | | + 5 J ,iot+ fiat — 
V2 v— v+i Uo + 1 f 
(vy — 1) j u 
: ; S4 ; 1 (v — Ilp = w*(o, wu, U, Vv) 19 


After a few more algebraic transformations the induced downwash appears as 


donor 1 jV(uto?+ +1)? V(ut+o)?+ ( — 1) 
-— << (x, ¥) = — > + D*(u, u, v (20 
y bo (u+a)\ (v + 1) (vy — 1) f 
where 
r ib d : 
D*(u, u,v) = | eo (1 7 2 log : 2 ; ) — iD(g,, b)| 7 
@tI1)L 2 i+ Vd + gi’) 


| th» Fo 2 
je —_ ( — — log = t ) — 1D (ge, bs) (21 
(v — 1) 2 . 2 Bt go?) § 


D(g, b) = e7"® $< br/2 + b-Im| F(b)} + B(gb, b) >i + <bV (1 + g?) + b[Re} F(b)! + 


log (2b) + w — 1] — A(gd, b) >} (29 
vob 
A(gb, b) = | (1 — cos u) } V (u? + b*)/u} du (23 
J0 le 
B(gb, b) = | (sin u) | V (uw? + b*)/u} du (24 
/ 0 


* The b appears here the first time as compared to the bo, p and g are the coordinates of the vortex center, and & is the local 


reduced frequency. 
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Re | F(b)} and Jm }| F(b){ .. . real and imaginary | 


parts of the Cicala Function, Eq. (15) (25) 
Y = 0.5772157 . . . Euler’s Constant (26) 
bh = ulvt+ i bo = plu — 1 (27) 
£i= u/jvut+ ll; go = u/|v — | (28) 


The Cicala Function has been tabulated by various 
Tables of the functions A(gb, b) and 
The reasons for these 


authors (Ref. 5). 
B(gb, b) are found in reference 4. 
purely mathematical transformations are evident: 
the ““D-Function D(g, b)’’ comprises all the principal 
transcendental terms of the downwash expressions. It 
is a function of two parameters only while covering the 
range of the three independent parameters yu, u, and v. 
And finally, it does not possess any singularity. 

Based on the foregoing mathematical developments 
downwash coefficient charts may now be developed for 
alternative (1) as follows. 

Consider the vorticity distribution of a_ tapered 
swept wing as shown in Fig. 2. There are two vortex 
lines, located at (c/4) and (c/2), each carrying eleven 
elementary vortices. There are two downwash lines 
located at (3c/4) and (3c/8). 
(Fig. 3) prepared for a specific value of the frequency 
parameter » = kr/b, containing numerous entries of 
D*(u, u, v) for various positive and negative values of 
uandyv = ., —10, —8S, —6, —4, —2, 0, +2, +4, 
+6, +8, +10, , the numbers D*(u, u, v) being 
located exactly at the geometrical points (uw, v). As- 
sume the transparent chart is placed above the wing 
vorticity diagram (Fig. 2) in such a manner that the 
0) coincides with the 


Assume there is a chart 


origin of the chart (wu = v = 
trailing edge point of the wing station of the vortex 
considered (Fig. 4). If 7 was chosen as the unit of 
length in both diagrams and if the » of the chart is the 
uw of the wing station considered, the major elements 
(i.e., D*) of the induced downwash due to vortices of 
this wing station may be read off the chart, directly 
and for all wing downwash points. Adding the few 
more terms of Eq. (20) the complete downwash ex- 
pressions are obtained. If the uw is not exactly the u 
of the station, interpolation is to be employed. 

It is evident from this that a great number for dif- 
ferent values of wu are to be prepared in order to render 
the process efficient. There are needed so many charts 
that linear interpolation will yield results of satisfac- 
tory accuracy. After the procedure is completed for 
one of the wing stations (and all vortices of this sta- 
tion) it is to be repeated for all other stations. It will 
be noted that the coefficients D*(u, u, v) are independ- 
ent of o and hence identical for all vortices of one and 
the same wing station. This is a necessary consequence 
of the fact that they represent wake contributions. 
Adding up terms belonging together, the total induced 
downwash at each of the downwash points is obtained. 
Equating induced downwash and geometrical down- 
wash, the latter being a mere consequence of the geom- 
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etry of the wing motion, the as yet unknown vortex 
intensities are found and from there the spanwise cir- 
culation distributions. Magnitudes and positions of 
the local circulatory lift forces may be deduced from 
there according to known relations of the two-dimen 
sional airfoil theory. The total aerodynamic forces 
and moments are obtained by adding the noncircula 
tory terms (“‘apparent mass and inertia terms’’) to the 
circulatory terms. Combining the aerodynamic char 
acteristics of the wing with the geometrical and me 
chanical parameters the flutter analysis of the wing 
may now be developed along conventional lines, finally 
leading to the flutter frequency and flutter speed for 
If the Mach Number of the flow 


is of moderate magnitude, results of satisfactory accu 


incompressible flow. 


racy may be deduced from there by application of the 
Prandtl-Glauert Rule or other adequate compressibility 
corrections. The computational effort thus seems to 
be reduced to a minimum, the advantages of the down 
wash coefficient charts are evident and it is felt that 
they will be welcomed by the general aircraft engineer. 
The development of the downwash coefficient charts 
for alternative (1) may now be considered as completed. 
The analysis will be continued by a brief discussion of 
the principal aerodynamic characteristics of the vor 
ticity model of alternative (2) which was developed in 
reference | and successfully applied in a set of extensive 
flutter 
theory and experiment as are reported in reference 2. 
Of course only these few aspects will be analyzed where 


calculations aimed at comparisons between 


the vorticity model of A(2) differs basically from the 
model of A(1). 

As was studied earlier it is an inherent property of 
the flow around an oscillating lifting surface that vor- 
ticity will be shed periodically at the trailing edge. In 
consequence, if a system of short-span horseshoe vor- 
tices is used to represent the wing vorticity, the wake 
of any of these vortices is to start at the trailing edge 
point of its station. This causes mathematical com- 
plications and it would be of computational advantage 
to shift the point where the wake is to start to the vor- 
tex center, i.e., to the forward quarter-chord point if a 
single lifting line located there is employed. The 
vorticity model thus obtained would not be adequate 
without further modifications compensating for the 
arbitrary change. The theory of reference | arrives 
at the conclusion that such a compensating modifica- 
tion may be accomplished along the lines of the follow- 
ing two principles: 

(a) As far as wake vorticity contributions are con- 
cerned the arbitrary shift of the origin of the wake may 
be compensated by a multiplication of all wake terms 


o (3/2)ik . . . 6“ . . 
by a factor e “ (a posteriori, i.e., “outside the inte- 
grals’’). 
(6) The bound vorticity terms should be multiplied 


(3/2)ik 


(a posteriori) by a correction factor n(k)-e corre- 
sponding to the requirement that the downwash due to 


a three-dimensional vortex of strength y at the three- 
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quarter chord point of its station should be equal to the 
downwash due to a two-dimensional vortex of same 
strength induced at the same point. 

One concludes then on the basis of known principles 
of the two-dimensional theory that 


- ; 
n(k) = e **/?. ~ mk: Ly) (k) — 1 (k)f + 


2 
: —ik jim - oF “2 ( 9 
ike 1 = + 2-Si(k) + Ci(k) ¢ (29) 
i j n(k) } 
(42/y)w'(x, v) (r/bo) = i 
(v - 1) le, Vv (1 +r 21") 
l £1 bh ( 1b, 
< +e “(] — log 
(v+1) SV + g,?) 2 
l Jo 1b. 
= + es ( _ “4 log 
(»— 1) SV (1 + g?) 2 


Downwash coefficient charts for the vorticity model 
of alternative (2) could now be developed along lines 
rather similar to the basic principles of the charts for 
alternative (1). The charts would now have to be 
prepared for the coefficients Q(u, wu, v) and, again, for 
various values of the frequency parameter 4. Assume 
now that a chart as shown in Fig. 7 contains numerous 
entries for Q(u, u, v) placed along the ‘‘v-lines,”’ i.e., for 
various values of « and even values of v. The vorticity 
diagram of the wing, on the other hand, should be as 
shown in Fig. 6. The transparent chart corresponding 
to a certain value of uw is now placed on top of the wing 
vorticity diagram in such a manner that the origin of 
the chart (« = v = 0) coincides with the forward 
quarter-chord point of the wing station of the con- 
sidered vortex, i.e., with the point (x = ~, y = q) (Fig. 
S). 

The downwash velocities due to this specific vortex 
induced at all the different downwash points (which are 
located at the rear quarter-chord line) may now be 
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Fic. 5. Bound vorticity correction factor n(k) of alternative (2) 


Assume now that the two relations (17) defining di- 


mensionless quantities u and v are changed into 


—- p= 4&1; y-—q=u-r (350 
while Eqs. (15) and (16) and Eqs. (22) to (28) remain 
unchanged. It can then be shown that the induced 
downwash of the vorticity model of alternative (2 


becomes 


O(n, u,v (31 


established by reading from the chart the coefficients 
Olu, u, v), adding n(k)-P(u, uv, v) and multiplying the 
two terms by e**’? in line with relation (31). In most 
cases linear interpolation with respect to u will be needed 
and many charts are to be prepared to render the 
method accurate and efficient. 

The development of the downwash coefficient charts 
has now been completed. The next two paragraphs 
are to contain a brief outline of the derivation of lift 
and moment coefficients for a wing of finite span os- 
cillating in bending and twist, as well as a complete 
flutter analysis based on these developments. 

Consider a swept wing of a geometry as shown in 
There are two Cartesian coordinate systems, 
The origin of the first is located 


Fig. 9. 
(x, y) and (xv, ¥,). 
at forward quarter-chord, the origin of the latter at 
the point of intersection between elastic axis ‘‘e.a.”’ 
(which does not have to be straight) and wing plane of 
symmetry. The actual wing oscillation may be 
thought of as superposition of bending oscillations 
about the hypothetical axis (A,B,) and pure torsional 
oscillations about the elastic axis and may approxima- 
tively be represented by the first normal modes of a 
cantilever beam in bending and in twist. Thus we 
have: , 


(a) Bending about hypothetical axis (A,B,) 


2(n,) = p(t)-bo-fr(n, 5, saga ae \ (29) 
2(n,) = 0 pe bis 

(b) Twist about elastic axis (e.a.) 
O(n.) = Q(t) -filne) ... Ig > 0) - 
O(n.) = 0 ne <0 a 


The nondimensional distances 7 and n, are related to 
each other by 


n, = n — (d-b/s) sin ¢g-cos ¢ (34) 
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DEVELOPMENT OF 


while the coordinates x, y and x,, y, are satisfying the 
equations 


x = y,‘sin g + x,-cos g + dy-bo \ 
y = y,cos g — xX,-sin ¢ 


(35) 
J 

dy stands for the value of d(n) at wing root. In case of 
a tapered wing of straight elastic axis d(n) is constant 
over the span. The first normal modes of a cantilever 


beam are 


(Bending) .. . f,(n,) = 0.5 (cosh on, — cos on,) — 
0.367 (sinh on, — sin on,) (36) 
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.) = sin | (n/2)n,} (37) 


¢ 


(Torsion) .. . fo(n 


where o = 1.8751. 

Translatory and angular displacements of a wing ele- 
ment PQ at the forward quarter-chord point R will be 
denoted by /(n) and a(n). his positive when directed 


downward, a is positive when nose-up. Introducing 
the amplitude functions 
p(t) = poe; q(t) = q-e™ (38) 


and the mode functions 


go(n) = foln.e = fotn — (db/s) sin g- cos gj (39) 
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Fic. 6. Wing vorticity diagram based on alternative (2) 
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Fic. 7. Downwash coefficient chart of alternative (2). 

















DEVELOPMENT OF CHARTS 


g(n) = filnde = fain — (dbs) sin ¢- cos gt | 10) 


the displacements /(n) and a(n) are found to be 


h(n) = p(t)bo-g,(n) — g(t)db- cos ¢-g4(n (41) 
a(n) = p(t) (b/s) sin ¢g-cos ¢-g,"(n) + 
q(t) cos ¢- g,_(n) (42 


where g,'(n) = dg,(n)/dn = dj) fp(n.){ ‘dn. According 
to Eqs. (41) and (42) the total wing oscillations may be 
conceived as superposition of four partial oscillations 
corresponding to translatory displacements /,(y) and 
Des- 


ignating the geometrical downwash corresponding to 


h(n) and rotatory displacements a,(y) and a,(n). 


a wing oscillation 2(x, y, t) by 


y) = Os/dt + I'(0s/0n) (43 


FOR 


DOWNWASH COEFFICIENTS 821 
and employing adequate subscripts, there are the fol- 
lowing basic relations. 

(a) Translatory oscillations corresponding to a linear 
caused by 


vertical displacement h,(n) = pobog,(nje™, 


wing bending about axis (A JB, 


= potko: g(n), 
(independent of chordwise coordinate x) 


Geometrical downwash wy, °'(.v, 


14 


(b) Translatory oscillations corresponding to a linear 


vertical displacement h,(n) = —gdb cos ¢-g,(ne™, 


caused by wing twist about elastic axis 


Geometrical downwash wy, °'(x, y) = 
—qyd(b by) cos ¢-tRo: go(n 


(independent of chordwise coordinate x) 15 
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Superposition of wing vorticity diagram (Fig. 6) and downwash coefficient chart (Fig. 7) for alternative (2). 
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Fic. 9. 


(c) Rotational oscillations corresponding to an angu- 
twt 


lar displacement a,(7) = po(bo/s) sin g-cos ¢-gy'(ne™, 
caused by wing bending about hypothetical axis (A ,B,): 


Geometrical downwash 
Wap” (x, Y) = po(bo/s) sin ¢- 
cos ¢: { iko(x — y-tan g)/bo + 1} -gp'(n) (46) 
Geometrical downwash at rear quarter-chord point 
Wap” (y tan g + b, y) = 
Po(bo/s) sin yg: cos g- (tk + 1)-g,'(n) (46a 
(d) Rotational oscillations corresponding to an angu- 
lar displacement a,(n) = qo cos og (ne™, caused 
by wing twist about elastic axis: 
Geometrical downwash 
Wag” (x, y) = 
go* COS go: fiko(x — ytan g)/b) + 1 -g(n) (47) 





s 

















Geometry of a swept wing oscillating in bending and twist. 


.” 
fo 
Geometrical downwash at rear quarter-chord point 
Wao (y tan g + b, y) = 
go:cos g: (tk + 1)-g,(m) (47a) 


If it is assumed for simplicity that a single lifting 
line located at forward quarter-chord represents the 
total wing vorticity, both the alternatives (1) and (2) 
apply directly to the present case. The downwash 
conditions then are to be satisfied at selected points of 
the rear quarter-chord line only. The distributions 
of geometrical downwash are given in Eqs. (44) to 
(47a). The distributions of induced downwash are 
found from Egs. (20) or (31). Equating geometrical 
and induced downwash at all of the downwash points, 
the vortex intensities Yr», Yq Yap: ANd Yoq are found. 
Using the general symbol y;; for these intensities, the 
wing circulation I';;(7) in each case is found from 
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DEVELOPMENT OF CHARTS FOR DOWNWASH 
Dij(n) = by: Vee + y5;(n) (48) a(n) -Lagln) = go (1/2 — t/k) cos g-g,(n) — 
: ; (i/mk)e™ Ky(tk) (bo/b) (Y¥aq/G) te" (55 
As known from the two-dimensional theory of the os- 
cillating airfoil, on the other hand, circulatory lift hy(n)>Mny,(n) = pobo(1/2)-gp(n)-e (96 
L ire() and circulation ['(y) are related to each other he(n)* Magn) = —qod(n)b(n) (1/2) cos o-ge(n)-e™ (57 
according to 
a(n) Map(n) = polbo/s) (3/8 — t/k) sin ¢- 
or ee a a ey | ry QO) ’ beat i 
Leire(n tk-e-Ky(tk)pV-T(n (49 cos 9+ g_'(n) -e™ 58 
where A, stands for the modified Bessel Function of the a(n): Man) = qo(3/8 — i/k)g.(n)-e- cos ¢ (59 


second kind. The total lift and moment distributions 
are obtained by adding the two-dimensional noncircu- 
latory contributions to the circulatory terms. The 
lift acts at the forward quarter-chord and is counted 
The pitching moment is counted 


Total lift and moment per 


positive upward. 
positive when nose-up. 
unit span (measured normal to direction of flow), re- 
ferred to the forward quarter-chord point, may be 


written 
L(n) = mpb?w? < hy(n)Liny(n) + hen) Lngn) + 

bj} a(n) Lapln) + ag(n)Lag(n)i> (50) 
M(n) = mpb a? < hy(n)May(n) + y(n) Mig(n) + 


b+} a(n) Map(n) + a(n) Mao(n {> (51) 
where lift-and-moment coefficients are determined by 


= poboi Zo(n) — 


(i/mk)eKy(ik) (Ynp/po)fe” (52 


hy(n)-Lip(n) 


h(n) Lng(n) = qobo) —d(n)b(n)/bo- cos ¢-g,(n) — 
(i/mk)eKy(tk) (¥ng/Qo)je” (53) 
ay(n)-Lap(n) = fo, (1/2 — t/k) (bo/s) sin ¢- 


cos ¢:g,'(n) — 
(i/mk)eKy(ik) (b0/b) (Yap/poje” (54) 


l 


The motion of the oscillating wing follows the La- 
grange equations of motion 
jot) 
1 ou I 


jor | 


1 du f 


dt { Ou) F,,, u= p,q 


(60 


where p and g are the general coordinates, 7 and U the 
kinetic and potential energies and F, and F, the general 
ized forces of the system. 


Using the definitions 


(a) M(n, . wing mass per unit span, measured 
parallel to elastic axis 

(b) S(n, . Wing static unbalance per unit span 
about elastic axis 

(c) L[(n, . wing moment of inertia per unit span 
about elastic axis 

(d) w, ... fundamental uncoupled bending frequency 
(bending about axis A,B, 

(e) w,... fundamental uncoupled torsion frequency 
(twist about elastic axis) 

. structural damping coefficients asso- 


(q) 


(f) g”’, 
ciated with the two degrees of freedom 


one has the fundamental relations 


T = (1/2) | p(t)}2? }sbo2/ cos ¢} [ M(n,) {fo(no)}? dn, + {p(t)} {q(t)}-{sb)/cos g} X 


/70 


1 l 
S(n, o(ne)falnddn, + (1/2) 1 q(t)} 2! s/cos yt { I(n,) lfe(ne)} dn, (61) 
J0 J0 


71 


U = (1 + ig”) (w,"/2) ) p(t)}? } by?s /cos ¢} | 


/7 0 


M(n,) iSp(ne) | P dn, oF 


7) 
| ling) {hiln,){? dng (62) 


d d fee 
Fio=s h,(n)L(n) + a,(n) A ( dn (63) 
J 0 \dp dp 
wa d 
Fi =s h,(n)L(n) + a(n) A\1(n) D dn (64) 
J0 dq dq 
But 
d\h,(n)} dp == by Zp(n (65) 
dy h(n t dq = —d(n)b(n) cos ¢-g,(n) (66) 
d\ a,(n)| ‘dp = (bo/s) sin ¢-cos ¢ g,'(n) (67) 
d} a(n)} /dq = cos ¢-g(n) (68) 
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and hence 


*1 
F, = rpw*shy | Lyn) (O*hpLny + WhgLng + BapLay + Balog dn + 
7) 
pw” by sin g cos ¢ [ £p'(n) Oh, My, + b%hy Mag + btayMay + btagM Qo} dn (69 
7) 

* 
Fy = — pw" s- cos e| An) b(n) gq(n) (OAL ny + Vhgling + BayLay + Baglag} dn + 

0 


| 
Tpw-sS* COS ¢ | £y(n) Oh Mny + b%hy Ming + btay May + bta Magi dn (70 
7/70 
After expressions for the generalized forces are combined with the terms derived from the kinetic and potential 
energies, a homogeneous system of two linear equations in the two unknowns p and g results. These may be 
written 


A-p+ B-¢ = 0 (71 
D-p+E-q=0 (72 
where 
*! ad | 
A -cos g-mpby? = ,1 — (1 + ig'”)§ (w/w)? | M(n)fo(n) f(a) dn + (b by)? Lay(n)go(n)go(n) dn 4 
/7 0 7 


7] | 
sin gy: COS ¢: } (b/ by)? (b/s) Lap (n)go(n)gp'(n) dn + sin ¢ ‘cos yg: } (b by)? (6/8) Map (n) p(n) gp"(n) dn 4 


7 WV 


sin”? ¢- cos? ¢ (b by)7(b's)? MW 0 (n) g,'(n)g,'(n) dn (73 
0 


>] | 
B- cos g:rph> = | S(m)fi(mfyln) dn — cos ¢: (b by) 3d (ng) Lig(n) gp(n g,(n) dn + 


SY 
0 / 0 


7! */ 
cos ¢: | (/by)* Lega) Sp (man) dn — sin g- cos? ¢° (b by) ®(b/s)d(n) Mig ()go(m)£p'(n) dy 
0 


e 7 0 


bd 
sill ¢- COS" ¢- (b by)3(b 8) Mao(n) go (ngp'(n) dn (7A 
Zi 
+! *] 
D+ cos g- mph? = } S(n)filmfiln) dn — cos ¢- | (b/by)8d(9) Liy(m)Zo(n) gan) dn 4 
0 { 


| | 
cos ¢: } (b/by)® May (n)gp(n)g,(n) dn — sin g- cos® ¢: (bby) 3(b 's\d(n) Lay (1) Zo (n)Ep'(n) dyn + 
0 


Sin ©: COS” ¢- b/b,)3(b/s)M (n)g, (ng ‘(n dn eer 
0) 


ad | >| 
E- cos g:- mpby' = ,1 — (1 + ig’) (w/w)? L(n)fi(n)filn) dn — cos® ¢ (b by) 4d(q) Lag a) ga(n) gala) 4 
7X v7 
>| >) 
cos” ¢ (b/ by) 4} d(n)t? Ling (n)gq(n)gy(n) dn — cos* ¢ (b/by)4) d(n)t Miag(n)g,(n)g,(n) dn 4 
7 7/7 


>] 


cos’ ¢: (bby) AM (mM go(n)go(n) dn (76 


The system (71), (72) can have a nontrivial solution Eq. 79 is a function of the wing root reduced frequency 


only if its determinant, the flutter frequency deter- ky) = dyw/V and may now be solved assigning different 
minant, vanishes: numerical values to kj). In this manner the root 2 cor 
Ae on BP nD . (77 responding to the preassumed value of g may be found. 

Eq. 78 then yields the flutter frequency w and the flutter 

One may put, within certain limits of accuracy, g'”) = speed V is found from I’ = ko/wb). This speed, of 
g? = g, introduce the symbol course, refers to incompressible flow and would have 


to be corrected by the Prandtl-Glauert Rule, V (compr. 
Q = (1 + ig) (a, w)* (7S) — (] 


— (Mach)?) “‘- V (incompr.), if it is considered 
; . — ; necessary to account for compressibility effects. Eq. 

and write the flutter frequency equation in the form fs : : ‘ ‘am ; 

79 has two different roots for each specific value of ko 


2 > a - Gok. 
2 + P-2 + G=0 (79) (Continued on page 834) 
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A Method for Analvzing the Heat Insulating 
Properties of the Laminar Compressible 
Boundary Laver’ 


PAUL A. LIBBY? ann ADRIAN PALLONE?{ 
Polytechnic Institute of Brooklyn 


SUMMARY 


In some cooling problems associated with high energy flows 
it may be convenient to localize strongly the cooling, as for ex 
ample by injecting a coolant through an upstream porous strip, 
and to depend on the insulating properties of the boundary layer 
to reduce, or to eliminate completely the need for further cooling 
on the surface downstream of the highly cooled section. This 
upstream cooling technique may be of interest in connection with 
optical windows in hypersonic wind tunnels, and on radomes, 
wings, and bodies of high-speed aircraft and missiles 

In this paper a method for investigating the insulating proper 
ties of a laminar compressible boundary layer on a two-dimen 
sional surface with zero heat transfer is presented. The physical 
situation considered thus corresponds to the case in which the 
heat transfer downstream of the strongly cooled section is com 
pletely eliminated. Of practical concern is how the temper 
ature of the uncooled surface varies in the downstream direction 
from its low initial value and thus how the low energy layer 
established by the upstream cooling insulates the downstream 
surface 

The Karman integral method extended to both the momentum 
ind energy partial differential equations of the boundary layer 
has been used. The station, at which cooling and/or injection 
ceases, corresponds to a discontinuity in boundary conditions 
ind thus in solutions. At this point the flux of mass, momentum, 
ind energy within the boundary layer has been made continuous 
by the introduction of three additional parameters in the velocity 
and stagnation enthalpy profiles. Thus the velocity and stag- 
nation enthalpy profiles have both been taken as sixth degree 
polynomials. The resulting two integral-differential equations 
are then solved for two unknown functions of the distance along 
These two functions are related to the boundary-layer 
Initial conditions corre 


the wall. 
thickness and to the wall temperature 
sponding to a given initial wall temperature and an initial bound 
ary-layer thickness are prescribed. Exact closed-form solutions 
for the case of zero axial pressure gradient are obtained. For 
flows with significant pressure gradients, numerical solutions 
are required in general. Several numerical examples of practical 
interest are presented 


SYMBOLS 


F,, 2, = integrals defined after Eq. (14 
G = function of 1) and A%, [ef. Eq. (34 
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and Applied Mechanics. 


H = stagnation enthalpy (u?/2) + c,7 

M = Mach Number 

p = static pressure 

R = gas constant 

S Sutherland constant, 216°R. for air 

T = an average temperature 

1 = absolute temperature 

i = transformed variable [cf. Eq. (11 

u = velocity component in x direction 

l = potential flow velocity parallel to the surface 

velocity component in the y direction 

iH" = ratio of stagnation enthalpy at the wall and in po 
tential flow 

x = coordinate along the surface 

y coordinate normal to the surface 

Fo = pressure gradient parameter;  [cf. Eqs. (19) and 
(43 


= ratio of specific heats 
= boundary-layer thickness 


0 

nN = nondimensional boundary-layer thickness ratio 
(A A/Ao) 

mM = coefficient of viscosity 

é = nondimensional space coordinate (£ pol ox/p 

p = mass density 

A = nondimensional boundary-layer thickness parameter 
(A = pl po lb? /pipo) 

T t/6 

y¥ = (p,7,,/polo)A 


Subscripts 
= value at the surface, y = 0, for any value of x 
value at the outer edge of the boundary layer 
( = valueatx = 0 
( = value in the x-/ or £-7 plane 
( = stagnation values 


INTRODUCTION 


. OF THE MOST important current problems in the 
aeronautical engineering field are those connected 
with preventing critical components of aircraft, pro 
pulsion units, and test equipment from being heated to 
temperatures greater than those permitted by struc 
While im 


proved materials can give some temporary relief, it 


tural, functional, or life considerations. 


appears that only by appropriate cooling can most of 
these heat transfer problems be solved. 

In many practical cases involving cooling, it may be 
convenient or even necessary from a functional or eco- 
nomic point of view to localize the cooling in a re- 
stricted region. For example, the sweat cooling tech 
nique (cf. for example references 1 and 2), which is of 
considerable interest, requires the use of porous mate- 
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rials which are expensive and which have, at the pres- 


ent, reduced strength properties. It may thus be 
desirable to localize their use. 


provided by the radome on high-speed missiles; the 


Another example is 


radiation requirements make cooling difficult and may 
even dictate the use of special materials, which have 
poor strength properties at elevated temperatures and 
which, therefore, increase the complexity of the heat 
transfer problem. Finally, the optical windows in 
hypersonic wind tunnels operating at elevated stagna- 
tion temperatures are difficult to cool without seri- 
ously affecting their use in connection with schlieren 
and shadowgraph systems. 

A method of preventing excessive temperatures from 
being reached or of realizing the advantages of localized 
cooling may be to cool the upstream portions of the 
critical section and to depend on the insulating prop- 
erties of the boundary layer to protect the downstream 
surface from the high-energy outside flow. The up- 
stream cooling results in the establishment of a con- 
tinuously replaced, low-energy layer of gas within the 
boundary layer and between the exposed surface and 
the high-energy flow. It might be mentioned that 
film cooling, which has been used in protecting rocket 
motor liners for some time, depends in a similar manner 
on the insulating and heat absorbing properties of an 
evaporating liquid layer. The efficacy of upstream 
cooling for any particular application is determined by 
whether the surface temperature does or does not in- 
crease from its initial low value at the upstream edge 
to above an acceptable value within the length of the 
critical section with reduced or no heat transfer. 

In this paper a method is presented for obtaining an 
indication of this variation of wall temperature when 
the boundary layer is laminar and two-dimensional 
and when there is no heat transfer downstream of the 
highly cooled section. The flow is assumed to be two- 
dimensional, steady and to be describable by the 
Prandtl boundary-layer equations. The fluid is as- 
sumed to be a perfect gas with a constant Prandtl 
Number of unity and with a linear viscosity-temper- 
The heat conductivity of the surface 
It is pointed out 


ature relation. 
iaterial is assumed negligibly small. 
that these assumptions have all been used extensively 
in boundary-layer theory in the past and have been 
considered satisfactory for many engineering analyses. 

The investigation is carried out along the following 
lines: It is assumed that the boundary-layer char- 
acteristics at the downstream edge of the cooled por- 
tion of the surface are known.* This edge represents 
a discontinuity in the solution to the boundary-layer 
equations since the boundary conditions are discon- 
tinuous there and, therefore, requires special treatment 


* The methods of analysis of references 2, 3, and 4+ may be of 
value in determining these characteristics. For completeness, 
in the Appendix, the case of cooling by injection over a constant 
pressure surface has been treated in terms of the variables t ved 


here. 
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as will be discussed presently. The momentum and 
energy equations of the laminar boundary layer are 
formulated as two integral-differential equations for 
the case of axial pressure gradient with no heat transfer 
After the method of von Karman profiles for the local 
velocity and local stagnation enthalpy in terms of the 
distance normal to the wall are assumed. These pro 
files satisfy certain boundary conditions at the wall 
and at the outer edge of the boundary layer. In addi- 
tion three parameters are introduced into these profiles 
and are determined by conditions at the leading edge 
discontinuity where the flux of the mass, momentum, 
and energy within the boundary layer are made con 
tinuous. It is shown that the discontinuities in the 
local velocity and local stagnation enthalpy profiles 
are comparatively small when the overall conditions 
are imposed. The downstream length in which an ad 
justment of the boundary-layer flow takes place can, 
therefore, be expected to be small compared to the 
length of interest in examining the efficacy of the in 
sulating ability of the boundary layer. 

With the assumption of these profiles the two in 
tegral-differential equations can be converted to or 
dinary differential equations which give essentially the 
downstream variation of the boundary-layer thickness 
and the wall stagnation enthalpy from some prescribed 
initial values. For the case of zero axial pressure 
gradient these two equations can be solved in closed 
general gradient 


form. However, for the 


case a numerical solution of the differential equations 


pressure 


must be employed. 

At first thought it may be disturbing that the results 
of this analysis indicate zero heat transfer on a surface 
with variable wall temperature. Ordinarily for Prandtl 
Number of unity the surface temperature for zero heat 
transfer is constant and equal to the stagnation tem 
perature of the flow. However, most boundary-layer 
analyses assume an isoenergetic flow at the leading 
edge whereas here the effect of the upstream cooling is 
to provide an energy variation at the leading edge of 
the surface under consideration. It is of interest to 
point out that Levy® has obtained from an exact nu 
merical integration of the constant gas properties, 
laminar flow equations, zero heat transfer, variable wall 
temperature solutions. As pointed out by Eckert' 
these solutions involve an infinite temperature difference 
at the leading edge of the surface and thus mathemati 
cally although not physically are analogous to those ob 
tained here. 

The authors are pleased to acknowledge that Professor 
Antonio Ferri suggested in private conversations the 
upstream cooling technique analyzed here and has been 
a source of helpful discussion throughout the research. 


Basic EQUATIONS 
The following equations describe the two-dimensional, 
steady laminar boundary-layer flow over a slightly 
curved surface of a perfect gas with a Prandtl Number 
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METHOD FOR ANALYZING 


of unity and with constant coefficients of specific heat 


Equation of motion in x direction: 
puu, + piu, = —pr t+ (uu,), (1) 
Equation of motion in y direction: 
Py = QO: —Dr = pUi(dU, ‘dx) (2 


Equation of mass conservation: 
(pu) + (pv), = 0 (3) 
Equation of energy conservation: 


pull, + pvll, = (uH,), (4) 


Equation of state: p = pRT (5) 
Viscosity-temperature relation: uw = Cyuo(7T/7; (6) 
where 

p = mass density of gas 

u = velocity component in x direction 

v = velocity component in y direction 

p = static pressure 

a = coefficient of viscosity 

H = (u*/2) + c,T, local stagnation enthalpy 


R = gas constant 

7 = static temperature 

C = (T,,/To)” (To + S)/(T« + S) 

S = Sutherland constant, 216°R. for air 
where the subscript 0 denotes conditions outside the 
boundary layer at x = 0, the subscript 1 conditions out- 
side the boundary layer at any value of x, and the sub- 
script w conditions at the surface (y = 0), and where 
is some average wall temperature used to modify the 
linear viscosity-temperature relation so that the exact 
Sutherland relation is closely approximated at the 
average wall temperature in the region of x of interest. 
This modified linear viscosity-temperature relation was 
first suggested by Chapman’ and has since been used 
extensively. 

Eqs. (1), (8), (4), and (5) with the accompanying 
equations and definitions may be considered to be four 
equations giving p, u, v, and // as functions of x and y for 
certain boundary and initial conditions. As is usual in 
boundary-layer analyses, the pressure and/or the poten- 
tial velocity distribution in the x direction is assumed 
known from potential flow considerations and the re- 
quired fluid constants R, S, and c, specified. 

For the investigation of the insulating properties of 
the laminar boundary layer given here, the following 
boundary and initial conditions are prescribed: The 
origin of the x axis is located at the upstream edge of 
the critical, uncooled section. At x» = 0, in particular, 
the boundary layer thickness and wall temperature 
are known while more generally the « velocity and 
stagnation enthalpy distributions determined from an 
analysis of the upstream boundary layer are pre- 
scribed. Along the positive x axis the velocity com- 
ponents « and v are as usual assumed zero, and the 
heat transfer is assumed zero corresponding to the as- 
sumption that the critical section of the surface (x > 0) 
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is thermally insulated. At the outer edge of the 
boundary layer or y ~ ~, the u velocity takes on pre- 
scribed values denoted by U, the density by p; and the 
stagnation enthalpy the constant value /. It is re- 
quired to find in particular for this investigation the 
distribution of stagnation enthalpy along the surface 
(x > 0), since from this the distribution of wall tem- 
perature can immediately be found. 

A convenient way to solve approximately the above 
formulated problem is to use the integral methods due 
to von Karman. Eqs. (1) and (4) are multiplied by dy 
integrated from y = 0 to y = 6(x) where 6 is the bound- 
ary layer thickness, an unknown function of x, and in 
this analysis is the value of y in which both the wu ve- 
locity and the stagnation enthalpy // are equal to their 


potential flow values.* Applying the boundary con 


ditions, 
u=-v=0'| 

at y = 0, (7) 
al H, = 0 f 
and at y = 6(x), 

se = U;, t = 0 | 

, (8) 
HW = #H,, H, = 0} 


and using Eq. (3), one can convert Eqs. (1) and (4) into 
the following integral-differential equations: 


6 
(d/dx) yal’ [ (p/p) (u/U;) fl — (u U1)] dyt + 


oi Ui(dU/dx) (p/p1) [(pi/p) — (u U;)] x 
J0 


dy = (uty). (9) 


/ 
(ed 


(d dx) yell (p pi) (u U;) [1 = (H H,)) dys = (0 
/70 


{ 
(10) 


Eqs. (9) and (10) can be handled more conveniently if 
the Dorodnitzyn transformation is introduced. A new 


variable ¢ is defined such that for a given value of x 
dt = (p/p) dy (11) 


This implies the introduction of a new boundary-layer 
thickness in the x-/ plane defined most conveniently by 
the inverse transformation 


§ = (p/p) dt (12) 
0 


If Eqs. (11) and (12) are applied, Eqs. (9) and (10) can 

* For the case of isoenergetic flow and Prandtl Number of 
unity the frequently distinguished thermal and velocity thick- 
nesses are exactly equal. For the flow with variable stagnation 
enthalpy treated here the two thicknesses cannot be expected 
to be greatly different Moreover, in the integral method there 
appear to be good mathematical reasons for adopting the point 
of view of one thickness within which all fluid properties change 
from their wall values to their local potential flow values. (Cf 


reference 4 
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be nondimensionalized and put in the following form :* 


(F,/2)A’ + F,Afin (FiG)]’ + 
FA{In (04/U)]' = [(a/U),]e (13) 


(GA “F;)’ = © (14) 
where 
*] 
i, = / (u/U,) [1 — (u/U))] dr 
JZ 0 
*] 
F, = l a (¥ — 1) (M,? 2)] } (77 /T,) — 
/7 
(u/U,)*] d1 
*] 
i; = | (u/U,) [1 — UT/M,)] dr 
/7 J 
NO = pil pol v6 /uim 
G? = (pi/po) (U1/ U5) (ur /mo) 
r = t/6, 
y = ratio of specific heats 


JJ; = potential flow Mach Number 


)’ denotes differentiation with respect to 
Eqs. (15) and (14) are in a form for 


and where ( 


E = pol w/w. 
approximate solution following the method of von Kar- 
man. Polynomials in 7 are assumed for u/l, and 
H/H,; the coetticients of these polynomials are selected 
so that the boundary conditions in the € — 7 plane cor- 
responding to those given by Eqs. (7) and (S) in the 
x — y plane are satisfied for any value of 6, and in this 
application of the method for any value of the wall 
temperature. Additional accuracy can be realized if 
the polynomials are taken of sufficiently high order so 
that they satisfy several additional boundary condi- 
tions corresponding to those which an exact solution 
of the partial differential equations would automatt- 
cally satisfy. Moreover, three additional parameters 
that are independent of x» and which are determined 
from conditions at the discontinuity existing at v = 0 
are introduced. t 

Sufficient accuracy for the purpose of the present 
investigation can be realized if sixth degree polynomials 
in the velocity and stagnation profiles are assumed. 
It may be noted that since three of the parameters are 
determined from leading edge conditions this corre 
sponds to assuming fourth and fifth degree polynomials 
for the velocity and stagnation enthalpy profiles respec- 
tively; such profiles have been found satisfactory for 
heat transfer and skin friction determination (cf. refer- 


ences 4 andS). Thus 


and 


* Although there are many possible choices of variables and 
methods of nondimensionalizing, the ones selected here seem 
especially convenient 

+ Professor Martin H. Bloom suggested the possibility of these 
additional parameters to the first author in a private conversa 


tion. 
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i 
A/ii, = > bjt (16 
7 0 
The boundary conditions on u/l, and ///Z/, in the &-r 
plane corresponding to those in the x-y plane given 
by Eqs. (7) and (S) are: 


t 0 u/ ly = 0 17 
a T= J ( ) 
(1/11), = 0 ‘ 
ind at as oh - a IS 
‘ atr = 1, : (1S) 
(iu l’)), == (II TT;) = 0 
Additional boundary conditions at y = 0 can be ob 
tained from Eqs. (1) and (4) with w = v = 0; they are 
—PD, = pil (dl dx) = —-(“wH, 
(ulT,), = ( 


which in the &-7 plane and in terms of A and € are at 


7 = U, 


(u/U,).-, = —[1 + (y — 1) (V,2/2)]WA[In 
(U4/U))’ 
(19) 
= —22 
(1/1h),, = 0 
where IV’ = //,,///,, at this point in the analysis an un 


known function of &, and where Z is evidently defined. 


Similarly, additional boundary conditions at y = 
6(x) can be obtained from Eqs. (1) and (4) with a 


IT, = 0; they are in the é-7 plane: 
at 7; = I, (u/U;).. = (47/f2,), = 0 (20) 
At the leading edge corresponding to x = O a discon 


tinuity in solution exists because the boundary condi 
tions are not continuous. Indeed, it is possible that in 
the immediate neighborhood thereof the boundary 
layer equations themselves may not be valid since the 
discontinuity corresponds to an abrupt change of the 
flow variables in the x direction; fundamentally the 
boundary-layer hypothesis precludes such changes. 
However, here the discontinuity has been handled in 
the customary manner; namely, the equations have been 
taken to be valid on either side of the discontinuity and 
finite conditions corresponding to conservation of 
mass, momentum, and energy across the discontinuity 
have been imposed. This approach can be expected 
to be valid provided the discontinuities of the local 
velocity and stagnation enthalpy within the boundary 
layer are relatively small since then it could be expected 
that a readjustment zone on the order of only several 
boundary-layer thicknesses would be required for the 
description of the boundary-layer properties presented 
here to be valid. It might be mentioned that in the 
past such discontinuities have been treated by making 
either the boundary-layer thickness or the momentum 
thickness continuous; both of these conditions are con- 
tained in the more complete treatment used here. 

The conditions at « = O will now be discussed ex- 
plicitly; the notation superscript (+) and (—) will be 
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used to denote values on the sides x > 0 and x < 0 of the 
discontinuity respectively. The finite potential flow 
quantities are taken to be continuous although dl\/dx 
may be discontinuous. Moreover, the boundary-layer 
thickness in the x-y plane is taken to be continuous; to 
do otherwise would be physically unappealing. Along 
with the other conditions of continuity of mass, mo- 
mentum, and energy within the boundary layer im- 


posed, the boundary-layer thickness in the x-/ plane 


and thus A also are made continuous. Thus 
\ = Att (21 
F; = F, (29 
F, nf (23 
} = Ff, (24 


where /; = | (u/U) dr. 


Eq. (21) has been discussed 


above; Eqs. (22), (23), and (24) result from making the 
mass, momentum, and energy, respectively, continuous 
and from Eq. (21). 

It is pointed out that Z defined by Eq. (19) is con- 


tinuous if dl’,/dx is continuous at x = O but with in- 
jection upstream this is not the same as making (u/ l)),, 
continuous. Actually, (#/l4),, is discontinuous at 


x = Oif injection exists upstream, or if a discontinuity 
in dl',/dx exists at x = 0. It is pointed out that this 
latter can occur frequently in supersonic flow since 
only a discontinuity in surface curvature is required 
to lead to a discontinuity in d1/,/dé. 

With Eqs. (17) 
the 14 coefficients of the velocity and stagnation en- 


through (20) and (22) through (24 


thalpy profiles can be evaluated in terms of Z, Il’, and 
of the jumps in the three parameters determined by the 
discontinuity conditions. These have been taken to 
be ds, ds, and bs, parameters which are zero for x < 0; 
thus the jumps are represented solely by a;‘~*’, a¢‘*', and 
bs‘*’, and can be expressed in terms of the parameter; 
ast = —Z(0), bo, bo, b . and 


ay » ae 2 
, which are determinable from initial con- 


b,‘+ = db 


ditions. Thus 


while as‘ *’ is determined from the equation 


deg. * (3a; . 1540) + (d»'* 2772) - 
(ag°*? /12012)| = (8/231) Ia; — a,‘*| + 
(1/277) [ae — ay‘*|] + (23/4620) X 
Cas" — @ Qe + (52/3465) |(a;°*)? — 
(ay 2] + (1/2310) [(ao'*)? — (ao‘—’)* (26) 
where 
Qy = 2 — (1/6) [a + ast] + (ag°*?/14) (27) 


Finally b.‘*? is determined from the equation 
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be* [(11/3080) + (3a, 30SO) + (ae'* /5544 
(a6°*?/12012)] = (1 — Wo) [(811/4620) X 
(a, — a,‘*’) + (281/27720) (a: — a,‘*’) + 
(llag‘*’/30S80)) — bd, (131/4620) + 
(52 ay 3465) + (23a. 9240 
by (181/27720) + (23a, 9240) + 
(Qs 2310 (28 
In summary then, at the discontinuity (£ = 0) one has 
at his disposal from the analysis of the boundary layer 


upstream (x < 0), a:\~’, a , b:'—, bo’. Moreover, 


if A and IV are made continuous across the discontinuity, 


Z(0) = —as and b, = by'*+’ = W, are known. 
From Eqs. (25) through (28) the constants a;‘*’, 


dst, and bg are computed. Thus the velocity 
and stagnation enthalpy profiles for x > O are specified 
except for Z and II’, which are determined by the solu 
tions of Eqs. (13) and (14 

Explicitly, the velocity and stagnation enthalpy 


profiles are: 


u/U, = [2r — 27° + 7] +4 
Z\ (7/3) — 7? + 7° — (7'7/3)) Ft 
as —(7r/3) + 2r* — (8r*/3) + 7°] 4+ 
as r+or® — ori t+ r' (29 


and 


T/T, = 1 — (1 — W)/1 — 107° + 1574 — Gr*! 4 


7 


Henceforth, the superscript (+) on a5\*', a 

will be dropped since it is understood that a 
It will be noted that Eqs. (29) and 
30) reduce to the customary profiles of fourth and 
a by 0 (cf. refer 


! == b, sp. 


fifth degree, respectively, if a 
ence 4 
The F integrals become with these profiles 


F, = (37/315) + (2/124740) X 
-264 + 2llas + 603a6) 4 
a5/415S0) (352 — 443a6) + 
(203a6/6930) — (27/2268) — 
D2 a5°/31185 1417 ag? /S3160 31 
F, = }(1 + (y — 1) (.M)2/2)]/27720; 
}11570 — 13860(1 — W) + (1159 — 259a6)a5 4 
3783a5 — (395 b¢/2 16a52 — 472 as? + 


(—521 + 14745 + 46a;)Z 110 9)Z7; 9 (32 


F; = (1 — W) [(17/70) — (311as5/13860) 4 
(31Z/2520) — (13la¢/1S84S8)] + b¢6/(17/3080 
(a;/3080) — (127a6/120120) + (2/6930 33 


The solution of Eqs. (13) and (14) may now be con- 
veniently discussed. It will be assumed that the po- 
tential flow quantities (denoted here by the subscript 
1) are expressed in terms of .\/, and that .J/,; is a known 


function of £. There is no loss in generality in making 
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this assumption, since if, for example, the potential 
flow characteristics are known in terms of U/l as a 
function of ~, U/l can be expressed in terms of 1, 
and the discussion presented here then applied. Apply- 
ing Eq. (6) and the usual isentropic relations in the 
potential flow, one finds 
G? = C(M,/M) {1 + (vy — 1) (AL2/2)] + |} 

(1 + (y — 1) (M,2/2)]}{ 87-9207! 
U,/Uo = (Mi/Mo) {01 + (vy — 1) (Mo?2/2))] + 


7 r (34) 
(1 + (y — 1) (M,2/2)]}'” 


| 


where J, is the potential flow Mach Number at & = 0. 

Eq. (14) can be integrated directly; if the initial 
conditions that A = Ay and W = JW, at — = O are im- 
posed, one obtains the following algebraic relation 
which relates implicitly A and W-: 


GoAo (F3)o => GA . F, (35) 


Note that Gp = C and thus that Eq. (35) is independent 
of C. Finally Eq. (13) can be written symbolically as 


A’ = f(A, W, Z, My, as, as, be) (36) 


With Eqs. (35) and (19) relating W and Z, respectively, 
to A, with 14, = A/,(&) known, and with as, ds, and bg 
determined by Eqs. (25) through (28), Eq. (36) can be 
integrated numerically to give A = A(é) subject to the 
initial condition A = Ay when = 0. The wall tem- 
perature which is of prime interest here can be found 
from the equation 

T/T, = W/Wo (37) 
where 7\,, is the initial wall temperature and where 
W is found as a function of € from Eq. (385) and from 
the solution to Eq. (36). 

The solution of Eq. (36) for the case of zero pressure 
gradient is of interest for the insulation of the walls of 
hypersonic wind tunnel test sections and on wedge- 
shaped airfoil profiles and will be discussed first. 


ZERO PRESSURE GRADIENT 


For the case with Z = 0, Eq. (86) becomes explicitly 


(F,/2)A’ = a'?) = 2 — (a;/3) — ae (38) 
as may be seen from Eqs. (13) and (29). Eq. (88) 
can be integrated directly; the solution with A = Ag 


até = Ois: 
A/Ay = 1+ (2/F,) [2 — (as5/3) — ag] (E/Ao) (39) 


Thus for this case 7,, = 7),(€) can be found readily 
since from Eq. (35), IV = IW(&) can be found upon sub- 
stitution of Eq. (39) and finally 7), is determined di- 


> 


rectly from Eq. (37). 


FLOWS WITH AXIAL PRESSURE GRADIENTS 
In the general case, in which 1/, is a prescribed func- 
tion of £, solution of Eq. (36) involves a numerical 
integration. However, the numerical work is reduced 


if a new dependent variable 
A = A/Ao (40 


is introduced and if JJ, is considered the independent 
variable. Eq. (13) becomes with these changes and 
with the right-hand side expressed in terms of Z by 
Eq. (29) 


(F, 2) (dr d\,) +X , F, (d dAI,) In (F,\G)] + 
F.[(d/dMy) In (U;/U»)]} = 
{2 + Z/3 - (a5 3) — ag dM, d(é Ao) ] ' (41 


If now the IV variable appearing in F, and implicitly 


given as a function of A by Eq. (35) is eliminated from 
Eq. (41), the final equation for numerical integration is 
obtained; it is 
(F\/2) (dd/dM) + (A/G) [m + 2] — 
(A ANG) 32 = 2 —- (a5 3) — a6] x 
aM,/d(E/Ao)|—! (42 


where 


m = m(M,) = (1 — y4q?)/2[1 + (y — 1) (M2/2)] 


2 = 2(Z) = (1/27720) [6952 + a;(1159 — 295a-6) + 
(3783a5) — (395b6/2) — (46a52) — (472a6562) + 
Z(—521 + 147a5 + 46a;)] — (22/2268) + 

(b6/3) [(17/3080) — (as/3080) — (127a5, + 
120120) + (2/6930)] [(17/70) — (8llas + 


13860) + (312/2520) — (131la¢/1848)]-! 


Zo = 3(Z) = [(F3)0o/3] [(17/70) — (311a5/13860) + 


(31Z/2520) — (131la5/1848) 


and where G = G(\/)) is given explicitly by Eq. (34). 

Eq. (42) can be integrated numerically by any of the 
standard techniques for first order, ordinary differ- 
ential equations, as for example Adam's method 
(reference 9), from the initial condition \ = 1 at 14, = 
My. The only auxiliary equation required is 


Z = (Wd/2M1,) [dA /d(E/Ao) | (43) 


J 


along with [d.1/,/d(é/Ay)| which must be determined 
from M, = M,(&) and Ao. With A = A(M) = A(é) 
given by Eq. (42) W = IV(é) and 7, = T7,,(&) are found 


from Eqs. (35) and (37), respectively. 


NUMERICAL EXAMPLES 


In this section several numerical examples based on 
the above analysis have been carried out. 

Example (1): Optical windows in hypersonic wind 
tunnel walls. As a numerical example of the zero 
axial pressure gradient analysis the protection of wind 
tunnel optical windows by coolant injected upstream 
according to the sweat cooling technique is considered. 
The following conditions are assumed: stagnation 
pressure, Ps = 500 psia, stagnation temperature, 7, = 
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2000°R., My = 4.0, W) = 0.30, the tunnel half-height 
h = 1 ft., and the coolant temperature 7, = 500°R. 

The wall temperature distribution was computed 
for the two methods of injection as indicated by the 
discussion of Eq. (A-14) in the Appendix: namely, 
constant injection with variable surface temperature, 
and constant surface temperature with nonuniform 
injection. In both methods the injection was dis- 
tributed over a 0.5-ft. length upstream of the region to 
be protected and the mass flow was adjusted so that 
I, ~ 0.3 in the first method and so that W = 0.3 for 

-&; < & < Oin the second method. It is pointed out 
that the values of a5, ds, and bs were 9.69, —2.38, — 14.89, 
respectively, for constant injection and 9.04, —3.35, 
and —8.29 for constant wall temperature. 

In Fig. 1 the variation of surface temperature with x 
is shown. Also shown are the total amounts of coolant 
mass flow in terms of the mass flowing per unit length 
per unit time in the nozzle half-height. The significant 
insulating effect of small amounts of injection is evi- 
dent. It is pointed out that with p,v,/pol’) on the 
order of 10~* throughout most of the length over which 
the injection is distributed and with the effect of cool- 
ing considered, laminar layers could be expected 
provided of course that the boundary layer is thin 
and laminar at é = —€&;. It will be noted that the 
boundary-layer thickness at this point has implicitly 
been assumed negligibly small compared to the bound- 
Undoubtedly, it would 


ary-layer thickness at — = 0. 
be of interest to study the insulating ability of the tur- 
bulent boundary-layer since circumstances wherein the 
boundary layer is already turbulent at & = —&, or be- 
comes turbulent due to large rates of injection will 
arise. However, at the present time not enough is 
known abgut turbulent boundary layers with injection 
to make this investigation possible. It is furthermore 
pointed out that the singularities at & = —E£, for both 
methods of injection, the first yielding a singularity 
in that IV = 1 and the second in that p,v,,/pol> > ©, 
are typical of the singularities encountered at the 
leading edge in boundary-layer theory. 

It is of interest to compare the velocity and stagna- 
tion enthalpy profiles on either side of the discontinuity 
at £ = 0. This has been done for the first method in 
injection; that is, with p,%/pol) constant; the re- 
sults are shown in Figs. 2and3. It will be seen that 
the velocity profiles on either side are almost identical 
indicating an insensitivity in the £7 plane to the 
change in heat transfer at the wall but that, as might 
be expected, the stagnation enthalpy profile is modi- 
fied somewhat by the discontinuity in wall heat trans- 
fer. Thus an adjustment zone downstream of & = 0 
wherein the flow description given here may not be 
quite accurate for determination of local quantities 
away from the boundaries, y = 0 and y = 4, is indi- 
cated. However, since the correct boundary condi- 
tions are prescribed it is believed that the existence of 
this zone does not alter significantly the distribution of 
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wall temperature which is based on the integrated 
values of the profiles and on the boundary conditions, 
both of which are correct. 

Example (2). Biconvex Airfoil: The insulating 
ability of the laminar boundary layer over a biconvex 
airfoil was studied by numerically integrating Eq. (42). 
The airfoil was taken as 4 per cent thick with a chord 
of + ft. It was assumed to be in steady flight at a 
Mach Number of 6 at 50,000 ft. The local Mach Num- 
ber distribution over the surfaces was determined by 
the shock wave expansion method neglecting reflected 
waves (cf. reference 10). The leading 10 per cent of 
each surface was considered porous and cooling air at 
SOO0°R. assumed available. The coolant was injected, 
at a nonuniform rate, so as to give a constant wall tem- 
perature of 0.37°, along the leading 10 per cent of each 
surface. The average rate of coolant injection was 
1.1 X& 10-* poU. 
the method described in the Appendix. 


The initial conditions were found by 
The param- 
eters d5, a6, bg were 13.77, —3.695, and —7.39, respec 
tively. 

Fig. 4 shows the distribution of surface temper- 
ature, in terms of stagnation temperature, on the sur- 
face of the airfoil. Also shown in Fig. 4, for purpose 
of comparison is the temperature distribution, as com- 
puted by neglecting the pressure gradient and using 
Eqs. (38) and (39) along a flat plate subject to the free- 
stream conditions used for the airfoil. It appears from 
a comparison of these two results that the influence of 
the favorable pressure gradient with its concomitant 
expansion of the outside flow is to keep the surface tem- 
perature lower than the equivalent constant pressure 
surface. 

From Fig. + the significant insulating ability of the 
laminar layer under these conditions may be seen. It 
may be pointed out that the actual surface temperature 
may be somewhat less than those indicated by this 
analysis since even at the altitude assumed in these com 
putations radiation will remove some heat and the 
equilibrium temperature is actually somewhat less 
than the stagnation temperature because the actual 
Prandtl Number for air is less than unity. 


CONCLUSIONS 


A method for analyzing the insulating ability of a 
laminar compressible boundary layer has been de- 
veloped. The flow situation studied involves an up- 
stream section which is cooled either by convective or 
sweat cooling; downstream of this section the surface 
is uncooled, that is adiabatic. Because of the insulat- 
ing ability of the laminar layer, the wall temperature is 
less than adiabatic and varies in the downstream direc- 
tion. 

The method is based on the integral approach. At 
the discontinuity where cooling and/or injection ceases 
the conditions of continuity of mass, momentum, 
energy, and of boundary-layer thickness are imposed 
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METHOD FOR ANALYZING 


by the addition of several parameters in the velocity 
and stagnation enthalpy profiles. 

Two numerical examples indicate strong insulating 
ability of the laminar layer when cooling upstream is 
accomplished by injection. The required rates of in- 
jection are so small that laminar layers may exist despite 


the destabilizing effect of injection. 
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APPENDIX 


For completeness and convenience an analysis of the 
laminar boundary layer with injection is presented in 
terms of the variables used in the downstream (x > 0)- 
region. The injection is assumed to occur on a con 
stant pressure strip and thus the equations correspond 
ing to Eqs. (13) and (14) for the case of constant pres- 


sure but with injection are* 


(F,/2)A’ + AFi! — (pute/poUo)A = [(u/U1), |e (A-1) 
(F3/2)A’ + AF3’ — (pute/polUs)AT — W) = 
(7/ENh),\~ (A-2) 


It should be noted that (p,v,/ py ly) is the ratio of masses 
injected and flowing in the potential stream per unit 
time and per unit area and may be a function of x. 

If fourth and fifth degree polynomials for the velocity 
and stagnation enthalpy profiles are assumed then 


4 


ef; = + 3 Q;T’ 


0 


(A-3) 


* Reference 2 contains an analysis of the case with axial pres- 


sure gradient, variable injection and constant wall temperature 


but the method differs from that used here in some details. 
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and 


H/H, = > byt A-4 


, and ///J/, in the &-r 


The boundary conditions on u/l 


plane are at 7 = 0, 


it U, = \ > 
I] IT, = W ; ss 
and at 7 = l, 
i U; = I] II, l ! AG 
= ( > 
(u/U (17/17), 0 : 


Additional boundary conditions at 7 0 can be ob 


tained from Eqs. (1) and (4) with y Q; in terms of 


A and pyt,/ pol); they are 


(u U;),, = { Pyt'y py l A’*(a/U; \-7 
(T/T), = (pute / polyA (H/T, ae 
Similarly, additional boundary conditions at y 6(4 
or tT = | can be obtained from Eqs. (1) and (4) with 

u, = 11, = 0; they are in the &-7 plane 
(u U1). (I IT, - 0 A-S 


With Eqs. (A-5)—(A-S) the coefficients for the ve 
locity and stagnation enthalpy profiles can be evalu 
ated in terms of a2, IV’, and an injection parameter. 


They are 


u/U, = (27 — 27° + 74) 4 
a —(7/3) + T r? + (7*/3 (A-9 
H/f, =1- (1-W) [i 1Or* + 157 67°] 4 
b[r — Gr? + Sr* — 37°] + (Wd,/2) X 
T i7* + OFT T (A-10 


where 


b, = I] IT; 1 
ado = y 1 + (yp 6) 
a= 2- (de 5 


y = ( Pul'y pol ‘ )A 


With Eqs. (A-9) and (A-10 
evaluated; they are 


the F integrals can be 


F, = (37/315) + (2a2/945 (do"/2268 (A-1 
F; = (1 — W) [(17/70) — (81a2/2520)] 4 
b,{ — (79/1260) + (1942/7560) } 4 


— (29/2520) + (a2/2520 A-12 


(wb,/2) 


In a practical application of sweat cooling the local 
heat transfer which is proportional to }; and the rate 
of injection p,%,/ pol’) are not independent but are re- 
lated through a local heat balance within the porous 
This introduces a coolant temperature 7’, which 


plate. 
5 indicates the 


is realistically taken constant. 
upstream injection zone being considered with the 
energy flux in an infinitesimal strip of porous plate. 
The heat balance of this element neglecting axial heat 


Fig. 
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conduction along the plate is 


(k dT dy) x = Commies a — T.) (A-13 
which in terms of the variables here becomes 
y = b,/(W — A./M,) (A-14 


where /7,=c,T,. It is pertinent to point out here that 
for constant static pressure G= 1. 

In practice with a homogeneous porous plate of con 
stant thickness and with a constant pressure difference, 
is constant, and thus Eq. (A-1I4) is satisfied 
On the other 


PwVu 
by having IW vary in the x direction. 
hand if, as may well be, it is desired to have IV’ constant, 
Eq. (A-14) is satisfied by a variable p,,v7,, which can be 
achieved in practice by either a nonhomogeneous 
porous plate or by a porous plate of variable thickness. 
These two different methods of sweat cooling require 
different solutions of Eqs. (A-1) and (A-2) and are 
discussed now. 

With p,v, constant and with Eq. (A-14+) substituted 


into Eqs. (A-1) and (A-2) one obtains 


A'?(d/dt) (FiA'”?) = 2 + W[(4 + W)/(6 + W)] (A-15) 


A?(d/dt) (FA?) = ofl — (H-/Ih)] (A-16) 


which define 
can be integrated directly; the solution with A = 0 at 


A and W’ as functions of x. Eq. (A-16) 


Cc = c iS 
— —2 
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CEMBER, 1954 


c 


FsA = Y(E + &) [1 — (1./fM))] (A-17 


Eqs. (A-12) and 


with the restriction —é; < & < 0. 
£). The integration of Eq 


(A-17) give W = W(A, 
(A-15) gives € = &(A) in the form 
(Putw/poUo)2(—E + &;) = 2.518 + 
0.888/(1 + ¥/6) + 0.3 
0.1903 In (1 + y/6)? 


12) —0.49498 tan 


0.11438y — 

336 In (1 + w/6) + 

+ 6y 4 
(A-IS 


~ 0.4475 In (y2 


(2y + (}) (12) 


Thus Eqs. (A-17) and (A-18) along with the equations 
for F; and y complete the solution for A and IW as func 
tions of x. 

For the case of IV constant, p,v, variable, the dil 
ferential equations and initial conditions are satisfied 
by taking yw and /,; constant, in which case Eqs. (A-1 
and (A-2) become 
dA/di = }2V [(4+y)/64+ WP]4+4!/F, = 

2y[1 — (1T./Ih)|/F A-19 
Eqs. (A-19) and (A-14) permit the deterinination from 
two algebraic equations of the constant values of y and 
and (/7,.//1,). With these 
values Eq. (A-19) can be integrated immediately to give 
A=A 
IN Pry 
well-known exact solutions of the laminar boundary 


b, for given values of IV 


£). It is pointed out that this solution results 


proportional to x and corresponds to the 


layer equations with such a distribution of injection. 





Development of Charts for Downwash Coefficients of Oscillating Wings of 


Finite Span and Arbitrary Plan Form 


(Continued from page 824) 


and hence two branches of roots. In most cases only 


one of the two branches has a root within the narrow 
permissible variation of g. If both branches yield 
roots of physical validity, the one with the lower flutter 
speed determines flutter conditions. 
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Three-Dimensional Laminar Boundary Layer 
with Small Cross-Flow’ 


ARTUR MAGER?# 
California Institute of Technology 


ABSTRACT 


linportant problems involving three-dimensional boundary 


layer occur in almost all internal and external aerodynamic 


configurations. For many of these, the flow outside the boundary 


layer may be resolved into a large principal component and a 


small crosswise velocity. In this paper, three-dimensional 


laminar boundary-layer flows over flat and curved surfaces art 
treated under such a simplification. For flat surfaces, the solu- 


tions demonstrate the effect of the free stream turning on the 
velocity profiles in the crosswise and primary flow directions 
When the surface curvature is large and varies so as to resemble a 
corner, the computed examples show the manner in which the 
isymmetric behavior of the boundary layer results from the 
cross-flow. The detailed examples are chosen to illustrate flows 


occurring on the casing and in the blade fillets of turboma 


chinery 
SYMBOLS 
( = cylinder radius 
Yee : Shae: = functions of z 
b = vector position 
( = curvature 
Q;t.,t, dt = 0, 1, 2, = functions of 7 
D,H constants 
fet das tis dauces Cae = functions of 7 
F = Blasius function 
Ga z 12.3 = body forces 
h(i = 1,2,3 = length stretching factors 
i = length of duct 
m,n, ?r = exponents of x 
M = Mach Number 
p . = pressure 
q = vector velocity 
q = vector velocity at y = 6 
R = cylinder path radius 
Re = Reynolds Number U*%:/p 
5 = width of duct 
u,v, Ww = velocities in the boundary 
layer 
U, V, W = velocities of external flow 
ees = Cartesian coordinates 
Fe Se 2 = coordinates in compressible 
flow 
a = yaw angle 
B = angle between tangent to 


¢ axis and ¢ axis 
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= ratio of specific heats 


boundary-layer thickness 


3 = curvilinear coordinates 
” similarity variable 
V l ™ vx also 
J l ve 
” cylindrical coordinate 
4 average turning at x/L 
O05 
A curvature parameter 
cl vt /U® 
v = kinematic viscosity 
p density 
x = perpendicular distance from 


axis of rotation 


vi¢ functions 


satisfying con 
tinuity 
w vector angular velocity 
2:3 components of angular ve 
locity 
0(1 = 2 angular velocity compo 


nents for cylinder 
v differential operator 
Superscripts 
() = basic flow also standard quantity 
= total value (basic + perturbation 
= differentiation with respect to 7 


= compressible flow 


Subscripts 


Coordinates indicate partial differentiation 


(I) INTRODUCTION 
6 bw IMPORTANT PROBLEM of three-dimensional 
- boundary-layer flow occurs in almost any real, 
Be 
cause of its complexity, this problem is generally 


external or internal, aerodynamic configuration. 


avoided, and either a two-dimensional boundary-layer 
flow analysis is made as an approximation, or none at 
all. 
the design is conservative and the fluid is asked to 
When, 
attempts to increase the performance by requiring the 


Such procedure is usually satisfactory as long as 


perform but moderate tasks. however, one 
fluid to flow against stronger pressure gradients, then 
the available theories completely break down, with 
rather undesirable results. This is particularly well 
illustrated by the design of turbomachinery, which is 
successful, provided the stages of these machines oper 
ate at extremely moderate pressure ratios, far below 
those obtainable in cascades and single-stage experi- 
ments. However, when high pressure ratio per stage 
machine is attempted, the real fluid effects become so 
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important that they overshadow the aims of the newly 
incorporated features and the designed performance 
cannot be realized. Consequently, the important engi- 
neering problem of increase in power output, per unit 
weight of the machinery, remains unsolved. 

The real fluid effects demonstrate themselves pri- 
marily in the so-called secondary flow and in the three- 
dimensional boundary-layer flow. The term secondary 
flow is currently used to designate flows arising from 
nonuniform total pressure or enthalpy, in regions not 
bordering on physical boundaries, where the effect of 
viscosity may be neglected. Such flows usually occur 
in internal aerodynamic configurations. On the other 
hand, by the term three-dimensional boundary-layer 
flow, we understand similar flows in the relatively thin 
regions, next to solid walls, where the viscosity plays 
an important part. Three-dimensional boundary- 
layer flows occur for both the internal, as well as ex- 
ternal, aerodynamic configurations, differing in the na- 
ture of the flow outside the boundary layer. In internal 


aerodynamic configurations, the flow outside the 
boundary layer may be secondary and thus rotational, 
while in the external aerodynamic configurations that 
flow is potential. 

A critical review of the three-dimensional boundary- 
layer literature reveals a definite division into three 
separate groups. The first group consists of the in- 
vestigations of the forms of the boundary-layer equa- 
tions and discussion of boundary-layer behavior as in- 
ferred from these forms, without actual solution of the 
equations. Here significant work has been done by 
Howarth,’ Hayes,” and Moore.* The second group 
consists of general, but approximate, treatment of 
three-dimensional boundary-layer flow by the momen- 
tum-integral method. The laminar case in this group 
has been investigated by Timman,‘ and the turbulent 
one by Mager.® The validity of these methods de- 
pends to a large extent on the choice of the velocity 
distributions in the boundary layer, and since the ex- 
perimental information is meager, attempts have been 
made to study the velocity distributions analytically. 
This leads us to the third group, which consists of 
solutions of some specialized cases of three-dimensional 
laminar boundary-layer flow. These specialized cases 
fall into two categories: (a) where the changes with 
respect to one of the independent coordinates are zero, 
or where there exists some sort of symmetry, and (b) 
where there is some characteristic parameter that is 
small. 
Bédewart,’ Howarth,’ Loos,’ and Carrier,'’ while in 


In the first category belongs the work of Sears,® 
the second one the publications of Moore''!? and 
Fogarty'* may be mentioned. 

The general treatment of three-dimensional bound- 
ary-layer flow is extremely difficult in view of the com- 
plexity of the equations. However, for many im- 
portant problems involving three-dimensional boundary 
layer, the flow outside the boundary layer may be re- 
solved into a large principal component and a small 
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cross-wise velocity. Since the pressure variation in 
the direction normal to the wall is usually negligible 
in the boundary layer, the flow there is affected by pres- 
sures existing outside the layer—that is, those govern- 
ing the external flow. Consequently, the velocities 
in the boundary layer will have the order of magnitude 
of their values in the external flow. When the cross- 
wise velocity is small as compared with the principal 
velocity, but otherwise not specified, the boundary 
layer will be only slightly three-dimensional, and the 
small perturbation procedure will be applicable to a 
great variety of external flows. The object of this 
paper is to treat the three-dimensional laminar bound- 
ary-layer flows over flat and sharply curved surfaces 
under just such simplification. 

SURFACE 


(II) BouNDARY-LAYER FLOW OVER FLAT 


Boundary-Layer Approximation 

The motion of an incompressible steady viscous 
fluid, flowing relative to a coordinate system rotating 
with a uniform angular velocity is described by the 
momentum equations, 


—qX(VXq) + oX(wXb) + 20Xq = 


i | 
-v (2+ = a°) + oI¥(¥-4) — ¥X(TXa)] (1 
p 


and the equation of continuity, 
V-q ee (2) 


Making the usual boundary-layer simplifications—that 
is, neglecting all terms of 0(6) as compared with those 
of O(1), 
ary laver equations, 


one obtains from Eqs. (1) and (2) the bound- 


u*u,* + v*u,* + w*u,* + 2ww* — w*xx; = 
* 
= Pr TT Vy” ° 
p 
») * * 2 
2 (wit =e) - xX, = = Py (la 
u*w,* + v*¥w,* + ww.* — 2wou* — w2xx, = 
* 
-_ p + VW yy 
p 
and 
u,* + v,* + w.* = 0 (2a) 


These equations differ from those used in two-dimen- 
sional theory, not only by the presence of an additional 
equation in the z direction (for notation see Fig. 1) 
but also by p, being of 0(1) instead of 0(6). This occurs 
because a pressure gradient is necessary to balance the 
effect of the centripetal and Coriolis forces. But the 
total change of pressure throughout the boundary layer 
along a normal to the wall is still of 0(6), and although 
this is larger than on a nonrotating flat plate [where it is 


of 0(6°)], it may still be neglected. A reasonable as- 
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sumption, therefore, is to consider p as a function of 


x and g only. 


External Flow 


Che external flow is defined as the flow where the 
effects of viscosity may be neglected, regardless of its 
character, which may be rotational or irrotational. 
This definition is different from that used in the con 
ventional mechanics of the boundary-layer flow, and 
it implies a division of the total flow field into a region 
where the effects of viscosity are neglected and a (still 
thinner) region where such effects are considered. If 
the external flow is rotational, the Eulerian equations 
must be used (instead of the conventional Bernoulli's 


relation) for the variation of pressure, 


wX(wXb) + 20 Xq = 


l | 
— Vv ( pt+ a’) (3 
p 2 


In the present analysis it will be assumed that the ex 


q X(¥ Xq 


ternal flow may be only slightly rotational, so that the 
derivatives (,* and JJ’,* are small as compared with 
(°.*. This simplification permits one to specify the 
pressure in terms of L*, II’,* and w alone, even if 1™ is 
W*, and 


the angular velocity w are assumed known. In _par- 


of 0(6). The external flow components l™*, 


ticular one takes: 


U* = U + U(x, 2); W* W(x, 2 
© = iw, + jo, + Rw (4 
ly < Uy 4 iW < { ai we SH @, Ww 


as describing the external flow and the angular velocity. 
The velocity ( and the angular velocity are assumed 
constant. The COMPONENT we 1S taken of the same order 


when compared to w, as l’ when compared with 


L 
Perturbed Flow Equations 


Corresponding to Eqs. (4), the flow in the boundary 


laver is expressed by 
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u” = u(x, VY) + U(X, Y, 2 
y* = (x, y) + v(x, ¥, 2 (9) 
w* = WX, V, 2 


where again the perturbation terms u, v, w are assumed 
small, when compared with the basic flow terms u° 
and v’. It is clear from the expressions (4) and (5), 
that, since lL” 
(la) and (2a) must result in the zeroth-order terms form 


is a constant, the substitution into Eqs 


ing the Blasius ‘‘flat-plate’’ equations. The first-order 
terms then give the relations, for the perturbation ve 


locities u, v, %, 


uu, + uu,” + vu, + vu,” = LU, + vu, 
uw, + vw, + 20(L" — vu) = LW, +rww,, (6 
U; + + WwW, = V0 
The solution for the basic flow is known: u°? = L°F’ 
v? = (1/2) (V »pU*/x) (nF’ — F),n = y V U"/ ox, with 


F(n) being the usual Blasius function tabulated in refer 
ence 14. Furthermore, proceeding in the manner of 
reference 3, the equation of continuity may be iden 
tically satisfied by introduction of two functions y and 
¢g, 


—(¥, + ¢; (7 


u v,; Ww Pus v 
One now assumes that the external flow velocity com 
ponents and aw. are expressible as power functions of the 


v coordinate l > A,,(z)x", W p B,,(2)x", 


m 


S* E,(s)x’, and puts the functions ¢ and y in a 


We = 


form that will later prove convenient, 

Vx " . OB 
y = | » hae, ae "iia 
a Oz 


« OE ; 
2 » "lta \ 
Os 


| px = - 
¢ Yj» [> Bee + 2>> En 'g 


Here f,,,, fn, fry Zn and g, are functions of n only 


Using 
Eqs. (8) and the expressions for the basic flow veloct 


ties, Eqs. (6) become 
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AXIS OF 


- ROTATION 


Fic. 1. Notation for flow over flat surface. 


Since U, W, and we will, in general, be different func- 
tions of x and 2, to satisfy Eqs. (9) for all values of these 
variables, it is necessary that the quantities in braces 
vanish separately. This is possible if all the coefficients 
of various powers of x vanish independently, or in other 
words, all the expressions in brackets in Eqs. (9), when 
set equal to zero, become the differential equations 
which functions fm, fn, fr, Zn, and g, must satisfy to 
make this true. Since the perturbation velocities 
(relative to the surface) are now given by 


OB,, + ~ 7 
> = >> Aux” my = > Oz * *<, i 
BE. se 
2 datas Ay 
» uC 
w= z. Bn"s.' + 2 ; > Ex" 4 
n Y 
I v (10) 


2 Ux 


1 A a {1 + 2m Vos — ie ‘n| + 


OB, , fee 
> yttl [(2n + Sin t+ 22n — Jn n| + 


n OZ 
OF, ,42 oe 7 
2 Sat? [r+ Sif. + 28, — Si'aly 


the appropriate boundary conditions are: 


fm(O) = fm'(O) = fr(O0) = fr'(0) = f,(0) = 
f,'(0) = gn(O0) = gn’(0) = g,(0) = g,'(0) = 0| bee 
m() = g,(~) = 1; f,'(~) =f,'(e) = 

g,/(~) = 0] 


As may be seen from Eqs. (9), the functions g are 
completely independent of the functions f, which means 
that, to the first order of approximation, the cross-wise 
flow w is completely independent of the principal flow 
perturbation u and is affected only by the gradients of 
the external flow W and the component of the angular 


velocity normal to the surface w,. When the cross-wise 
velocity w is not a function of the z coordinate, the 
equations for the principal flow perturbation and the 
cross-wise flow perturbation are completely uncoupled. 
This means, that according to first-order theory, the 
changes in the principal flow velocity profile, and in 
particular the tendency toward the separation, are 
completely independent of the existence of cross-flow 
and will occur as if the flow were quasi-two-dimensional. 
This result is similar to the one obtained by Sears® and 
applies even when the « velocity varies in the z direc 
tion. The cross-wise flow w brings about changes in 
the principal velocity distribution in the form of func 
tions f,’ and f,’, but only if that cross-wise flow is a 
function of the cross-wise coordinate z. 


Two types of body forces were considered originally 
in Eq. (1): the centripetal force, and the Coriolis 
force. Of these, the centripetal force was assumed 
independent of the y coordinate and consequently 
has been completely canceled out of the equations. 
The Coriolis force, which did depend on the y 
coordinate, remained in the equations, and its effect 
demonstrates itself through functions g,’ and f,’. In 
internal aerodynamics, particularily in turbomachinery, 
it is often convenient to represent the external flow 
field as caused by some imaginary force field. It 
should be clear, from the above discussion, that such 
representation will have no effect on the boundary 
layer flow except that caused by external flow veloci 
ties, provided the force field is not a function of the y 


coordinate. 


Solutions of Eqs. (9) for the various functions /f’ 
and g’ were obtained by the relaxation method and are 
shown in Figs. 2(a-d). The inspection of these figures 
and Eqs. (9) reveals that the first-order solution does 
not affect the total boundary-layer thickness, which 
remains the same as that determined by the basic flow 
solution. This, of course, is not true for the momen 
tum and displacement thicknesses, which will change, 
because of changes in the velocity profiles. We note 
also, that the maximum values of the functions f,,’, f,’, 
and g,’ are considerably lower than those of functions 


fm’ and g,’. This indicates that the effect of the 


Sn 
Coriolis acceleration and of the variation of the w ve 
locity with z coordinate is small when compared with 
the changes in velocity distribution brought about 
by the gradients of the external flow in the x direc- 


tion. 


Effect cf Compressibility 


To estimate the effect of compressibility on three 
dimensional boundary-layer flow, one considers all quan- 
tities relative to some standard quantities, the vis- 
cosity as varying linearly with temperature, no heat 
transfer at the surface, and Prandtl Number equal to 
unity. Using Howarth’s transformation, 
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The functions Y and ¢@ now written (for a nonrotating 


surface) as 
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t= % ux | F + > AX" f+ 
| 
. 
57 OB: “ 
lo = az 


x 5 
Qo = [0 P a Ga é| 


yield again the Blasius equation as the zeroth-order 


equation. The first-order equations now are 


, —_ ; 
23 Ax” oe + In —~ MF Tn + (; + m) F’'f, +m X 


») 


' : a mre 7 An Ald tt 
[ a g Mere [a 


re ie r . - 3 is 


; a 
aS F 


n 


Defining a new auxiliary function satisfying 


with boundary conditions 


i, - 0) =f, (0) = Ju. -'(*) = 0 


m « 
it is easy to show that the solution of the compressible 
equations is 

:< 


Fim = $m(a) + M? [fm(a) — F(a)] — 


wer y= it, 
2 m(zy¥ —_ 1) + ) Ton ‘ 


baal Y ee . (15) 
fn = Sfn(n) + - M?| fi (7) —fm cl 
(with m = n+ 1) 
_ . i os ‘ _ ivi 
Bn = 8nli) + —>— M*len(a) — F(a) 


The Mach Number is thus removed from the differen- 
tial equations and the solutions for the cross-flow of com- 
pressible fluid are constructed from those of incom- 
pressible fluid [note also that: f,,-; .' = fy=o'(%)]. 
From these solutions it is apparent that compressibility 
of the fluid causes not only the characteristic ‘‘stretch- 
ing’’ of the y coordinate but also accentuates the cross- 
flow. 


Examples 


In order to illustrate how these results are applied to 
actual cases of three-dimensional boundary-layer flow, 


. vre —~ ff Ald , : I 9 mayo 
>, Ba X ) Bn + — % — afF'h, + #11 4+ M?(1 — F’?)| = 0 


two examples were computed. The first is taken from 
external aerodynamics and concerns a thin cylindrical 
shell of circular cross section, flying at a small angle 
of yaw along a circular path, and simultaneously spin 
ning about its own axis of revolution (Fig. 3). The 
problem is to compute the boundary-layer flow on such 
a body. As has been shown by Howarth,! all the re 
sults derived for a flat plate apply to the circular cylin 
der as well. For small angle of yaw, the external flow 
velocities are given by 


Y=9,R; U=0: W = 2a0.R sin 0 — Qa 


The angular velocity 2, has no component in the direc 
tion normal to the surface of the cylinder, so that 





NUR=U° 

















Fic. 3. Notation for flow over cylinder. 
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Fic. 4 Flow in simulated s = shaped duct 
w Q; sin 6. Remembering that z = a@ and that for 


this case n = 0, since Il’ and a». are independent 
of x, we write immediately, from Eqs. (10), the ex- 





pressions for the “* and w* velocities in the boundary 
layer as 
* ea Xx . , 
a: 0Q,RE PaQ,R cos 0 j,.9 4 
a 
2x; cos Of, o 
a 

Ww" (2aQ,R sin 86 — Qea)g,—9' + 2xQ) sin 6 g,-0' 


The second example concerns the boundary-layer 
flow over a flat surface caused by a continuous force 


field simulating an s-shaped duct. Such flow is of in- 


terest, aside from the applications to internal ducting 
design, because it bears some resemblance to the flow 
over a casing of turbomachinery. In such machines, 
the successive rotor and stator stages serve to deflect 
the flow from the direction of the mean velocity, so that 


it oscillates about it in a periodic manner. In line with 


this analogy, the flow external to the boundary layer in 
this example is influenced by a body force in the cross- 
wise direction, which is proportional to the square of the 
principal velocity, and the local curvature of the stream 
lines. This curvature, to simulate the double bend, 
is taken as 


6s (: — cos 


27 





G 


In addition, the velocity in the principal direction is 


assumed to be a constant. After cross differentiation, 


to eliminate the pressure, the Euler's equations yield 


U°(Wr, — Uz) = (Gs)z — (Gi); 
so that (since LU) = Gy. = O) the external velocity is 
given by 
U* = U°- W* = W = 16 tan dl X 


r 


I 


v \ 


eee ia ai 


The integration constants determining the initial value 


Dr 2{- 
§ 


s 


of the body force and the inclination of the simulated 
duct have been arbitrarily set to zero. The picture of 
the stream lines (with the average deflection 3 equal to 
9°, and the width of duct s, equal to half its length / ) 
is shown in Fig. 4. 

Since the body force varies only in x and z direction, 


the boundary-layer velocities can be written, using 


Eqs. (10), as 


as Pia 


yee 


* = 16 tan dl (: — cos 2r 
/ 


[esx -2 (7) 


( 

Ss) 
From the knowledge of the velocity distributions, one 
can compute the shape of the ‘‘limiting stream lines” 
at the wall (y = 0) 


~) 


, and these are also shown in Fig. 4. 
It can be seen in this figure, that the boundary layer 


flows toward the ‘‘wall of the duct’ nearest to the center 


If the duct is to simulate two rows of 
blades, this means boundary layer flows toward the 
suction surface of the first row of the blades. 


of curvature. 


By study 
ing the velocity distributions at points A and B (in Fig 
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4 the scale of W* is enlarged out of proportion), we note 
that the shift of the boundary-layer flow implies also 
a redistribution of the momentum of the principal flow, 
since the principal flow at point A slows down and 
tends toward separation, while that at point B acceler- 
ates. In fact, since in the vicinity of point A the cross- 
wise flow also changes direction, somewhere near that 
point may occur a true separation of three-dimensional 
boundary layer—that is, all of the flow will leave the 
surface. The change of direction of the limiting stream 
lines occurs further downstream than the inflection 
point of the external flow. This means, that in the 
second row of blades the boundary layer will first con- 
tinue toward the pressure surface (where it may sepa- 
rate) and only afterwards, deflect again toward the suc- 
tion surface. 


(III) BOUNDARY-LAYER FLOW OVER SHARPLY CURVED 
SURFACE 


Perturbed Equations in Curvilinear System 


In the last example, the ‘‘sidewalls’’ of the duct were 
straight, starting from the leading edge, and therefore, 
no cross-flow occurred there. In actuality, however, 
the sides would be curved, and, since the experimental 
evidence points strongly to the corners between these 
walls and the bottom of the duct as regions of high 
total pressure loss, the flow in such corners should also 
be studied. 
corners as filleted, so that the boundary layer flows over 
a sharply (but continuously) curved surface. 


To simplify matters, one considers the 


For the purpose of this analysis it is convenient to use 
an orthogonal curvilinear coordinate system shown in 
Fig. 5. The transformation from the Cartesian system 
x, y, z to to the curvilinear system &, \, ¢ is given by 


$ 
x= y=d.+Acoss; Y= - f sin 6B d¢; 


0 
t 
z= cos8d¢+Asing (16) 
0 
Restricting the analysis to curvature of the ¢ axis 
c = dB/d¢, being a function of ¢ coordinate only, the 


length stretching factors used in the transformation of 
the equations of motion become 


k= 1; b=1; k= 1+ea (17) 


It should be pointed out that in this coordinate system 
the ¢ axis plays the role of a reference line and will be 
used to describe the shape of the wall over which the 
fluid is flowing. Coordinate \ will then measure the 
distance perpendicular to this wall, while the & axis is 
pointed in the direction of the principal flow, so that § = 
A = O represent the leading edge of the configuration. 
Using the standard transformation formulas (veloci- 
ties u*, v*, w* are taken here in the é, XA, ¢ directions, 


, Ww 


respectively) and applying the boundary-layer approxi- 
mations—that is, considering v* and Xd of 0(6), c of 
0(1/6) or larger, and neglecting all terms of 0(6)- 
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one obtains the boundary-layer momentum equations, 


l 
u*u.* + v*¥u,* + wn. = 
1+ cj 
* : x 
- Pe + v( u)* 4 Uy 
p 1+ cv 
[c/(1 + cad)] (w*)? = —(1/p)py 
(18 
u*w; : v*w,* + ww * - 
1 + cA 
( ae l l 
vw = — pe tux 
1+ cA pl+a°* 
€ c : 
Wy ,* + w,* — ( w* 
1 + cA 1+ cA 
and the equation of continuity, 
* * c * 
Ue + v,* + w,* + vy =0 (19 


i+a 


‘ wy 
1+ cA 


One now considers the flow as only slightly perturbed 


from the £ direction, so that the zeroth-order terms 
yield again the basic flow equations, 


Wy? + VU? = v iUy,® + [c/(1 + cA) Jury} | 


u,® + 2° [c/(1 + cA) ]v® = 0 


(20 


/ 


and the first-order terms give the perturbation equa- 


tions, 
l 
U"u, + UU,” + UU, + VU," + ee WU? = 
l ri 
— - PT (ms ita us) 
Oe, Oe c Oe 
UW, + VW) + ita vw = | 
Lo} = r (21 
= ol tA Pp tv ES + i+ 2 
C . | 
cay 
Pp, = 0(w?) 
' 
u, + Vy, + 1+ - W, + lta v=0 


Solution of Basic Flow Equations 


The basic flow equations may be solved by defining a 

stream function °(&, A, ¢), and an additional parameter 
y l/e e 

x = c(vt/U°)” related to the local curvature of the 


wall. The stream function is first expanded, 


] (22) 


y= Vv ‘Uv [F(n) + xd(n) + x2di(n) +... 


and introduced in this form into the equations of basic 
flow. 


the coefficients of the various powers of x are set sep- 


Again, since the d’s must be independent of ¢, 


arately to zero. The maximum permissible value of 
x will be small [of 0(1/n;)], otherwise the curvilinear 
system would lead to double values, and, consequently, 
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FIGURE 6.- VELOCITY PERTURBATION DUE TO 
WALL CURVATURE 
a) IN THE PRINCIPAL DIRECTION 


it is sufficient to consider only the first two terms in 
this expansion. One thus obtains the Blasius equation 
for F(m) and an equation for d as follows: 


2d’** + Qdd"’ — (d')? + (nF — 2)F" + 
F'(F'n — F) = 0 (23) 
Since the basic flow velocities, neglecting higher 
powers of x, are: 
u® = [1/(1 +xn)]U® (F’ + xd’) 
) : 
+ Kn 


v 


((F — F'n) + » (24) 


ys — - 


dole 


0 
g 


x(2d — d’n)] } 


d must satisfy d(0) = d'(0) = Oandd’'(~) =n. A 
numerical solution for d’ (reference 15) shows that 
d’ = n everywhere, and therefore a fair approximation 


for the u° velocity is 
u® = U°[(F’ + wn)/(1 + xy)] (25) 


This approximation has been used subsequently. 
Since, for small values of n, F’ may be approximated by 
F’ = 0.332n, expression (25) shows that, for small 
values of n, u° = (0.332 + x)n + O(m?). In other 
words, the slope of the velocity distribution near the 
surface depends to a great extent on the local curva- 
ture of this surface. On concave surfaces, « will be 
negative, and thus the velocity profile will be nearer 
separation (will have a smaller slope), and on convex 
surfaces the basic flow velocity profile improves. 


Solution of the Perturbed Equations 
To solve the first-order equations, one proceeds in 


exactly the same manner” as for basic flow. To sim- 
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FIGURE 6.-CONT'D. 
b) IN THE CROSSWISE DIRECTION 


plify matters but still retain a significant model, it is 
sufficient to consider the external flow to be of special 
type: U = 0, W = Dé". The appropriate expres- 
sions for ¥ and ¢ may now be written immediately, as 
a result of the flat surface investigation, 


A P Ox W ane 
y(&, $7 ™ & ar [pe V vél xX | 
(1,(n) + xh p(n) +... | | 
> (26) 
W 


0&5) = 7 V vEU® [gn(n) + 


kt, (n) + «tp n(n) t... | 
It is clear that the equation for g, must be the same 


as for the flat surface, and the equations for /, and f, 
become 


F 3 
1!" + 1,’’ - ( + ”) F'l,’ + 
9 9 


(2+ n)F'l, + Ft, — (d' — F'n)g,' = 0 | 
see F sv l nL , sre 
~ += . -t— +a F't,’ — 12n.° + 
F g } (27) 
(« a te i)e n= | mF + d’) + 


l 
=o F) | a. + 2nn + 
a ws ” , 
lim (F'n — F) = 0 


—> 


Since the perturbation velocity components are 
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d(log ec)  cé l ; 
v= Woe L, 
de 1+ AWVRe 
(28) 
| 
w= W B + c& —= h | 
V Re 
the appropriate boundary conditions are: 
/,(0) = 1,'(0) = 1,'(~) = t,(0) = ¢,'(0) = 
t, (2) = 0 
Equations (27) were solved numerically for n = 0 and 


n = 1. These solutions are shown in Figs. 6a and 6b. 
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Example 

To illustrate the manner in which the derived results 
are applied to an actual boundary-layer flow over a 
sharply curved wall, such flow has been computed for 


U = 0, W = Dé", with nm = Oandn = 1. The wall 
was assumed to have a hyperbolic shape. Assuming 
ns = 8.0, the maximum permissible value of « has 


been set at 0.125, so that at any given Reynolds Num- 
ber Ll” ve,,,,, there exists a maximum value of cé: 


Cond tax. = TOs 0") weq.. 


for which these considerations hold. In this example 
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FIGURE 7.- CONT'D 
c) EXTERNAL FLOW: U"=U°; W*=D¢ 
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FIGURE 7.- CONT'D 
b) EXTERNAL FLOW: U*=U°; W*=D 
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THREE-DIMENSIONAL 


U°/vemez, = 634, so that £,,, = 10/Cmax. The results 
of the computations are shown in Figs. 7 (a—d), where 
are plotted the velocity ratios, and the deflection from 
the external flow direction. It should be noted that 
at the computed position & = &,,,,, the cross-wise com- 
ponent of the external flow WW* is the same in both 
cases, so that the differences can be traced directly to 
the cross-wise pressure gradient. Physically, such an 
external flow may be thought of as a vortex of constant 
strength for » = O, and a vortex of an increasing 
strength for m = 1.0. Examination of Figs. 7a and 7b 
shows that, for m = 0, an accumulation of low momen- 
tum air occurs on the “‘upstream”’ wall of the corner. 
When x = 1.0, however, Figs. 7c and 7d reveal that now 
the large accumulation of the low momentum air occurs 
wall, with corresponding large 
We thus see that, when there 


on the “downstream” 
deflections in that region. 
is no pressure gradient in the external flow (n = 0), 
the wall curvature causes a weak secondary flow as 
shown in Fig. 7b. When, on the other hand, the ex- 
ternal flow has a cross-wise pressure gradient (” = 1), 
the secondary flow due to that pressure gradient ap- 
parently (Fig. 7d) is in such a direction as to hinder 
the principal flow on the “downstream” wall, so that 
almost all of the flow ‘‘climbs’’ the wall in that region. 


(IV) CONCLUDING REMARKS 


The analysis of three-dimensional, steady laminar 
boundary layer with small cross-flow has indicated that 
many important problems in this class can be solved 
with relatively little effort, thus affording a greater 
understanding of real fluid flow. Although the accu- 
rate determination of the separation point is certainly 
beyond the scope of small perturbation analysis (and 
also beyond the standard boundary-layer theory), the 
factors that affect the separation of three-dimensional 
boundary layers are clearly pointed out. Among these 
factors, the cross-wise variation of flow is indicated 
as important in affecting the separation. In_par- 
ticular, the examples show the strong tendency toward 
the separation as a result of cross-wise variation of 


LAMINAR 


BOUNDARY LAYER 845 


cross-flow, both on a cylindrical shell and in a simulated 
s-shaped duct. Similarily, the sharp, laterally varying, 
concave curvature of the surface tends to separate 


even the basic flow profile. 
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Note on Photoelastic Study of a Swept-Wing 
Model 


George Gerard and Ralph Papirno 
Research Division, College of Engineering, New York University 
July 10, 1954 


Sew PURPOSE of this note is to report briefly on preliminary 
photoelastic results obtained on a cemented swept-wing 
model using a ‘‘back-silvering’’ technique 

The extension of photoelastic techniques to relatively complex 
built-up structures has been hampered until recently by the lack 
of suitable photoelastic materials which can be cemented without 
introducing stresses at the bond and are not subject to slip 
Becker! used cemented models in an airframe bulkhead investi 
gation but found that although the technique gave a satisfactory 
first approximation, creep of the cement under load was a serious 
problem 
which is suitable for photoelastic work and which can be satis- 
factorily bonded. 

We have been intrigued by the possibility of investigating com- 
plex built-up structural components by use of cemented model 
With the availability of Paraplex, an exploratory 
investigation has been conducted on an idealized swept-wing 
model. 


Recently, Williams? referred to the resin Paraplex 


techniques. 


Fig. 1 shows the cemented swept-wing model and Fig. 2 
shows the model installed in a loading jig. 
If this model were loaded and viewed in the usua: transmission 


polariscope, no fringe pattern would be observed. Since the 


upper cover is in tension and the lower is in compression the 
respective fringes annul each other. To surmount this diffi- 
culty, the inside surfaces of the model have been silvered and a 
doubling polariscope arrangement has been used to determine 
the fringe patterns. By use of the silv. cing technique, it is thus 
possible to investigate selectively the areas of interest in the 
model in addition to obtaining the increased sensitivity of a 
3, the fringe pattern is shown for 


doubling polariscope. In Fig 


the complete tension cover of the model. Fig. 4 shows the pro 


gressive increase in fringe order with load for one half of the ten- 


sion cover 
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Fic. 4. Swept-wing model fringe photographs, light field retiec 
tion polariscope. 


Further development of the techniques noted here can possibly 
lead to the widespread use of photoelastic analysis for complex 
built-up structural components 
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Rectangular Plates with Stringers and Ribs 


J. Irving and N. Mullineux 

Senior Lecturer in Applied Methematics and Lecturer in 
Mathematics, Respectively, Department of Mathematics, The 
University, Southampton, England 


July 11, 1954 


NOTATION 


1,6 dimensions of plate 

h thickness of plate 

x, ¥ = coordinates 

q(x, 3 intensity of lateral loading 

w(x, ¥ lateral deflection of plate 

w(m, n sine transform of lateral deflection 
D flexural rigidity of plate 


xj position of rib 





position of stringer 
\ Vy normal compressive forces per unit length of plate edge alon, 
the x- and y-axes, respectively 
i area of stringer 
B flexural stiffness of stringer 
1 area of rib 
B’ flexural stiffness of ril 
FP effective force at point due to stringer 
(@) effective force at point luetor 
x x b(y Dirac delta function 


Poisson's ratio 


(1) INTRODUCTION 


_— USE OF THE DIRAC DELTA function in conjunction with the 
finite sine transform simplifies the treatment of problems 
involving the loading of rectangular plates, for which at least 
Recently 


one pair of opposite edges are simply supported.' 


Radok? has discussed the problem of a simply supported rec 


tangular plate strengthened with ribs and stringers by the 


usual tentative methods Moreover, he considers a stringer 


or rib as being equivalent to the limit of a certain triangular 


loading along a strip of finite width It is, however, more 


direct to represent the forces due to the stringer (or rib by 
a Dirac delta function and to carry out the subsequent analysis 
by means of sine transforms as described in the following para 
graphs. 

> 


(2) Isorropic PLATE WITH ONE STRINGER OR RIB 


rhe deflection w(x, y) of a simply supported rectangular plate, 
subject to uniform compressions N, and N, in the x- and y 
directions, respectively, which are taken along adjacent edges of 
lengths a and 3, satisfies the equation 
' rw 07% 
Viw + (Nz/D - + (N,/D D 4 
Ox? OV 
P(x) (9 y)/D (2.1 


where g(x, ¥) is the intensity of the applied lateral load and P;(x 
is the effective load at the point (x, y,) due to the stringer at 
y= ¥; 

Using the double finite sine transform! (or Sneddon‘), and 
taking account of the boundary conditions for a simply sup 


ported plate, Eq. (2.1) becomes 


L™ Yar + ef ~~ aD ~ wD 
b 
. marx . nay 
q(x, ¥) sin sin —— dxdy + 
J0/70 a D 
_ nry : _ mr i 
sin . P(x) sin dx (2.2 
b Jo , 
where 
2a (*b 
_ Marx . nry 
w(m, n) w(x, ¥) sin sin dx dy 
7 a b 
J0/0 1 


The well-known property of the delta-function, 
_ nny ny 
sin 5(y ydy = sin 

b ; b 

0 

99 


has been assumed in deriving Eq. (2.2) 
By the usual inversion for 


Eq. (2.2) determines 
w(m, n) for any lateral load q(x, y) 
mula, the deflection is given directly by 


} mu Mi qwx uwy 
“w(x, y) = > S w( mi, m) sin sin (2.3 
- } —s — } 
@O m=1 n=1 a, , 
For g(x, y) = 0, Ny = 0, w(m, n) from Eq. (2.2) is equivalent 


to Cmy' as given by Radok? [Eq. (2.3.1 


The case of a rib is treated in a similar way 
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(3) DISCRETE STRINGERS 
Consider now the case of stringers at y; where 7 = 1, 2, ot 
For comparison with Radok’s results,? we take g(x, y) and NV, 


to be zero. 


The generalization of Eq. (2.2) for f stringers involves the un 
known quantities 
fra 
: marx 
P(x) sin dx (g=1,2 
/J0 a 


which are required to satisfy the internal boundary conditions 


aty=y These are simply 
Ow A,N, Ow 
B, 1 = = —DP;(x) lyxy 3.1 
ox? h Ox? 


& = I, 2, e 


j This form of the boundary condition avoids 
the analytical difficulties of Radok.? 


From Eqs. (2.2) and (2.3 


generalized for f stringers 


‘ ; 
w(x, y) = >: 2: y x 
- 1 — — oe 
ao » l 1 i ] 
ra ‘ 
_ nny Max 
sin P(x’) sin dx 
b Jo u _ max . nry = 
- sin sin 3.2 
fm? n2)2 mex? N, a 0 
ta2 bz { a? D 
re * 7 ? , cs : WAX 
Substituting Eq. (3.2) in Eq. (3.1), multiplying by sin and 


integrating from 0 toa 


1 1 
_ Anyi. ny 
sin sin ea 
D / M1 aX 
: P(%’) sin dx = 
[ jm? n*| m?r?N Jo a 
Tr’ . + 
| la? D ( a®?D 
ea , 
; NAN 
D P(x) sin ax’. e= Ii? 
J a 


This is a system of f homogeneous linear equations in the 
a 


. . max’ 
variables P(x’) sin dxv',i = 1, 2, a 
0 a 


The elimi 


nant is identical with Eq. (3.1.3) of Radok 

Eq. (3.2.1) of Radok’s paper is readily obtained by the abov« 
methods in the case of discrete stringers and ribs 

The case of distributed stringers which is equivalent to an 
orthotropic plate may be treated in the same way whether or 


not ribs are present. The internal boundary condition for a 


rib atx = x; is 
, Ow A’ N, Ow 
B’) ee : _= DQi(y) |x x, (1 = 1, 2, , 2) 
oy! h ody? 
and the right-hand side of Eq. (2.1) contains the term 


D for each rib 


The above calculations illustrate the concise and constructive 


Oi y)i(x — x 


nature of the transform approach in conjunction with the 6-func 
tion technique 
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On Motion of Submerged Cylinder 


H.S. Tan 

Department of Engineering Mechanics, University of Notre Dame 
Notre Dame, Ind 

July 17, 1954 


—_— TWO-DIMENSIONAL problem of a submerged cylinder moy 

ing with constant forward speed c, under the assumption of 
irrotationality of fluid motion, and in coordinate system moving 
with the cylinder, reduces to the determination of a velocity po- 


tential ¢(x, y), satisfying the differential equation: 


with boundary conditions: 


(1) on free surface, y = 0) cp r go, = UV ~ 


Il) on surface of cylinder gn = 0 3 
and the physical limitation that the upstream disturbance must 
vanish at infinity —i.e 
Lim [¢,? + ¢,7] = 0 j 
_— 
The solution for a deeply submerged circular cylinder was first 
obtained by Lamb! in 1913, in the form: 


QoQ = cx( 1 + a*/r* T X\X, ¥ ” 


the depth of submer 
The effect of x on the 


Lamb's solution is accordingly 


where a denotes the radius of cylinder, 
sion; a< h, and r = [x? + (y + h)? 
surface of cylinder was neglected 
a first approximation in the sense that the cylinder was replaced 
by a single doublet Later, in 1936, Havelock? solved the sam« 
problem exactly. He expressed the solution in the form of an in 
finite series, the coefficients of which can be so determined that 
both boundary conditions are taken care of. Havelock’s solution 
amounts to replacing the submerged cylinder by a system of dis 
tributed sources and sinks of zero total strength inside the cyl 
inder, and replacing the free surface by a properly adjusted imagt 
system 

In 1935 Keldish, 
tial 


by introducing the complex velocity poten 


w=-ao-+iy, =Nx+1Y 6 


and rewriting the boundary condition (1) in the form 


aw du 
a CU 


was able to utilize Schwarz’s reflection principle for constructing 


solutions involving free surface. Indeed, he obtained in this 


way the solution for a submerged moving eddy. Generalizing 
Keldish’s result, and making extensive use of Cauchy's integral 
theorem, Kotchin,! in 1987 was then able to establish an exact 
solution for submerged cylinder of arbitrary shape in the form 
of an integral equation. A close examination of Kotchin’s solu 
tion reveals that here the cylinder was replaced by a distribution 
of sources, sinks, and vortices. Thus in addition to drag force, 
the lifting force due to circulation can also be calculated 

An entirely different approach to the solution of the problem of 
submerged circular cylinder was given in 1938 by Sretensky 
Here the disturbance stream function yY was considered. On 
taking account of the facts that: (1) upstream disturbance 
vanishes at infinity, (2) free surface is given by vanishing stream 
function y = 0, (3) cylinder surface is given by constant stream- 
the two 


line ¥ = a, and (4) the relation ¥ = —cy + wy holds; 


boundary conditions can be put in integrated form as follows: 
(I’) free surface, y = 0 vy = (g/oy¥ (8) 
(II’) surface of cylinder y=a+tcy (9) 


Sretensky further succeeded in transforming the entire flow 


field into a region enclosed between two concentric circles in ¢ 
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ie by a transformation: 


>—wX 
i= : A = h? — q? ’ «6 (10) 

+ XN 
where ¢ = 1 corresponds to the free surface, and ¢) = [(A 4+ 
\ = p the surface of cylinder. The corresponding trans 


formed boundary conditions expressed in polar coordinates 


’ of ¢ are 
vr = py /(1 COS w)|r=1, (uw = —Ag/c- (11 
y¥ oa @ a T cX(1 - p~ (1 2p COS w™ 7 p* (12 


hese, upon further introducing the complex potential w, take 
on following form: 


(1”) T\2uw + we (1 — &)? cj} =U, co; = 1 (13 
II”) Iw _ cs] = a, ¢ =p (14 


Although w here is a holomorphic function in the ring p 


¢ <1, it is no longer single valued. Indeed, a solution was 


obtained by assuming following series form: 


w=gqgin¢g + =. Pm =” (15 


where p,, are determined by substituting (15) into conditions 
P Z ; oy sin 
I”), (II1”), and then equating to zero the coefficients of ma 

c 


os 


It is to be noted that Sretensky’s solution as presented in refer- 
ence 5 is not correct; this is because his boundary condition cor- 
responding to (I”) was incorrectly transformed. Indeed, it will 
be valuable if a complete solution can be carried out with the 


correct boundary conditions (I”) and (II”) arrived at above 
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Note on the Propagation of Infinitesimal 
Disturbances in Gases According to the 
Navier-Stokes Equations 


M. Lessen 

Professor of Applied Mechanics, University of Pennsylvania 
Philadelphia, Pa. 

July 14, 1954 


Tp TRUESDELL’s! comprehensive memoir on the propagation 

of plane sound waves according to the Navier-Stokes equa- 
tions, the dispersion and damping of longitudinal waves is re 
viewed and studied. The analysis, however, may be generalized 
to any stress disturbance and the results presented in convenient 


form. 
NOMENCLATURE 
Xi space coordinate vector 
i time coordinate 
ui velocity vector 
p = density 


r pressure 
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Ti stress tensor 
temperature 
K thermal conductivity 





Mi; Mw viscosity coefficients (shear, bulk 
specific heat at constant volume 


specific heat at constant pressure 


perfect adiabatic velocity of sound 
q ig 
q Og/ Ox 
Repeated indices connote summaticn 


ANALYSIS 


The Navier-Stokes equations for the motion of a viscous fluid 


may be written as follows: 


Continuity + Uke = 0 


Momentum pi; = ry, 


Stress—Rate of Strain 


r 4 pb, = (2/3)(p ay) 6, + pay ra 
Energy pe | = (AK7 TT ‘ 
State 
f 
— = R7 
p 
K = K(7 
mw = wil 7 
be = wo(7 


The equations may be linearized by considering each of the 


dependent variables to be of the form 


where g'< @ and @ is a particular (constant in x, and /) solution 
of the equations. Since it is assumed that the gas is grossly at 
rest, 

u i) 


The linearized equations for the disturbance therefore ar¢ 


; 1 Op’ 
Continuity + trp = O 
p Ol 
’ 
On j 
Vomentum p Fy 
dt 
Stress— Rate of Strain 
tri + p'bei = (7/3) (hs fi)ui’. 16 Tr Mil iy TM 
Energy oT’ in 
p ¢ = AT, pu 
ol 
ty [ / 
State — - : 
p p / 


If the stress—rate of strain equation is substituted into the 
momentum equation, the resulting equation differentiated by 
x, and summed over the repeated ‘7’ index, the resulting velocity 
divergence terms may then be eliminated from all the equations 
using the continuity equation. The temperature and pressure 
disturbance terms may be eliminated in combination with the 
energy and state equations. After rearranging, the resulting 


equation for the generalized density disturbance is 


re) 07 p’ a a K (<= a? 
AG 7 a9’ a) pte \ ate oy” ‘) 
$ji; + 2g. d (2 K 
3p or \ ot pe 


where a? yRT. 
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The equation for the density disturbance may be put into 
dimensionless form by the following transformation and group- 


ings. 
a 
Y= “" 
WwW 
l 
{= -?° 
Ww 
3p a? " - 
Re = “Reynolds Number 
(44, T 2») Ww 
(4; + 2f2)cp 7 ” 
= — Z “Prandtl Number 
3K 
where w = characteristic disturbance circular frequency. 


The resulting dimensionless equation for the density disturb 
ance (with all primes and asterisks dropped) is 


ro) (2° ) Y (2 1 ) 
P. kk = 2” 6 sé eke = 
ot \ of? aRe \ Oot? y oat 


1 Oo (2 ¥ ) 
sie P, kk 
Re ot \ ot aRe ; u 
or in vector form 


0 oO” j oO? 1 
( _v-V p= tv (% - vv) e4 
ot \ Or? oRe ot? Y 
1 Oo ° / 
vv ( en vr), 
Re ot ot aRe 
The form of the final, dimensionless equation is of interest. 
For o ~ Oand Re > ~, the density disturbance is given by the 
three-dimensional wave equation where the velocity of propaga- 
tion of a longitudinal wave is the usually considered ‘‘perfect”’ 


adiabatic velocity of sound 


oO? 
= T29 h pp =O 
ot? 


The foregoing corresponds to the first correct calculation of 
the inviscid velocity of sound by Laplace? who assumed that the 
gas went through an isentropic state change. 

For the case of o ~ 0 and Re ~ o@ such that oRe — 0, the 
density disturbance may be given by the three-dimensional wave 
equation wherein the velocity of propagation of a longitudinal 
wave is the frictionless, isothermal velocity of sound. 


oO? ] 
- V-V p=0 
ot? ¥ 


This case corresponds to the first calculation of the velocity of 
sound by Newton who incorrectly assumed that the gas went 
through an isothermal state change. 

All intermediate cases are given by the full equations and 
correspond to the investigation of Kirchhoff‘ for the plane, longi- 
tudinal wave. Detailed numerical results for the case of a plane 
longitudinal wave are given by Truesdell.! 

It is interesting to note, from the linearized disturbance equa- 
tions, that a pure shear (transverse) mode of disturbance does 
not interact with a longitudinal disturbance. This dilatation- 
free disturbance therefore behaves as if in an incompressible, 
viscous fluid, for which this problem has already been completely 


solved. 
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The Effect of Heat Transfer on the Behavior of 
a Compressible Laminar Boundary Layer 
Under Favorable and Adverse Pressure 
Gradients 


Ting-Yi Li and H. T. Nagamatsu 
California Institute of Technology, Pasadena, Calif. 
July 17, 1954 


— EMPLOYED the differential analyzer to evaluate 
solutions of the Falkner-Skan? equation 


K'" = — KK" + 6[K"” — 1] (1 


with boundary conditions K = K’ = 0Oatyn = 0, K’~1lasy— 
©. It was shown that to make the solutions for 8B < 0 unique 
and appropriate to the boundary-layer applications, the condi 


tion at infinity was replaced by 
K’ — 1 from below, and as rapidly as possible, as7 ~ ~ (2 


The value of K”(0), required to satisfy the above condition for 8 
negative, was found to decrease for increasing (8 and to become 
zero at a value 8 = By) = — 0.199. The value of %) is significant 
because it determines the point at which an incompressible lami 
nar boundary layer breaks away from the surface. The aim of 
the present note is to extend Hartree’s computations into the 
compressible flow régime and to present some general results re 
garding the effects of heat transfer on the behavior of a laminar 
boundary layer with free-stream pressure gradient. 

In a recent paper,*® * we have shown that if the free-stream 
velocity distribution satisfies a certain differential equation and 
if the surface temperature assumes a constant value, the com 
pressible laminar boundary-layer equations can be reduced to 


K'" = —KK" + B[K” — G] (3 
G" + KG’ =0 (4) 


with boundary conditions 
n = 0, K = K'=0 G =.= (5 


n> « K'=> 1 G—>!1 (6 


Here AK’ = u/m, G = H/M, and 7 is the similarity variable (cf 
reference 3). H and u denote, respectively, the specific total 
energy and the velocity component parallel to the surface. 7 
denotes the absolute temperature. Subscript 1 refers to the free 
stream; subscript 0, the stagnation condition in the free stream; 
subscript w, the condition on the surface. In the deviation of 
Eqs. (3) and (4), the simplifying assumptions that (1) Prandtl 
Number is unity, (2) the viscosity coefficient is a linear function 
of the absolute temperature, and (3) the specific heats are con- 
stants, are used. 

If the surface is insulated, we must have G = 1. Such a case 
corresponds exactly to the case analyzed by Hartree. The 
present formulation permits us to discuss the effect of heat trans 
fer on the laminar separation in the sense that when G varies 
from 1—i.e., the surface is noninsulated—we can find the cor- 
responding variation in 8. If we specify that A”’(0) must 
change continuously with G, we can derive the following rela 
tion® 


AB = —f) AG (7) 


In the deviation of the above equation we have assumed that (1) 
for G = 1, K"(0) vanishes when B = By) = —0.199, (2) for G = 
1 + AG, K''(0) would vanish when B = 6 + A8, where AG< 1, 
and A8 << Bo = O(1). In other words, Eq. (7) strictly applies 
only to the case of a surface with a very small amount of heat 
transfer. It is seen from Eq. (7) that the boundary layer on a 
heated surface will separate under a comparatively weaker ad- 
verse pressure gradient than in the case of an insulated surface; 
and that the boundary layer on a cooled surface can sustain a 
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stronger adverse pressure gradient. Recent numerical analysis 
of Eqs. (3) and (4) on a REAC* confirms these deductions. 

In Fig. 1, the values of A”’(0) of the numerical solutions ob- 
tained on the REAC are plotted against 8 for various 7,’s. Since 
the local skin friction can be shown to be directly proportional to 
K"(0), the following observations can be made with regard to the 
boundary layer under a favorable pressure gradient in the free 
stream: 

(1) For a fixed 8 > 0 the value of A”(0) decreases as Ty, de- 
creases; that is, in the presence of a favorable free-stream pressure 
gradient, cooling the surface tends to reduce the skin friction. 

(2) For a fixed T,,, the value of A”(0) increases with 8 and the 
rate of increase is higher at an elevated surface temperature 
Therefore, on a surface with fixed temperature, a larger favor 
able pressure gradient tends to increase the skin friction. 

Besides these interesting observations, we note that an increase 
in the adverse pressure gradient in the free stream tends to de- 
crease the skin friction with A”(0) eventually becoming zero as 
8 = B. We have thus obtained four values of o's corresponding 


to the following cases: 


ao 
(1) insulated surface, = = 1, By = —0.196 
Vw ; 
(2) cooled surface, r = 0.6, Bo = —0.240 
le 
—- = 0) Bo = —0.313 
7 


* Sincere thanks are due to Dr. C. Gazley, Jr. of the RAND Corporation 


for his kind help in obtaining these results 


FORUM 851 


3) heated surface, r, =2 Bp = —0.124 


Therefore, we conclude the effect of surface cooling is to enable 
the boundary layer to sustain a stronger adverse pressure gra 
dient. The effect of surface heating is just the opposite 

It is perhaps interesting to note that A”(O) = 0.47 for 6 = 0 
with various 7),’s.¢ This is significant because when 8 = 0, we 


have by Reynolds analogy 
2C¢, G'(0 
Cr, (Tw/To) — 1])K"(0 


where C;, and C;, represent, respectively, the local coefficients of 
heat transfer and skin friction based on the free-stream quantities 
Therefore, we anticipate that the value of G’(0)/[(7y/Ty) — 1 
for 8 = O should not be affected by changes in 7,. This fact is 
actually verified from our data obtained on the REAC, and can 
easily be proved analytically 

Illingworth‘ recently argued that when the surface is heated 
the fluid near it is lighter and is therefore more readily accelerated 
or retarded by a pressure gradient. In other words, the effect of 
pressure gradient is enhanced when the surface is heated and our 


data seem to substantiate this argument 
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+ For such a case, Eq. (3) becomes the Blasius equation 


Second-Order Methods in Inviscid Supersonic 
Theory 


J. Fell and D. C. M. Leslie 

Armaments Division, Sir W. G. Armstrong Whitworth Aircraft, 
Limited, Baginton, near Coventry, England 

July 25, 1954 


Se FULL EQUATION for the inviscid supersonic motion of a 
polytropic gas is very complex, and can only be solved for 
a few special cases. For general problems, the usual procedure 
is to assume the existence of a velocity potential 

Q = U(x + ¢) (1 


where U is the free stream velocity. If it is assumed that @ 
is so small that terms of order | ¢{? may be neglected, one arrives 
at the usual linearized equation of motion. In the second-order 
procedure associated with the name of Van Dyke,! one sets 

2 = Ux+ort ¢o) (2 
assumes that @ satisfies the linearized equation, and neglects 
terms of order {| ¢}?. The resulting equation for ¢ is given by 
Van Dyke: it contains both double and triple products of ¢ 
The triple products can be ignored in some problems, and treated 
approximately in others, so that we need only integrate precisely 
the equation involving the double products, 


[ }*¢ = Puy ez ~~ Berger: = 2.M!?{( N - 1) 8*6.@22 + 


PyPry Tt P2Pz (3) 


(the notation is the same as Van Dyke's). 
The first step is to find a particular integral (p.i.) of this equa- 
tion; the complementary function is then chosen to satisfy the 
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boundary conditions of the problem. Van Dyke has been able 
to do this for two-dimensional and axisymmetric flow; an ex- 
tension of this work to the general case is obviously most desir- 
able. 

Now in all practical problems, the sources of @ will lie on some 


surface 2, so that 


1 | ; f dxdy, 
>= - (Xa, ¥ 
" w settee R-R, 


where we are using the ‘‘hyperbolic vector notation” 
A-B = A.B, — B%( A,B, + A-B;) 
The first step in the problem is to find a p.i. for 
] 


4 
R—R||R — R,| sind 


Cy = 
Using two different systematic procedures, the authors have 
found four related p.i.’s of Eq. (4) 


h=t . log ||R: — Ro| + |R—R,| + |R — R,|| 
Any combination of signs may be taken in the argument of the 
logarithm: it will then be obvious which sign should be taken 
for the whole expression. The most useful combination of these 
p.i.’s is 

1 ' (R—R)-(R—-—R) —|R—R,||R—-R, 

» = 49: SR _—R)(R-R)+|R-R||R—-R|| 
It will be noted that A, vanishes on the Mach cones |R — R,; = 
0; |R — R.| = 0. Using this expression, one may obtain the 
known p.i.’s for axisymmetric, two-dimensional, and yawed 
conical flows. 

Eq. (5) has been applied to the calculation of aerodynamic 
derivatives of wings of finite aspect ratio and thickness and the 
lift of bodies of revolution. Results and a full description of this 
work will be published later. This work forms part of the re- 
search program of the Armaments Division of Sir W. G. Arm- 
strong Whitworth Aircraft, Limited, and is published with the 
permission of the firm. 
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Buckling Stress of Thin Cylindrical Clamped 
Shells Subject to Hydrostatic Pressure 


P. P. Bijlaard 
Consultant, Bell Aircraft Corporation, and Professor, Cornell 
University, Ithaca, N.Y 


July 26,1954 


I A RECENT NOTE! Professor W. A. Nash refers to the existing, 
rather involved methods for calculating the buckling stresses 
of cylindrical shells with various boundary conditions subjected 
to hydrostatic pressure. He then gives a rather simple deriva- 
tion of an approximate formula for the buckling stress of a thin 
cylindrical shell that is clamped at the edges. 

In this connection it may be of interest to present two other 
simple methods which lead to explicit formulas for the critical 
hydrostatic pressure of clamped cylindrical shells. These 
methods have also been worked out by the writer for the case of 
elastic end restraints. 

In reference 2 the writer has considered the deflecting and re- 


straining forces that act upon a longitudinal strip of a shell at the 
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Fic. 1. Clamped Cylindrical Shell. 
Fic. 2. Simply Supported Cylindrical Shell 











moment of incipient buckling. The deflecting forces are caused 
by the compressive membrane stresses in the shell due to the 
changes of curvature. The restraining forces are due to the 
excess membrane forces, transverse shear forces and twisting 
moments along the longitudinal edges of the strip. He proves 
that for shells where the ratio of the half wave length of buckles 
in the axial direction and that in the circumferential direction is 
sufficiently large, the resulting deflecting forces (deflecting minus 
restraining forces) acting on such a strip at any point are propor- 
tional to its deflection w. Hence, these resulting deflecting forces 
can be represented as Cw per unit length, where C is a constant 
for a given shell, dependent on the external pressure, but inde- 
pendent of the boundary conditions at x = +c (Fig. la). These 
resulting deflecting forces, acting on an element of length dx of 
the strip, have to be in equilibrium with the restraining forces 

dQ, where Q, is the transverse shear force in the strip. Since 


QO. = — Nd*w/dx*, where N is the flexural rigidity of the shell, the 
equilibrium equation —dQ, = Cw dx leads to the differential 
equation of the strip 
Nd*w/dx* = Cw (1) 
This equation may also be obtained from Donnell’s equation® 
for the present case, where o, = '/20, (Fig. 1), 
; Et 0'w Ow 1 Ow 
NViw + + to,V+4 + = 0 (2) 
a? Ox‘ oy? 2 Ox? 


by taking into account that with the above assumption the de- 
rivatives in the axial direction can be neglected with respect to 
those in the circumferential direction, so that Eq. (2) reduces to 


_ Ow Et o'w Ofw y 
N + — — +t, — =0 (3) 
dys a? Ox! oyé 
With w = wo cos ny/a, where wy is a function of x only, this 
becomes, with N = £Ef'/[{12(1 vy?) | 
d*w n*{? Cy 
+ »6 = a Wy =O (4) 
dx! 12(1 va‘ Fa‘ 


Equations similar to Eqs. (3) and (4) are derived in reference 4 
from Fliigge’s three differential equations by order of magnitude 
considerations 

Since Eq. (1) also obtains for a simply supported shell, with 
length / (Fig. 2), which is known to buckle in lobes with only one 
half wave in the axial direction, in that case w = wy cos 7x /l, 
so that from Eq. (1) C = x*N/I4 


as 


Hence Eq. (1) may be written 


d4w/dx* = (r'/l4)w (5) 


where / is the length of the simply supported shell that buckles 
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it the same critical pressure as the clamped shell. For the latter 


t 


case, owing to symmetry, the general solution of Eq. (5) is 


w = A cosh rx/l + Boos rx/l (6) 
Inserting this equation into the boundary conditions, vw = 0 
and dw/dx = 0, at x = xc (Fig. la), one obtains two homo- 


geneous linear equations, which for deflections different from 


zero, are only satisfied if the determinant of the system is zero, 


yielding 
tanh rc/] + tan rc/l = 0 7) 


Chis condition is satisfied if / = 1.38¢ = 0.665L => (2/3)L 
Hence, a fully clamped shell of length Z has the same critical 
external pressure as a simply supported shell of length (2/3)L 
As shown in reference 2, for steel shells and aluminum shells 
having L/a ratios larger than 1.5 and 2.25, respectively, that 
is, for shells having ratios of the half waves in the axial and 
circumferential direction that are sufficiently large (larger than 
about 4), excellent agreement exists with results obtained in 
reference 5. For shorter shells, reference 2 presents another 
solution. 

From reference 6 the von Mises equation for the critical pres 
sure p-, of simply supported shells can be reduced to a formula 
which may be written in the form 

O.92E(t/a)? 
il (8) 
(l/a) (a/t) /*? — 0.636 


Hence, for clamped shells of length LZ, with / = (2/3)L, 


1.38E(t/a)? 


— 0.954 


(9) 
(L/a) (a/t) 


A practically similar equation to the subject problem can be 
derived using the method employed by Professor P. Lardy in 
reference 7 in calculating the bending moments in clamped slabs 


Assume 
w = COS NELWyU (10) 
where the functions 


Cosh a,x COS &»,X 
u = _ 11) 


cosh an COS a» 


are the normal functions that automatically satisfy the boundary 
conditions, which for clamped shells require that 


tanh a,c + tan a,c = VO (12) 


For the symmetrical case the solutions of Eq. (12) are a, = 


2.365/c and a», = (4m 1)x/(4c) for m 2 2 As a matter of 
fact, for m = 1, Eqs. (11) and (12) are identical to Eqs. (6) and (7) 
The functions u,, are orthogonal Further 0*u,,/O0x' = a,,'u 
ind 071,,/Ox? = amt», where 
cosh a,,X COS AX 
ul + 13 
cosh a» COS Ami 
However, due to the orthogonality of the functions u,,, @,, can 
be developed into a series of u 
thm = Rimtty + Remtle + Ramtt3 4+ 14) 
From reference 7, for example, ky; = —O0.54984 


The solution converges rapidly. Using only the first terms of 
Eqs. (10) and (14) the writer found for a compressed plate that 


is clamped at the unloaded edges a buckling stress coefficient 


k = 7.04, which is in excellent agreement with the exact value 
k = 7.00. Also in the present case, due to the identity of Eqs 
(6) and (11), the latter for m = 1, a good convergence may be 


expected. Insertion of Eq. (10) into Eq. (2) and using only the 
first terms of Eqs. (10) and (14) yields, after solution for o, and 


t 
Pers 


FORUM 85. 


N(n58° + 12.3n°3° + 188.4n'84 + 386n282 + 985) + 31.4Et 7784 


B7a5(n*B° + 7.70n'B8* + 62.8n28? + 48.34 


15 
where 8 = c/a = L/2a and y = a/t. With relatively small half 
waves A = xa/n in the circumferential direction the ratio c/X = 
nc/ra = nB/x, and certainly n?8?, is large with respect to one, so 


that approximately 


Nun's? + 31.4Et 728% 
Poe = = 


B*a*n®3* 
n* 31.4 /fa\? ft 
= E (16) 
12(1 — »? n°B* \t a 


The condition Op,,/On = 0 gives n = 2.38 (¥/8?)'/* Inserting 
this into Eq. (15) and dividing the numerator by the denomi 


nator, leads to 


Per = [0.6918 y + 0.190 + 1.43/(By _ 
1.40/(By ‘*)? 4 E/(B*y*) (17) 
Taking only the first two terms of Eq. (17) gives 
1.382(t/a)? + 0.76(t/a)’*"(a/L) 
>. = E (18) 


(L/a) (a/t) 


which is practically identical to Eq. (9). 
For example, for Model 2 in Table 1 of reference 1, Eqs. (9), 


(17), and (18) give p., = 32.5, 32.4, and 32.3 psi, respectively, 
against p.- = 34.0 and 35.2 in columns C and D of Table 1 in 


reference | 
The writer wishes to thank the Bell Aircraft Corporation for 
permission to publish these results 
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A Supplementary Note on the Spanwise 
Loadings over Slender Airfoils 


Hsien K. Cheng 

Aerodynamicist, Aerodynamic 
Buffalo 5, N.Y 

July 30, 1954 


’ i SHIS NOTE concerns mainly the aerodynamic load distribution 
over a slender lifting surface. It is written with the aim to 


generalize a seemingly obscure conclusion arrived from load 


Section, Bell Aircraft Corporatior 


prediction for a slender delta wing with symmetrically linear and 
parabolic twists 
At a distance sufficiently far downstream of the ‘leading edge”’ 


of a flat plate of low aspect ratio—i.e Ww?— 1 (A.R.) < 1, 
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the spanwise lift distribution or the potential difference across 
the wing surface Ay, according to Jones’s slender-wing theory, is 


Ag = 2UaV a? — y? (1) 


where U is the free stream velocity, a = a(x) the local half-span, 
x and y the chordwise and the spanwise rectangular coordinates, 
and a@ the angle of attack. The cross-flow field of the lifting 
airfoil is thereby completely determined by the upwash distribu- 
tion along the local span, it follows that the spanwise lift dis 
tribution or Ag for a cambered slender airfoil of zero twist, i.e., 
0a/dy = 0, will also be elliptic and given by Eq. (1), in which 
both the angle of attack @ and the span 2a are functions of x 
This has been pointed out by Jones.! 

For twisted slender airfoils with zero camber, i.e., 0a/Ox = 0, 
a general feature in the load distribution somewhat analogous 
to the previous can be found. The cross-flow field in the neigh 
borhood of a slender plane lifting surface is? 


o¢ o¢ l "a VV a? — stu ds 

=i = g(f) = : (2) 
oy 02 Vat — f° J —a 3 
where ¢ = y + is, g(¢) the complex velocity in the cross-flow 


plane, i.e., ¢-plane, and wo the upwash distribution over the lifting 
surface at s = 0. For the airfoil surface with zero camber (as- 


suming a(y) remains finite), 

W=- Ualy), — 5, 7 < é, (3) 
the side-wash on the top surface of the wing, according to Eq 
(2), is 


te U a /a? — st a(s) dS 


=e ——— mae P oe ee (4) 
oy rV a? — J —2 s—-y, 
where the principal value of the integral is to be taken. The 
linearized pressure distribution or 0¢g/Ox on the top surface may 
be computed from Eq. (4) as follows: 


09 da ov 


Oyox ~ dx da 


Ua da ° as) ds 
nw dx as-y 


/ 
V a? — s? l 
i o\3 T / 
a? — ¥) Va? — y2 Va? 52 
which gives 
oO? y l 
¢ = A(x) - nr fe) a a<y<a 
Oyox (a? — y?) a? — y?)’/? 
(5) 
where 
Ua da “a ay) dy) 
A(x) . . 7 
° Ee J Age — y? 
(6) 
Ua da a yay) dy 
B(x) . 7 
x dx J_, Vat — y? 


Integrating Eq. (5) with respect to y along the top surface, one 


has 
o¢g l B(x) y 
" = A(x) - + ee + F(x), —a<y<a 
Ox Va? — y? a? Va? — y? 
(7) 


But in approaching the leading edge—i.e., | y| — @ the value of 
O¢g/Odx or the linearized pressure difference should either tend to 
infinity as a pure inverse square root or vanish. In other words, 
the constant of integration appeared on the right-hand side of 
the above equation, F(x), must be identically zero; otherwise, a 
logarithmic singularity in the upwash distribution will occur at 
the leading edges, which will certainly violate the boundary con- 
dition on the wing surface, namely, Eq. (3). Therefore, one 
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arrives, under condition Eq. (3) at 
o¢ 1 

A = 2A(x) - + ig : 
Ox V a? y? a* Va? — y? 


-a<y<a (8 


where A(x) and B(x) have been given by Eq. (6). 

Hence, within the scope of the slender-airfoil theory, the follow 
ing conclusions can be drawn: 

(I) If the plane lifting surface has no twist at the station 4 
i.e., a = a(x), the spanwise distribution of the lift (pUAg¢) at 
that station is the same as over a flat plate at angle of attack 
which is the elliptic distribution given by Eq. (1), as pointed out 
by Jones 

(II) If the surface has 
upwash is prescribed symmetrically with respect to the midspan, 


sero camber, 1.e., a = a(y), and if the 


i.e., a = +a(—y), the spanwise distribution in the pressure dif- 
ference (Ap ~ —pUAd¢g/Odx) will be found the same as over a 
flat plate at angle of attack, which is given by the first term of 
the right-hand member of Eq. (8). If, instead, the upwash 
distribution is prescribed antisymmetrically with respect to the 
midspan, i.e., a(y) = —a(—y), the spanwise distribution in the 
pressure difference will then be found the same as for a plate in 
rolling, which is given by the second term of the right-hand mem- 
ber of Eq. (8). 

It is of interest to note that the ‘‘zero cambered”’ airfoils men 
tioned above include not only the case of symmetrically linear 
twist as well as the case of differentially deflected triangular 
plate, but cases of control surface deflections in general, provided 
the ‘‘gaps’’ between the surfaces are of zero width and run across 
the wing’s plane from the leading edge to the trailing edge in the 
streamwise direction. 

In arriving at the above conclusions, it has also been assumed 
that the flow around the leading edge is ‘‘nonseparated”’ and that 
the planform is symmetrical, whose shape is, however, arbi 
trary. It is worth noticing that when the airfoil becomes not 
so-slender in the sense of reference 3, the first conclusion stated 
above can be shown still valid, while the second is no longer true 

In a private communication, H. R. Lawrence of Cornell Aero 
nautical Laboratory kindly informs the writer that he has also 
arrived at the same result of Eq. (8) in a different manner 
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A Semiempirical Method for Determining 
Delta-Wing Pressure Distributions in an 
Incompressible Flow 


William R. Laidlaw and P. T. Hsu 
Massachusetts Institute of Technology, Cambridge, Mass. 
August 4, 1954 


SYMBOLS 


b = root semichord 
root chord 


NW 


pressure correction function 
k = (whb/U) = reduced frequency of oscillation 
Ap(x, y, t) = lift per unit area 
s(x) = local semispan 
U = airspeed 
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x, y, o) downwash at camber surface 
chordwise coordinate (measured aft from vertex) 
spanwise coordinate 
= vertical coordinate 
x, y, bt) = vertical deflection of camber surface 
= air density 
perturbation potential in plane z 0 
circular frequency of oscillation 


DISCUSSION 


wr THE EVER-INCREASING popularity of the delta wing in 
the design of modern high-speed aircraft, there has come a 
demand for methods of analysis which are sufficiently rapid for 
use in the design office but at the same time may be relied upon 
to give aerodynamic information of reasonable accuracy 

[he delta-wing problem is by no means a new one. In the 
past, Jones! has considered the case of an extremely slender wing 
in a steady incompressible flow, whereas Garrick? and Voss* have 
extended this work to cover the unsteady case. Lawrence* and 
Lawrence and Gerber® have made worth-while improvements on 
the slender-wing theory to cover delta wings in a wide range of 
aspect ratio. The latter theory is, however, considerably more 

omplicated than the previous slender-wing analysis; and, at the 
same time, it is no longer possible to obtain the pressure distribu- 
tion. 

The object of the present discussion is, therefore, to propose a 
semiempirical approach to the problem, which retains the 
simplicity of the slender-wing theory and yet gives results com 
parable to those of Lawrence and Gerber. 

In general, the lift per unit area at any point on a thin lifting 
surface may be written: 

O¢ _ Og 


Ap(x, y, t) = 2p (x,y, + Ll 
aie ot 


Ox 


According to the theory of Garrick? and Voss,* the velocity po- 


tential g(x, y, t) is given by the relation: 
Ie 1 W(x, n, LV (x)? — 7? 
(x,y, t) = - dn 

Jy eV s(x)? — y? J -s =e 

(2) 
where the downwash W/,(x, y, ¢) is prescribed in terms of the wing 
motion Z,(x, y, 4) by the expression: 

0Z4 Za, 


(x, y, t) + l (x, y, t) (3 


W(x, y, t) = 
a(%, 3 Ox - 


Eq. (2) is valid for any arbitrary motion of the wing camber 
surface, and therefore includes, as a special case, the steady flow 
The only requirement attached to the above expression for 
g(x, y. t) is that the ratio of wing span to wing chord be very 
small In the past, there has been much conjecture concerning 
the consequences of this aspect ratio limitation when applying 
the theory to practical problems, and further studies were def 
initely indicated. 

Tests and calculations were conducted at M.I.T.® to determine 
pressure distributions on two rigid oscillating delta wings of aspect 
ratios 1.07 and 2.31, and the following conclusions have been 
stated: 

(1) The main fault of the slender-wing theory is that it fails to 
provide for the Kutta condition of zero pressure at the trailing 
edge. This causes an overestimation of the pressures on the rear 
portion of the wing, and gives rise to extremely large errors in lift 
and pitching moment even for the relatively slender aspect ratio 
1.07 wing. 

(2) The slender-wing theory predicts pressures near the vertex 
with good accuracy for both delta wings. 

(3) The shapes of the theoretical spanwise pressure distributions 
at fixed chordwise stations are in good qualitative agreement with 
the experimental results, although errors in ‘‘scale’’ exist which 
increase as the trailing edge is approached 
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(4) The pressure phase angles, in all cases, are in good agree 
ment with experiment over the complete plan form 

In view of the above statements, a possible means of improving 
the theory immediately presents itself: namely, the use of a 
multiplicative function which will permit a more realistic repre 
sentation of pressures near the trailing edge, and at the same 
1 


time leave the theoretical results unaltered in all other respects 


For this purpose, the function: 
F(x/c) = V1 — (x/c)? 


is suggested. It will be noted that F(x/c) is a real function of 
only the chordwise variable; hence, although it alters the pres 
sure amplitudes, it leaves the pressure-phase angles and the 
spanwise shape of the pressure distribution unchanged. Further 
more, conditions in the vicinity of the vertex remain identical to 
those of the unmodified theory; whereas at the trailing edge, it 
provides for the Kutta condition and gives the pressure distribu 
tion a slope singularity resembling that of two-dimensional thin 
airfoil theory 
If we then denote Ap, as the pressure at a point (x, y) as com 
puted from the slender-wing theory of Jones, Garrick, and Voss, 
the modified pressure Ap,, at the same point is given by the rela 
tion: 
App(x, y) = V1 — (x/c)*Ap, (x, y 4 


Eq. (4) was used to compute pressure distributions for the two 
delta wings described above. Comparisons were then made with 
the previous experimental and theoretical results. Excellent 
agreement with the experimental data was obtained for a wide 
range of reduced frequency [k = (wb/l’)] for both rigid body 
pitching and vertical translation motions. This was shown not 
only by the pressure distributions, but also by the chordwise lift 
distributions and the total lift and pitching moments 


theoretical and experimental comparisons may be found in refer 


Complete 


ence 6. 

In order to determine the validity of this pressure correction 
function when applied to the air forces which are exerted on a wing 
performing higher order elastic motions, a flutter calculation was 
carried out on a delta wing which had previously been tested at 
M.1.T 


vious slender-wing flutter calculations are thoroughly described 


The details of this wing, the flutter tests, and the pre 


in reference Briefly, the flutter model consisted of a delta 


‘ 
wing of 25.75° vertex angle (aspect ratio 0.915), fabricated of 
“rigid’’ spanwise segments but free to deform in the chordwise 
The model was restrained in the wind tunnel by a 
The ob 
served flutter speed and frequency respectively are given in 
Table 1 

The original flutter calculations employed an 


direction. 
simple support located at 32.15% of the root chord 


‘assumed mode” 
technique, which incorporated rigid body pitching about the sup 
port point, and the first two still air elastic modes, in the flutter 
mode. 

The calculations of reference 7 included a very accurate analy 
sis of the inertial and elastic parameters of the wing, and hence 
the large errors which appear in the flutter speed and frequency 
(Table 1) may be attributed primarily to the inadequacy of the 
unmodified slender-wing theory to provide a satisfactory rep 
resentation of the aerodynamic forces 


TABLE | 


Theoretical and Experimental Flutter Data 
0.915 Delta Wing 


Aspect Ratio 


Flutter Speed, Flutter Reduced 
Uy Frequency, w; Frequency, 
Case (M.P.H (Rad./Sec.) ke = (wr-b)/U; 

Experiment 39.1 14.6 0.557 
Unmodified slender 

wing theory 65.7 22.0 0.500 
Modified slender 

wing theory 33.5 11.5 0.512 
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The new set of flutter calculations, which have recently been 
completed, utilize the original structural analysis, but incorporate 
aerodynamic forces modified according to Eq. (4). Specifically, 
this amounted to a re-evaluation of the integrals which lead to 
the generalized forces in the flutter equations. The results of the 
revised analysis are presented in Table 1, and show a remarkable 
improvement over the previous theoretical speed and frequency 

The improvements resulting from the use of this empirical cor 
rection function, may be attributed to the following two factors 

First, it provides a means of satisfying the Kutta condition, 
thus giving a more realistic representation of the pressures in the 
vicinity of the trailing edge. 

Secondly, it has been observed that, according to the unmodi- 
fied theory, the pressure distribution on a delta wing is affected 
strongly by the slope and curvature of the camber surface 
Flutter results, therefore, become extremely sensitive to the shape 
of the wing, and since pressures near the trailing edge are greatly 
overestimated, this region assumes a high degree of importance. 
When these trailing edge pressures are adequately represented, 
the accuracy of a given flutter calculation becomes less dependent 
on the nature of the assumed modes. 

With the experimental evidence of references 6 and 7 asa basis, 
we may conclude that the semiempirical technique described 
above, for obtaining delta-wing pressure distributions, is satis- 
factory. It should be noted, however, that since the procedure is 
based essentially on a combination of experimental observations 
and physical reasoning, it should be used with discretion when 
applied to delta wings of large aspect ratio or to other than trian- 
gular wings. 

The method does, however, yield results which compare favor 
ably with the more exact Lawrence and Gerber theory, and at the 
same time it retains the simplicity of the slender-wing analysis 
These features should give it much appeal from the engineering 
point of view, particularly since it is valid for both the steady and 
unsteady cases regardless of the shape of the camber surface or the 
complexity of the oscillation mode. 
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Addendum—‘‘Characteristic Solution for 
Axially Symmetric Transonic Flow’”' 


Erik Mollo-Christensen 
Research Engineer, Department of Aeronautical Engineering, 
Massachusetts Institute of Technology, Cambridge, Mass. 


August 4, 1954 


T° A LETTER from Kelsey Walker, Jr., Project Engineer, 
Municipal University of Wichita, Kansas, he points out that 
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the solution of the transonic equation given by me! had _ pre 
viously been found by E. W. Graham in 1949.2 

At the same time, Prof. Leon Trilling, M.I.T., found that 
is included as a special case in a solution given by G. Guderley 
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Shock-Wave Yaw for 9 !4° Semiapex Cone 


Eugene D. Boyer 

Ballistics Research Laboratories, Aberdeen Proving Ground, 
Maryland 

August 5, 1954 


HE SEMIAPEX ANGLE of the shock wave and the ratio of the 
yj pentians yaw to the yaw of a projectile for small yaw is 
theoretically predicted by A. H. Stone.' Wind-tunnel verifica 
tion of this theory is described in reference 6. An experimental 
verification of the theory has also been obtained from cones and 
cone cylinders’ fired in the Free Flight Aerodynamics Range of 
the Ballistic Research Laboratories. (See Figs. 1 and 2.) 

The theoretical values as predicted by Stone are obtained di 
rectly from Kopal’s tables.” In order to determine thes« 
values for a cone, semiapex angle 6, of 9! »°, we must first deter 
mine the surface velocity component l’,, in terms of the Mach 
Number J/. Knowing l., we can then obtain the shock-to 





Fic. 1. Shadowgraph pure cone 








Fic. 2. Yawing cone. 
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iG. 3. Ratio of shock-wave yaw to yaw of projectile vs. Mach 
Number 


cone-yaw ratio 6/e, by two-way linear interpolation. Fig. 3 
shows the predicted behavior of the ratio 6/e for various Mach 
Numbers. 

For the empirical values a number of rounds with different 
Mach Numbers were selected for a distribution of points on the 
5/evs. Mcurve. Since the change in .1/ is less than 0.16 through 
the range, each round is considered to have a constant Mach 
Number. The yaws were measured only in one plane, vertical or 
horizontal. Only those stations were selected for which the yaw 
in one plane was very small, less than ! , and the yaw in the 
measured plane was fairly large, at least six times as great 

The projected angular positions 6* and e* were measured on a 
grid. From these measured angles the true yaw is obtained.‘ 

Since we are considering a body of revolution it is reasonable to 
assume that for zero yaw the yaw of the shock wave is also zero 
rhe ratio of the two yaws can then be determined by fitting the 
function 6 = be, by least squares, to the yaw angle data. These 
values along with the theoretica! values are seen in Fig. 3 plotted 
against Mach Number. The agreement is good with the dif- 
ferences of the experimental and theoretical values being less than 


the statistical error of the experimental values. 
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Effect of Spin in Aerodynamic Properties of 
Bodies of Revolution 


L. E. Schmidt and C. H. Murphy 

Ballistics Research Laboratories, Aberdeen Proving Ground, 
Maryland 

August 6, 1954 


EHNERT AND HASTINGS" ? have measured the effect of spin on 
L base pressure of cone-cylinders at a Mach Number of 2.86 
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and for Reynolds Numbers from 2.5 XK 10® to 4.5 & 10° They 


found that the character of the boundary layer at the base was an 
important factor. An increase in spin would cause either a con 
siderable decrease in base pressure if the flow was laminar or a 
slight increase if it was turbulent. Recent firings* in one of 
BRL’s spark ranges indicate a different situation at R, = 3.0 
108 and M = 1.3 

The primary purpose of these firings was to study the effect of 
spin on the aerodynamic coefficients of spin-stabilized missiles 
rhe model was a 5-caliber long body of revolution with 2-caliber 
secant ogival nose and a cylindrical body. The radius of the 
ogive was 8.5 calibers Measurements were made of total 
drag, Cp, normal force, Cy» overturning moment, Cu damping 
and spin decel 


moment, Cy. 4 Cu, Magnus moment, Cy 
pa 


erating moment, Ci», coeflicients at different values of spin 
pd/2V, between 0.079 and 0.323 (d is model diameter Since 
Cp is quite sensitive to the magnitude of yaw, only those rounds 
with less than 3° of yaw were considered for drag determination 
Only the total drag and the Magnus moment coefficients showed 
a measurable variation with spin. 

In reference 5 it is shown that for configurations possessing a 
plane of mirror symmetry the coefficients are even functions of 
spin. For this reason Cp and Cup, were plotted against 
(pd/2V)*?. A linear relationship in spin squared was assumed 
and fitted by least squares to the Cp and Cup, data (see Figs. 1 


and 2). The slopes, Cp, Cy, , and intercepts Cp(0), Cy, (0) 
pa pa 
were computed to be: 
C)/(0) = 0.415 + 0.002 Cup (0 = 0.64 + 0.05 
Cp = 020 + 0.02 Cy = —5.1 + 0.8 
pa 


Although CMe is poorly determined, the trend is statistically 
significant. 

In all firings the boundary layer at the base was turbulent and 
transition occurred between one and two calibers forward of the 
base. An upper bound to the contribution of the skin friction 
variation with spin to Cp was estimated* to be 0.05. The dif- 
ference of @.15 must then be caused by a decrease in base pressure 
for increasing spin. For comparison purposes Cp fcr base pres 
sure variation was computed from the Lehnert and Hastings 





data For laminar flow at R. = 3.0 X 10%, a = ().44 and for 
Cy 
45 + 


40 + / 
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Fic. 1. Zero-yaw drag coefficient vs. spin squared. 
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Fic. 2. Magnus moment coefficient vs. spin squared 
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turbulent flow of about two calibers extent, forward of the base, 
at R, = 3.8 X 108, Cp = —0.14. 
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Ballistic 


Aerodynamic Symmetry of Projectiles, 


Tip or Leading-Edge Temperatures on Pointed 
Heat Conducting Bodies at High Supersonic 
Speeds* 


M. Richard Denison 
North American Aviation, Inc., Downey, Calif 
August 8, 1954 

SYMBOLS 


A = cross-sectional area at distance x 

AL cross-sectional area at distance L 

I 32hALSLL/kAL dimensionless parameter for cone 
G 16h_pS_LL/kAL, dimensionless parameter for wedge 
h heat-transfer coefficient at distance x 

hy heat-transfer coefficient at distance L 

Ii modified Bessel function of first kind 

k thermal conductivity 

Ky modified Bessel function of second kind 

L = distance from tip or leading edge to base 

Q = heat absorbed by sink 

5S = surface area from tip to distance x 

SL = surface area from tip to distance L 

7 = body temperature at distance x 

Tiw = insulated wall or recovery temperature 

To = tip or leading edge temperature 

u = (phLSLL/NRA t)'/2 (x/L)%, dimensionless variable 
x = distance from tip or leading edge 

i) Tiw — T 


DISCUSSION 


F A BODY WERE TO move through the atmosphere at a fast 
I enough speed, recovery temperatures high enough to melt 
any known material could be reached. It is known that over 
most portions of the body surface temperatures could be re- 
duced by cooling the boundary layer, but there seems to be 
widespread belief that the temperature at the leading edge or tip 
of a structure must necessarily be at the recovery temperature 
For example, see references 1 and 2. This belief probably arises 
because in these regions the boundary layer is so thin that the 
local heat-transfer coefficient becomes very large 

At the suggestion of Dr. E. R. Van Driest of North American 
Aviation’s Aerophysics Department, the following study was 
carried out to determine the validity of the above concept. The 
results show that the concept is overconservative even with a 
heat transfer coefficient that approaches infinity at the leading 
edge. ‘ 

Consider steady-state heat conduction in a solid pointed body 
at zero angle of attack with a heat sink at the base. Included 
among such bodies are cones and wedges which are important 
shapes for aircraft design. It will be shown that for both the cone 
and wedge with a heat sink, the tip or leading edge temperature is 
less than the recovery temperature. If it is true in these simple 
cases, then it must also be true in more complex cases that the tip 


Fox of the Aerophysics Depart 


* The author wishes to thank Dr. J. L 
ment for his advice in this study 
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SCIENCES 


or leading edge temperature depends on the particular heat con 
duction conditions which exist 

Assume that the body is sufficiently thin so that the internal 
temperature gradients normal to the surface are negligible. Thus 
the body temperature is a function of the distance from the tip 
only. This will be designated the x direction. Assume, also, that 
the cross section of the body changes with a power of x, and neg 
lect radiation from the surface 

A heat balance for an element gives, neglecting higher order 


terms, 
d?T dA di ‘is. . 
kA +k + h (Tt <a 2) ae | 
dx? dx dx dx 
Let@= 7;, — T. Then 
d*6 dA dé dS 
; 4 — f @=0 2 


+ 1 
dx? dx dx dx 


Assume that A, S, and #4 are power functions of x and express 


Tt 


them in terms of their values at the base of the body, x = L 
For convenience in the final result introduce constants A and 


in the powers of x such that 


A = A,(x/L)'-?™ 2a 
Hence, 
1 dA 
= (1 — 2dy)x™! (2b 
A dx 
Also let 
dS PhiSi ( =n() . 1 
h = ) (2 
dx £ s 


where p is another constant. When Eqs. (2b) and (2c) are sub 
stituted, Eq. (2) becomes 


d*6 de ( phiS1 


+ (1 — 2dy)x7! VAY = 0 3 
kA, Lr} 


dx? dx 


(Fete . ( A 
“= 

kA yd? : # 
In terms of u, Eq. (3) can be written as 


d*6 (1 — 2v) do 
du2 u du 
rhe solution of Eq. (4) is 
6 = u’[q I,(u) + @K,(u)] (5 
where J, and K, are modified Bessel functions of the first and 


second kind respectively. 


SOLUTION FOR CONE 


The cross-sectional area of a cone is proportional to x?. Hence 

A = Ax(s/1)* 

Comparison of this with Eq. (2a) gives 
2=1— 2dr (6) 


The surface area is also proportional to x? so that 
S = S1(x/L)? 


and 


dS/dx = 28,/L(x/L) 


For compressible laminar flow with constant-wall temperature 


and static pressure, Crocco‘ has shown that 

+ Gardner? has solved this equation with hk constant, A and S power 
functions of x He also used different boundary conditions than will be used 
here. 
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con 


h = hi(x/L . 100 -————__ SS SS 
ernal Although the temperature gradient will be large and the Prandtl 50 |- | | 
rhus boundary layer equations are not applicable in the vicinity of the 
e tip cone tip, this relationship will be used for the entire region. Then 
that dS 2S thy x /2 
neg h— = 
dx L L 10 a 
rder Comparison of this with Eq. (2c) gives p = 2 and 5 
1/2 = 2\(1 — ») - 1 (7 | 
— a ee Tin To 
l Simultaneous solution of Eq. (6) and Eq. (7) gives ( Iw o 
QL/KA 
y=-2, A#=1/4 70M, 
| + — 
rherefore, the solution for a cone is by Eq. (5) 
2 @= u 2{col-(u) + oK—2(u)] (8) 05 


where 

A B2ALSLL t f/x \'/6 
ee RA I 

di EL af Ol 1 4 > 4 4 

Let F = 32h,S_L/kAt. Then u = F/*%x/L)/* Now the 











Ja boundary conditions must be satisfied. These are: 0.05 | | ] T 
atx = Ooru = 0, @ is finite 
'b atx = Loru = F*, kAt (). =o 
dx 0.01 
0 02 04 O06 O08 1.0 1.2 1.4 


where Q is the amount of heat absorbed by the heat sink at the 
base of the cone 








( Since K-.(0) = o@ the first boundary condition gives c = O 
In order to find ¢,, d0/dx must be obtained a , 
, Fic. 1 lip or leading-edge temperatures for cones and wedges 
From Eq. (8) with a heat sink at the base 
dé dI , ‘i — 
= 2u3J-2 + u~? be seen that the temperature gradient d7/dx goes to infinity at 
5 os 2 
du du the origin as x7’? 
= ¢u~3[uJ_3]* 
H WEDGE SOLUTION 
ence 
P k For a wedge of infinite span a similar procedure gives 
dé v2 ; 
i “Gg u>*]_3(u) (9 »=0, \=1/4, p=! 
ad 4 
r : i From Eq. (5) the solution for the wedge 1 
Atx = Loru = F/? Eq. (9) and the second boundary condition . —* 
give 0 = O,'To(u) + Co’ Ko(u) (14 
1OL Fi with 
ie i/o 10 1 
kA, I-3( F7/*) (=) (&) ‘ 
“= 
Finally kA1 L 
1i0L F? Let G = 16h,Si_L/kAx. Then with the same boundary condi 
ali ie ‘ kA, I-3( F*) u~*T-2(#) I tions as for the cone 
. ' ‘=(0 
From the properties of Bessel functions 7 
(1/2u)? + 2™ ’ 1OL I 
I_,u) = Iu) = +? a ; ” ~ kl G*I,(G/*) 
mao MT(v + m+ 1 a Sa 
‘ - . Since /,(0) = 
Hence from Eq. (8) ince [ 
— a oe. > 
(1 /94)2m 0o= liw = : —— (lo 
ee oe (1/2u) (12 kA} i 1G 2) 
| ro 
| } miT(3 + m) = ' . 
fea Therefore, the leading edge temperature for the wedge is also less 
| Atx=O0Oars =0 than 7;,. if Q is positive. It again turns out that d7°/dx goes to in 
T.. - T, = a/4X1/2 finity at the leading edge as x72, 
is The tip or leading edge temperatures for cones and wedges with 
so that a heat sink at the base are presented graphically in dimensionless 
OL Fi/: form in Fig. 1 
7, = Ty.= F "ye (13 
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On Panel Flutter in Supersonic Flow 


W. C. Randels 
Lockheed Aircraft Corporation, Burbank, Calif 
August 9, 1954 


HE WORK IN reference 1 analyzes the possiblity of flutter of 
gpm with multiple simple supports. The authors seem dis- 
turbed by the results which do not offer guide lines for the design 
of skins for supersonic vehicles. It is suggested that this may be 
due to the inapplicability of the model used. In reference 2 the 
more plausible model of a panel clamped at both ends is analyzed 


giving a criterion for design. 
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Throat Station Concept in Transonic Flow 
over Two-Dimensional Wedges 


Robert Wesley Truitt 

Professor of Aeronautical Engineering, Virginia Polytechnic 
Institute, Blacksburg, Va. 

August 11, 1954 


— FLOW OVER THE front surface of a wedge at zero angle of 
attack is particularly adaptable to the concept of the con- 
vergent portion of a de Laval nozzle. The existence of a stagna- 
tion point at the leading edge of the wedge serves as the ‘‘reser- 
The favorable pressure gradient over the front surface 


voir.” 
accelerates the flow from zero speed at the leading edge ‘‘reser- 
voir’ to a maximum speed at the wedge shoulder (throat station). 
Under the assumption of isentropic flow along the front wedge 
surface, the area-pressure relation may be written as:! 


2 (y + 1)/(4 1) 
any 3 er, 
= 2 (2) (2)" “4 (1) 
Pi pi 


where A, is the local cross-sectional area of the stream-tube cor- 


Ar 


responding to the local static pressure pz; A,r is the cross-sec- 
tional area of the stream-tube at the minimum throat (wedge 
shoulder); pp is the pressure in the ‘‘reservoir’’ (stagnation 
point ). 

Experimental pressure distributions for a 10 per cent thick 
diamond airfoil? in the form of p,/p; have been used in Eq. (1) and 
the results plotted in Fig. 1. It is seen that in the subcritical flow 
range—i.e., where 4,/A7 > 1.0 at x/c = 1.0, the area ratio at 
each x/c-station decreases with increasing free-stream Mach 


Number. However, at the critical Mach Number—i.e., when 


A,/Ar = A,/A* = 1 at x/c = 1.0, any further increase in free- 
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stream Mach Number results in practically no change in the 
area ratio at a given x/c-station. In a fixed-wall two-dimensional 
converging channel, if the Mach Number reaches 1.0 in the 
throat, the local Mach Number at all other sections becomes fixed 
at a definite subsonic value. This is evident from Eq. (1) since, for 
isentropic flow, 


fo (ren ee | (2 
Po 2 


where j/;, is the local Mach Number his throat station con 
cept was used previously by Tsien and Fejer.* The advantage 
of the present application is that the throat station is known to 
occur near the wedge shoulder, whereas, in reference 3, its location 
necessitated a trial-and-error method for a conventional subsonic 


airfoil 


CONSTANT LocAL MACH NUMBER RANGE 


From the above, the beginning of the constant local Mach 
Number distribution can be found approximately by a com 
pressibility correction method. For example, applying the 
Karman-Tsien method to low-speed experimental data for a wedge 
of semiangle 7.5°, the critical Mach Number corresponding to 
the peak negative pressure near the shoulder was found to be W,, 
= O.S15 
of the constant local Mach Number range is approximately 


Thus, in transonic similarity notation, the lower limit 


Ey = * —().72 (3 

Due to “symmetry” about free-stream Mach Number one, the 
approximate upper limit of the constant local Mach Number range 
can be found from the exact normal shock relation:! 


1 + [(y — 1)/2] Af? 


M,? = 
yl? — [(y — 1)/2 


where J/; and \/. are the Mach Numbers before and after the 
normal shock, respectively. The upper limit Mach Number, 
Moy.p. = Ah, is found from Ep. (4) by substituting Me= Moris 
= 0.815. This gives My, = A, = 1,246, which corresponds 
to the upper limit reduced free-stream Mach Numbers: 
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It is assumed that the detached shock is approximately a nor- 
mal one. It is known, of course, that for & ,-;, = 1.2 it is some- 
what curved and, as such, the flow behind the shock is slightly 
nonuniform. However, neglecting shock curvature, the local 
Mach Number distribution is approximately constant in the 


range —0.72 < & < 1.2 


PRESSURE DRAG 


The pressure distribution and hence the pressure drag can be 
found by applying, say, the Karman-Tsien method to low-speed, 
(M, =~ 0) experimental pressure data for a wedge up to and in- 


cluding the critical Mach Number—i.e., £o,,;,, = —0.72. Now 
the local Mach Number distribution, corresponding to the 
Karman-Tsien solution at é,,;,, = —0.72, will remain virtually 


constant in the range —0.72 < &< 1.2. Using a method similar 
to that used first by Maccoll and Codd‘ and later by Griffith,’ it is 
a simple matter to find the pressure distribution and pressure drag 
through the transonic range —0.72 < & < 1.2. 

As described, the Karman-Tsien method was applied to low- 
speed experimental data for a 7.5° semiangle wedge with straight 
afterbody. The results are plotted in transonic similarity form 
along with previous experimental and theoretical data in Fig. 2 
The agreement of the present method with other findings is 


evident. 
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Effect of Helicopter Blade Stall on Blade 
Flapping 


K. H. Hohenemser 
McDonnell Aircraft Corporation, St. Louis, Mo 
August 9, 1954 


N REFERENCE | an equation is derived giving the flapping angle 
I of a helicopter blade between azimuth wr and 27” under the 
assumption that the blade is completely stalled and that its lift 
coefficient is constant in this azimuth angle range. Assuming 
further zero flapping velocity at azimuth 7 the flapping angle at 
2x becomes a function of the flapping angle at 7. For a typical 
numerical example the flapping angle at 27 (downwind) is down- 
ward when the up flapping angle at x (upwind) is more than 10° 

The same concept as to region of blade stall had already been 
used in reference 2 to derive a complete solution of the problem 
of blade flapping in conditions of extreme blade stall which is free 
of specific assumptions as to flapping angle and flapping velocity 
at the boundary of the stalled region. Reference 2 theory was 
checked against results of wind tunnel tests with a 10-ft. inter- 
meshing rotor model, conducted in 1939, and was found satisfac 
tory in the applicable region. These tests included nose-up condi 
tions of extreme blade stall, leading to very large flapping angles 
The droop stop of —3° did not prevent the blade tips from assum- 
ing downwind positions corresponding to down flapping angles of 
much more than —3° 

These tests, in agreement with reference 2 theory, have shown 
that large clearances between blades and fuselage or boom are 
necessary (in spite of limiting droop stops) in order to avoid 
blade-fuselage or blade-boom interference in conditions of high 


nose-up altitudes producing severe blade stall 
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mittee at ¢Brown University, Providence 12, R.I.: Prof. D. C. Drucker, Secretary; Prof. E. H 
or Prof. W. Prager, Chairman. 


Lee, Treasurer; 

















—— a wast + A rer aa gt a onl ang ace Se Sonim lin mcrae tine 0 Ys Sie ih Ra habit hee PTs 





and in 
use of 


quent 

makes 
partic 
service 


name 
should 
which 


beginn 
as an? 
major 
stitute 
referer 
cially 

papers 
the mz 
contail 
paper. 


Some fF 
nomen 
Please 





SUGGESTIONS for CONTRIBUTORS 
to the Publications of the 
INSTITUTE of thee AERONAUTICAL SCIENCES 


The Institute of the Aeronautical Sciences invites both members and nonmembers from any 
country to submit papers for publication in the JOURNAL OF THE AERONAUTICAL SCIENCES and the 
AERONAUTICAL ENGINEERING REviEW. The Institute, following the practice of other societies, does 


not pay for contributions. 


The following directions for the preparation of papers, if followed by authors, will save corre- 
spondence, avoid the return of papers for changes, minimize the work of preparation for the printer, 
and save the expense due to the charges made for ‘“‘author’s corrections.” 





Manuscripts: The original typewritten copy of the paper is 
desired, double or triple spaced on one side of white paper sheets, 
consecutively numbered. There should be wide margins to 
allow for the marking of directions to the printer. Correcting, 
changing, or adding to papers after they are in type is costly. It 
is, therefore, imperative that papers submitted be in final form. 
Typographical errors may be corrected on proofs, but if authors 
wish to add material, they may do so at their own expense. In 
mailing, drawings may be rolled but manuscripts should be sent 
flat. Send by first-class mail (register if you wish for your own 
protection) to the Editorial Office, Institute of the Aeronautical 
Sciences, 2 East 64th Street, New York 21, N.Y. All manu- 
scripts will be examined by the Editorial Committee and by the 
Editor. Authors will be advised as promptly as possible whether 
the paper is acceptable for publication. 


TitLtEs: The title of the paper should be brief. The name 
and initials of the author should be written as he prefers. The 
use of the full name of an author is advocated because of the fre- 
quent duplication of initials and surnames which sometimes 
makes it difficult to establish the identity of the author. This is 
particularly important for large annual indexing and abstracting 
services. All titles amd degrees or honors are omitted. The 
name of the organization with which the author is associated 
should be placed after his name on a separate line. The date on 
which the paper is received will be inserted by the Editor. The 
author’s title should be indicated in a footnote. 


SUMMARIES OR ABSTRACTS: An abstract to be printed at the 
beginning should accompany each article. It should be adequate 
as an index and asa summary. It should contain a statement of 
major conclusions reached, since summaries in many cases con- 
stitute the only source of information used in compiling scientific 
reference indexes. Abstracts printed in other journals, espe- 
cially foreign, in most cases, consist of summaries from printed 
papers. The summary should explain as adequately as possible 
the major conclusions to a nonspecialist in the subject and should 
contain from 100 to 300 words, depending on the length of the 
paper. 


SuBHEADINGS: Subheadings should be inserted by the author 
at frequent intervals. The work of editorial preparation will be 
simplified by the author’s provision of many subheadings. 


Matrer Usvatty Devetep: Photographs or illustrations 
of little technical interest and not showing advances in general 
practice. Too detailed tabular matter (general results of such 
tables may be included in the text). Lengthy descriptions of 
materials or processes or of preliminary experiments or theories 
that preceded final results; salient features only are of interest. 


REFERENCES AND Footnotes: References should appear as 
footnotes only, numbered consecutively and grouped together at 
the end of the manuscript. The arrangement should be as fol- 
lows. For books: ! Durand, W. F., Aerodynamic Theory, 1st 
Ed., Vol. 1, p. 23; Julius Springer, Berlin, 1934. For maga- 
zines: ! England, C. R., Crawford, A. B., and Mumford, W. W., 
Some Results of a Study of Ultra-Short-Wave Transmission Phe- 
nomenon, Proc. I.R.E., Vol. 20, No. 12, pp. 481-482, March, 1933. 
Please give author, title, edition, volume, page, publisher, and 


date of publication as indicated. Omission of one required fact 
causes much extra editorial work and possible inaccuracies. 


ILLUSTRATIONS: Illustrations should accompany manuscripts, 
and each should always be referred to in the text by number. 
Drawings or graphs should not be larger than 12 X 16 inches, 
and must be made with jet black India ink on white paper or 
tracing cloth, the latter being preferred. Do not use typewriter 
for lettering. The smallest lettering on 8 X 10-inch figures should 
be no less than '!/, inch high. Cross-section paper (white with 
black lines) may be used but it should not have more than 4 lines 
per inch. If finer ruled paper is used, the major division lines 
should be drawn in with black ink, omitting the finer divisions. 
In the case of finely ruled paper, only blue-lined paper can be 
accepted. Tracing paper and blueprints are not acceptable. 
Lettering and all markings must be large enough to be readable 
after reduction. Mail rolled or flat; never fold. Drawings that 
cannot be reproduced (including pencil drawings) will be returned 
to the author for redrawing thus delaying publication of the 
paper. Photographs should be distinct and show clear black 
and white contrasts. They must be on glossy white paper. 
Avoid round and oval photographs. 


CAPTIONS AND LEGENDS: Legends or captions must accompany 
each drawing or photograph submitted. If written on the draw- 
ing or photograph, they should be placed below and well outside 
the part to be reproduced. Each table should have a caption 
such as Tabie 1, Table 2, Table 3, etc. Captions should be com- 
plete in themselves so as to make the data intelligible to the 
reader without reference to the text. A duplicate list of captions 
for figures should be included as the last page of the manuscript. 
Use “Fig. 1’”’ (not Figure 1), “Figs. 3 and 4,” etc., in both the text 
and the numbering of illustrations. In the text, ‘‘Eq. (1)” or 
“Eqs. (1) and (2)” is used, not “Equation (1).”’ In captions 
and legends, except for ‘“‘Fig.’’ and ‘‘Eq.,”” and in table headings 
write all words in full; do not abbreviate. 


MATHEMATICAL WorK: Only the simplest formulas should be 
typewritten; all others should be carefully written in pen and 
ink, the writing to be large enough so that ample room is provided 
to mark mathematical matter for the printer. A considerable 
spacing for marking should be allowed above and below all equa- 
tions. All complicated equations should be repeated on separate 
sheets with plenty of space left for marking. The solidus should 
be used for simple fractions appearing within the text. Make 
all expressions clear to the typesetter. Greek letters used in 
formulas should be clearly designated by name on the margin of 
the manuscript. All symbols should be clearly written and care- 
fully checked. The difference between capital and lower-case 
letters should be clearly distinguished and care taken to avoid 
confusion between zero (0) and the letter (0), between the 
numeral (one) and the letter (ell) and the prime (’), between 
alpha and a, kappa and k, u and mu, v and nu, n and eta. All 
subscripts and exponents should be clearly marked, and dots and 
bars over letters or mathematical expressions should be avoided. 
Avoid complicated exponents and subscripts. When it is neces- 
sary to repeat a complicated expression, it should be represented 
by some convenient symbol. 

NOMENCLATURE AND ABBREVIATIONS: Standard abbreviations 


should be used, and it should be noted that most abbreviations 
are lower case such as m.p.h., b.m.e_p., i.h.p., b-hp., hp., . . .ete. 
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